#### 1. INTRODUCTION

Deposition of ultrathin (island) layers of metal on the dielectric surface significantly changes its properties, incl. and optical. Theoretically, this can be related to the behavior of the electron gas of nanoparticles in external electric and magnetic fields. Thus, in a homogeneous electromagnetic field of the optical frequency, plasmon oscillations of the electron density are excited in the conducting nanoparticles, as a result of which the near field near the surface of the particles changes [1-3]. It is obvious that in a system with a large number of particles, it is necessary to take into account their mutual influence on each other [4-5]. In the case of a regular spatial arrangement of particles, this calculation is quite simple.

Nanolattices have a number of properties that make them promising for use as elements of a new generation of miniature biosensors and compact optical radiation control elements. Anomalous transmission [6], anomalous opacity (non-transparency) [7], artificial optical and magneto-optical [8] activity, as well as the possibility of plasmon focusing [9-13] are described in the literature. The development of compact optoelectronic sensors is closely related to the need to obtain reproducible resonance characteristics in the interesting part of the spectrum.

In the first part of the paper, the field structure was studied in a quasistatic approximation and the probabilities of induced transitions of molecules located in a regular lattice of long parallel conducting rods of nanometer radius were calculated.

#### 2. REGULAR GRATINGS COMPOSED FROM LONG CONDUCTIVE NANORODS (QUASISTATIC APPROXIMATION)

We orient the z axis of the Cartesian coordinate system along the axis of the lattice nanorods, and choose the direction of the external electric field **E
** along the x axis perpendicular to z (Fig. 1), then the resulting field

We introduce the vector
*dipole moment* of a unit length of a single nanorod of radius R by the relation

Here

– is polarizability of unit length of nanorod, the dielectric constant and the plasma frequency of the metal, respectively.

Potential of the summary near field outside the nanorod is given by expression

and the field strength (Fig. 1) is given by expression

In the case of regular spatial arrangement of parallel nanorods in the form of one-dimensional chain or two-dimensional grating, the effective dipole moment of unit length of each rod can be represented as

(1)
and the dyadic quasistatic Green's function

(2)
*i* and *j*, **e
** and

**e**2
are the elementary vectors of two-dimensional periodic nanorods lattice.

The field potential in such lattice is given by

(3)

It can be seen from (1) and (2) that more distant conductors contribute less to both the effective polarizability of each rod and the total field potential at a given point. Therefore, in practical calculations it is sufficient to take into account the influence of several nearest conductors.

For quasi-one-dimensional chain of parallel rods, one of the elementary vectors becomes infinite and summation over the corresponding index is not carried out. Numerical calculations show that the characteristics of the field of parallel chain of nanorods depend on the orientation of the vector **E
**. Therefore, the local distortions of the field have the largest amplitude with orientation

**E**0
across the axis of chain (Fig.2).

In quasi-two-dimensional geometry, calculations were performed for gratings with square

Because of the known symmetry of these gratings, the resulting field in them depends periodically on vector **E
** orientation. For the case of rectangular package, this dependence is repeated through an angle

**E**0
is oriented along vectors

**e**1
and

**e**2
of structure.

Simulation also shows that as the size of the system increases, the amplitude values of the local field distortion decrease, and the spatial dimensions of such regions increase. In this case, along the direction of the vector **E
**, co-directed with the lattice vector

**e**1
, the energy density of the field is concentrated near the nanorods, and in the transverse direction the maximum of the field is located between them.

It is shown that in calculating the probability of induced transitions for a molecule in nanorods lattice, it is necessary to take into account the orientation of its transition dipole moment
**E
** and relative to the principal axes of the quasi-two-dimensional structure. The case of a coincident orientation of the vectors

#### 3. POLARIZATION PROPERTIES OF NANOLATTICES COMPOSED FROM SPHERICAL LAYER COMPOSITES WITH EXCITONOGENEOUS SHELL AND MAGNETIZED METAL CORE

#### Polarizability of a particle in a planar lattice in the approximation of the first two coordination spheres

Let us now consider the aggregate of nanoparticles forming an infinite planar lattice of hexagonal symmetry. We will calculate the interaction of polarized nanoparticles in the approximation of the first two coordination spheres, as shown in Figure 6. The effect of particles in other coordination spheres is not taken into account, because their contribution to the polarizability of the particle in question will be negligibly small.

Since such a system is symmetric, the induced polarizability will be the same on all particles of the infinite lattice. Then taking into account the six nearest neighbors of the first coordination sphere and the twelve dipoles of the second coordination sphere in the lattice gives the following expression for the dipole moment vector of the isolated nanoparticle

Here

Otherwise, by regrouping the terms in (4), we obtain

Then the effective polarizability tensor of a spherical nanocomposite in a lattice in the approximation of the two first coordination spheres can be written in the form

For the dipole polarizability tensor
**B**, the following expression was obtained previously [3]

Here

The components of the tensor

Below, in Figures 7 and 8, the spectra of imaginary and real parts of the polarizabilities of a separate layered nanoparticle and a layered nanoparticle that forms part of a planar lattice are presented. It is seen from the graphs that for the transverse polarizability of the nanoparticle in the lattice, the resonances will shift to the low-frequency region, and for the longitudinal polarization to the high-frequency region of the spectrum relative to the polarizability resonances of the isolated nanocomposite. In this case, the amplitudes of the exciton resonances of the transverse polarizability

In a magnetic field, the amplitudes of all three resonances of the polarizability of a layered nanoparticle in a planar lattice decrease, and plasmon resonance, in addition, splits into two spectral components that expand with increasing magnetic induction *B*, just as in the case of a separate nanocomposite, as demonstrated in Fig. 9 and 10.

Thus, for a layered nanoparticle with an excitonogenic shell and a metal core located in a planar lattice of the same particles, the exciton and plasmon resonances are displaced, and the effect of the external magnetic field is analogous to the case of a solitary layered nanoparticle.

#### Comparison of the dipole polarizabilities of a separate layered nanocomposite with its lattice analog

Next, we compare the polarizability spectra of a separate layered nanocomposite

As can be seen from Figures 11 and 12, resonances of the transverse polarizability

For resonances of the longitudinal polarizability of particles in a plane

Thus, a comparison of the polarizability spectra of a separate layered nanocomposite and nanocomposites included in a planar and bulk lattice shows that the symmetry and dimension of such lattices cause a displacement of the resonances, which is more pronounced in the case of a planar lattice.

#### 4. FDTD MODELING OF ELECTROMAGNETIC FIELD IN THE SQUARE LATTICE FROM NANOCYLINDERS

An alternative way to calculate the electromagnetic field can be a direct numerical solution of Maxwell's equations by the FDTD method [14]. One of the advantages of this method is the automatic accounting of the delay effect.

In this section, the spatial distribution of the electromagnetic field in periodic nanostructures composed of metal nanorods with activated plasmon modes was calculated by the FDTD method. Simulation was performed using the MEEP software package [14].

For a numerical experiment, a two-dimensional countable domain was constructed. In this domain was a fragment of a square grating of nanorods. The scheme of the numerical experiment is shown in Figure 15. The axes of the nanorods were perpendicular to the counting domain. The radius of the nanocylinder was r = 25 nm, and its length is equal to infinity. The distance between the axes of two adjacent nanocylinders in the X direction was 3r, in the Y direction was 4r.

In the simulation, two types of nanocylinders were considered: from silver (Ag) and from silicon dioxide (SiO

Calculations were made for two different external media with dielectric permittivity

A planar linearly polarized monochromatic wave with a wavelength λ = 390 nm (maximum of plasmon resonance in a silver nanocylinder in vacuum) was modeled as the initiating field. The plane of polarization of the electromagnetic wave coincides with the plane of the countable region (Fig. 15).

Figure 16 shows two-dimensional maps of the electric field average intensity spatial distribution in the interaction of an electromagnetic wave with nanocylinders in a square grating.

It can be seen from Fig. 16a that the field structure in a medium with metallic nanorods is highly inhomogeneous. We see a lot of "hot spots", where the field is locally amplified.

On the other hand, in Figure 16b for comparison, the distribution of the field is presented in the case when the grating is composed of nanocylinders without plasmon properties (SiO

Thus, it can be stated that the strong inhomogeneity of the field in the case of a metal grating is a manifestation of the plasmon properties of its constituent elements.

For a more detailed analysis of the field structure, we considered the average electric field intensity in two sections: along the vector **E
** (blue dotted line in Fig. 15) and perpendicular to the vector

**E**0
(red dotted line in Fig. 15).

Figure 17 shows the dependence of the electric field intensity in a rectangular grating of nanocylinders in the cross section along the vector **E
**. It can be seen from the graphs that in the vacuum and in water the interaction of the electromagnetic wave with the grating occurs in different ways. This is evident from the structure of the field. It should be noted that the field is concentrated near the surface of the metal nanocylinder and is enhanced in comparison with the case when the grating consists of dielectric nanorods. In our case, the field intensity near the surface of the metal cylinder is approximately 5 times greater than the field intensity near the surface of the dielectric cylinder.

Figure 18 shows the dependence of the electric field intensity in a rectangular grating of nanocylinders in a section perpendicular to the vector **E
**

It should be noted that, in contrast to the quasi-static model, the wave decays. In this figure, as in the previous, you can observe field amplification near the surface of the conductive nanocylinders. In a grating of nanorods that do not possess plasmonic properties, the wave weakly decays and does not experience amplification.

It is seen from these dependences that the structure of a local field in a system of regular conducting nanocylinders calculated by different methods is qualitatively identical.

#### Acknowledgements

The authors acknowledge the support of The Ministry of Education and Science of Russia (№ 3.7758.2017/BP); the Russian Foundation for Basic Research (15-08-04132/17_а); Council President of the Russian Federation on grants for state support of young Russian scientists (SP-1340.2015.1). The authors are grateful to L.V. Terenina for assistance in the calculation of part 2.