Next Article in Journal
Yb to Tb Cooperative Upconversion in Supramolecularly Assembled Complexes in a Solution
Next Article in Special Issue
The B2 Structural Motif as a Tool for Modulating Ring Currents in Monocyclic Li Clusters
Previous Article in Journal / Special Issue
Spatial Contributions to 1H NMR Chemical Shifts of Free-Base Porphyrinoids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electronic Currents and Anapolar Response Induced in Molecules by Monochromatic Light

by
Francesco Ferdinando Summa
*,† and
Paolo Lazzeretti
Dipartimento di Chimica e Biologia “A. Zambelli”, Università Degli Studi di Salerno, Via Giovanni Paolo II 132, 84084 Fisciano, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 9 August 2021 / Revised: 30 August 2021 / Accepted: 1 September 2021 / Published: 5 September 2021

Abstract

:
It is shown that the electric dipole- and electric quadrupole–anapole polarizabilities, denoted respectively by f α β and g α , β γ , and the anapole magnetizability a α β , are intrinsic properties of the electron cloud of molecules responding to optical fields. f α β is a nonvanishing property of chiral and achiral compounds, whereas a α β is suitable for enantiomer discrimination of chiral species. They can conveniently be evaluated by numerical integration, employing a formulation complementary to that provided by perturbation theory and relying on the preliminary computation of electronic current density tensors all over the molecular domain. The origin dependence of the dynamic anapolar response is rationalized via related computational techniques employing numerical integration, as well as definitions of molecular property tensors, for example, electric dipole and electric quadrupole polarizabilties and magnetizability. A preliminary application of the theory is reported for the R a enantiomer of the hydrogen peroxide molecule, evaluating tensor components of electric dipole-anapole polarizability and anapole magnetizability as functions of the dihedral angle ϕ H-O-O-H in the range 0 ϕ 180 .

1. Introduction

A time-independent magnetic field induces toroidal electron flow in chiral molecules [1,2,3,4,5,6,7,8], giving rise to a static toroidal moment, which can be represented by a vector odd under charge conjugation C, parity P and time reversal T. The optical fields associated with a beam of monochromatic light, oscillating with frequency ω , also induce the anapolar response, which can be rationalized in terms of dynamic anapole polarizabilities and magnetizabilities [9].
The relevance of the anapolar response cannot be overemphasized after the advent of metamaterials, containing meta-atoms endowed with strong toroidal moments, used in a wide series of promising technological applications [10].
Within the electric quadrupole approximation, two second- and one third-rank tensors, depending on ω , have been introduced [9], namely anapole-electric dipole polarizability f α β , anapole magnetizability a α β , and anapole-electric quadrupole polarizability g α , β γ . Corresponding quantum-mechanical definitions have been reported within the framework of time-dependent perturbation theory [11], and their magnitude, as well as fundamental properties under C, P, and T symmetry operations [12], have been discussed. In particular, the origin dependence of these tensors has been investigated, determining equations that connect values corresponding to a passive translation of the coordinate system [9].
The present paper sets out to propose a derivation of f α β , a α β , and g α , β γ alternative to that arrived at via perturbation theory, allowing for optional definitions of quantum-mechanical anapole moments in terms of electronic current densities induced by monochromatic light shining on a molecule. Instead of going along an earlier trodden path [9], our investigation aims to reconceptualize a few aspects of the anapolar response. In fact, this novel approach is expected to throw light on fundamental physical mechanisms that come into play in determining the anapolar feedback, as it preludes to future visualization of current density vector fields, toroidal flow in particular.
The paper is organized as follows. An outline of notation is presented in Section 2 and a short description of the theoretical methods employed to compute anapolar response properties is provided in Section 3, introducing current density tensors which depend on the frequency ω of a beam of monochromatic light radiating on a probe molecule. An application of theoretical methods relying on induced current densities is discussed in Section 4 and concluding remarks are made in Section 5.

2. Outline of Notation and Theoretical Methods

Within the Born–Oppenheimer (BO) approximation [13], for a molecule with n electrons and N clamped nuclei, charge, mass, position and canonical momentum of the k-th electron are indicated, in the configuration space, by e , m e , r k , p ^ k , k = 1 , 2 , , n , using boldface letters for electronic vector operators. Analogous quantities for nucleus I are Z I e , M I , R I , and so forth, for I = 1 , 2 , , N .
The imaginary unit is represented by a Roman i . Throughout this paper, SI units are used and standard tensor formalism is employed; for example, the Einstein convention of implicit summation over two repeated Greek indices is applied. The third-rank Levi–Civita pseudotensor is indicated by ϵ α β γ . Capitals denote n-electron vector operators; for example, for position, canonical and angular momentum,
R ^ = k = 1 n r k , P ^ = k = 1 n p ^ k , L ^ = k = 1 n l ^ k ,
except for the electric and magnetic dipole,
(1) μ ^ = e R ^ , (2) m ^ = e 2 m e L ^ .
Electronic tensor operators are expressed, specifying the components; for example, for the second-rank electric quadrupole operator
μ ^ α β = e 2 k = 1 n ( r α r β ) k
for the diamagnetic contribution to magnetizability,
ξ ^ α β d = e 2 4 m e k = 1 n ( r 2 δ α β r α r β ) k
for the magnetic quadrupole operator in Hermitian form
m ^ α β = e 6 m e k = 1 n ( l ^ α r β + r β l ^ α ) k
and for the paramagnetic contribution to the anapole operator
a ^ γ = 1 2 ϵ α β γ m ^ α β = e 6 m e k = 1 n r 2 δ β γ r β r γ p ^ β + i r γ k = e 12 m e k = 1 n r 2 δ β γ r β r γ p ^ β + p ^ β r 2 δ β γ r β r γ k .
A “perturbed” operator, linearly dependent on an external magnetic field B , is defined via the relation,
a ^ α = a ^ α + a ^ α β d B β ,
where the diamagnetic contribution is given by
a ^ α β d = e 2 12 m e ϵ α β γ k = 1 n r 2 r γ k .
Expressions for the polarization charge density and current density induced in the electrons of a molecule by a monochromatic plane wave are obtained by time-dependent quantum mechanical perturbation theory [11], assuming that the eigenvalue problem for the time-independent BO electronic Hamiltonian H ^ ( 0 ) Ψ j ( 0 ) = E j ( 0 ) Ψ j ( 0 ) has been solved to determine a set of eigenfunctions Ψ j ( 0 ) and corresponding energy eigenvalues E j ( 0 ) . The reference (ground) state is indicated by Ψ a ( 0 ) and the natural transition frequencies are ω j a = E j ( 0 ) E a ( 0 ) / .
A general definition of n-electron density matrices is used throughout this work allowing for the McWeeny normalization [14],
γ x 1 ; x 1 = n Ψ x 1 , X 1 Ψ x 1 , X 1 d X 1 .
They are functions of electronic space-spin coordinates x k = r k s k , k = 1 , 2 , , n , where
X 1 { x 2 , , x n } , X = { x 1 , X 1 } , d X 1 { d x 2 , , d x n } ,
thus, integrating over d s 1 , one gets from Equation (9),
γ ( 0 ) ( r ) γ ( 0 ) ( r ; r ) = n Ψ a ( 0 ) r , X 1 Ψ a ( 0 ) r , X 1 d X 1
for the reference (ground) state Ψ a ( 0 ) of the molecule. The probability current density [14] is obtained from Equations (9)–(11) for the density matrix,
j ( r ) = 1 m e π ^ γ r ; r r = r ,
in which one puts r = r after operating with the electronic mechanical momentum
π ^ = p ^ + e A ,
adopting the Bloch gauge [15] for the vector potential A. The corresponding charge current density is obtained multiplying by e , that is, J = e j . The interaction Hamiltonian considered in the present work does not contain terms depending on electron spin, therefore the probability current density Equation (12) includes only orbital contributions.

3. Current Density Approach to Anapolar Response Properties and Their Origin Dependence

The time derivative of the electric field associated with a beam of monochromatic light impinging on a molecule induces an electronic current density J E ˙ ( r ) [16]. A corresponding current density tensor (CDT),
J α E ˙ β ( r , ω ) = n e m e j a ω j a 2 ω 2 1 × a μ ^ β j Ψ j ( 0 ) r , X 1 p ^ α Ψ a ( 0 ) r , X 1 d X 1 Ψ a ( 0 ) r , X 1 p ^ α Ψ j ( 0 ) r , X 1 d X 1 j μ ^ β a ,
to first order in the electric field, is obtained by differentiation,
J α E ˙ β ( r , ω ) = J α E ˙ ( r , ω ) E ˙ β .
Allowing for the long-wavelength approximation, only electric dipole terms are retained within the Goeppert–Mayer interaction Hamiltonian [17,18]. The electric field of the monochromatic wave is assumed spatially uniform, and Equation (14) defines a quantity independent of the origin of the coordinate system.
The mixed anapole-electric dipole polarizability induced by the time derivative of the electric field of monochromatic radiation can be cast in the form [9]
f α β = 1 j a 2 ω ω j a 2 ω 2 a a ^ α j j μ ^ β a = 1 6 ω r 2 δ α γ r α r γ J γ E ˙ β d 3 r .
The interesting feature of Equation (16), and of analogous relationships (34) and (41), is that they provide computational recipes alternative to the definitions determined by Rayleigh–Schrödinger perturbation theory [9] by expressing them as integrals of current density tensors. Thus anapole polarizabilities and magnetizability can be calculated by numerical integration using standard methods [19,20]. In this work, spatial integration of the density functions has been performed using the Becke scheme, Ref. [21] adopting 131 angular points for the Lebedev’s quadrature of 59th order of accuracy [22] and 131 radial points for the Gauss–Chebyshev radial quadrature of the second-kind [23].
According to Equation (16), the behavior of the f α β polarizability under the fundamental symmetries is different from that of the anapole magnetizability a α β [5]. Whereas the former is even under C, P and T, the latter is C- and T-even and P-odd, therefore only a α β can be used for chiral discrimination [9].
Since the CDT in Equation (14) is invariant under a translation of coordinates, the origin dependence of f α β is only due to the prefactor within brackets in Equation (16), which will change in a passive transformation of reference system,
r r = r + d ,
corresponding to an arbitrary shift d of the origin. In fact, a simple calculation allowing for the second identity of (16) shows that two different observers, at the origin of reference systems K and K , would calculate, or measure, different anapole-electric dipole polarizabilities related by
f α β ( r ) = f α β ( r ) + ω 3 d δ ( r δ r δ ) J α E ˙ β d 3 r ω 6 d α ( r γ r γ ) J γ E ˙ β d 3 r ω 6 ( d 2 δ α γ d α d γ ) J γ E ˙ β d 3 r ω 6 d γ ( r α r α ) J γ E ˙ β d 3 r .
The fourth term on the r.h.s. of Equation (18) is soon evaluated via Equation (90) of Ref. [24], as the argument (14) of the integral
J α E ˙ β ( r , ω ) d 3 r = α β α ( R , P ) = e 2 m e j a 2 ω j a 2 ω 2 a R ^ β j j P ^ α a
defines an origin-independent density of dynamic electric dipole polarizability in the ( R , P ) gauge, identical to that in the ( R , R ) and ( P , P ) gauges
α α β ( R , R ) = 1 j a 2 ω j a ω j a 2 ω 2 a μ ^ α j j μ ^ β a ,
α α β ( P , P ) = e 2 m e 2 j a 2 ω j a ( ω j a 2 ω 2 ) a P ^ α j j P ^ β a ,
if the hypervirial theorem [25,26,27]
a P ^ α j = i m e ω j a a R ^ α j , j P ^ α a = i m e ω j a j R ^ α a ,
is satisfied [28].
Calculation of the other integrals appearing in Equation (18) is greatly facilitated by employing a general and quite effective relationship for the moments of total electronic current density [24]. According to Equation (82) of that paper, to first-order in the perturbing fields of the monochromatic wave,
r α J β ( 1 ) d 3 r = d d t Δ μ ^ α β + ϵ α β γ Δ m ^ γ .
Within the algebraic approximation [29], right- and left-hand sides of Equation (23) are expected to yield different numerical results, unless the basis set is close to completeness, or flexible enough to provide accurate representation of R ^ , P ^ and L ^ at the same time, and to satisfy hypervirial theorems [25,26,27] to a good extent. In fact, the current density tensors in Equations (14), (26) and (27) are expressed in terms of p ^ , whereas Δ μ ^ α β is usually evaluated within the dipole length gauge in Equation (28). Thus the degree to which the equality (23) is fulfilled provides a measure of basis set quality.
This relationship is applied in Equation (18) by taking into account terms that depend on the same perturbing fields at both sides. Therefore, within the electric quadrupole approximation [16], it is properly used relaxing the long wave-length assumption, that is, supposing that J β E ˙ ( r , ω ) is not origin independent, since it varies according to the change of electric field at different points within the molecular region,
E β r = E β r + d γ E γ β .
At any rate, in Equations (23) and (24), both the electric field gradient γ E β E γ β and the magnetic field of the impinging radiation are assumed to be spatially uniform all over the molecular domain. They are coupled with the electric quadrupole and magnetic dipole operators, respectively, within the interaction Hamiltonian. This amounts to considering also the current density vector fields induced by the gradient E ˙ and the magnetic field B , that is,
J α ( 1 ) = J α E ˙ + J α E ˙ + J α B ,
in agreement with the electric quadrupole approximation.
The current density terms appearing in Equation (25) are taken into account introducing the corresponding CDTs [16],
J α E ˙ β γ ( r , ω ) = n e m e j a ω j a 2 ω 2 1 a μ ^ β γ j × Ψ j ( 0 ) r , X 1 p ^ α Ψ a ( 0 ) r , X 1 d X 1 Ψ a ( 0 ) r , X 1 p ^ α Ψ j ( 0 ) r , X 1 d X 1 j μ ^ β γ a ,
and
J α B β ( r , ω ) = n e m e j a ω j a ω j a 2 ω 2 × a m ^ β j Ψ j ( 0 ) r , X 1 p ^ α Ψ a ( 0 ) r , X 1 d X 1 + Ψ a ( 0 ) r , X 1 p ^ α Ψ j ( 0 ) r , X 1 d X 1 j m ^ β a e 2 2 m e ϵ α β γ r γ γ ( 0 ) ( r ) .
On the r.h.s. of Equation (23), the electric quadrupole induced in the electronic cloud is expressed, within the electric quadrupole approximation, in terms of intrinsic tensor properties; see Appendix A,
Δ μ ^ α β = α γ , α β E γ + α α β , γ δ E δ γ D γ , α β B ˙ γ ω 1 .
An analogous expression is used for the induced magnetic dipole:
Δ m ^ α = κ β α E ˙ β ω 1 + D α , β γ E ˙ γ β ω 1 + ξ α β B β .
From Equations (18) and (23), one finds
(30) ( r δ r δ ) J α E ˙ β d 3 r = α β , δ α ( r ) ϵ δ α γ κ β γ ( r ) , (31) ( r γ r γ ) J γ E ˙ β d 3 r = α β , γ γ ( r ) , (32) ( r α r α ) J γ E ˙ β d 3 r = α β , γ α ( r ) ϵ α γ λ κ β λ ( r ) ,
hence the relationship connecting the anapole polarizability f α β for the different origins becomes
f α β r = f α β r ω 6 α β γ d 2 δ α γ d α d γ + 1 2 κ β δ r ϵ α γ δ d γ + ω 6 α β , α γ r d γ α β , γ γ r d α ,
the same as that arrived at via a more complicated procedure [9], allowing for the translation of operators (6) and (7). Thus the origin dependence of f α β is fully determined by the prefactor in round brackets in Equation (16).
The mixed anapole-electric quadrupole polarizability induced by the time derivative of the gradient E γ β associated with the monochromatic beam [9] can be expressed in the form:
g α , β γ = 1 j a 2 ω ω j a 2 ω 2 a a ^ α j j μ ^ β γ a = 1 6 ω r 2 δ α δ r α r δ J δ E ˙ β γ d 3 r ,
even under C and T and odd under P, hence suitable for chiral recognition. A procedure allowing for the second identity in Equation (34), analogous to that leading to (18), gives
g α , β γ ( r ) = g α , β γ ( r ) + ω 3 d ϵ ( r ϵ r ϵ ) J α E ˙ β γ d 3 r ω 6 d α ( r δ r δ ) J δ E ˙ β γ d 3 r ω 6 ( d 2 δ α δ d α d δ ) J δ E ˙ β γ d 3 r ω 6 d δ ( r α r α ) J δ E ˙ β γ d 3 r .
The integral in the fourth term on the r.h.s. of Equation (35) corresponds to the electric dipole-electric quadrupole polarizability in the dipole velocity gauge [24]:
J α E ˙ γ β ( r , ω ) d 3 r = α α , β γ ( P ) = e m e j a 2 ω j a 2 ω 2 a μ ^ β γ j j P ^ α a ,
identical to that in the dipole length gauge, reported in the Appendix A, Equation (A1), if the hypervirial theorem (22) is fulfilled.
From Equations (23), (28) and (29) one finds
(37) ( r ϵ r ϵ ) J α E ˙ β γ d 3 r = α β γ , ϵ α ( r ) ϵ ϵ δ α D δ , β γ ( r ) , (38) ( r δ r δ ) J δ E ˙ β γ d 3 r = α β γ , δ δ ( r ) , (39) ( r α r α ) J δ E ˙ β γ d 3 r = α β γ , α δ ( r ) ϵ α ϵ δ D δ , β γ ( r ) ;
then the relationship connecting the anapole-electric quadrupole polarizability g α , β γ for the different origins becomes:
g α , β γ r = g α , β γ r ω 6 α δ , β γ r d 2 δ α δ d α d δ + ω 6 α β γ , α δ r d δ α β γ , δ δ r d α + 1 2 ϵ α δ ϵ D δ , β γ r d ϵ .
Thus the prefactor within round brackets in Equation (34) yields only a few terms linear and quadratic in the d shift, as shown by comparison with Equation (96) of Ref. [9]. Missing contributions would be recovered by coupling the prefactor to terms arising from a change of origin of the J α E ˙ β γ CDT. In fact, only the total first order current density (25) is origin independent, whereas the three contributions on the r.h.s. interchange in a passive transformation of the coordinate system [16].
The anapole magnetizability is evaluated via a magnetic-field induced current density tensor (27),
(41) a γ δ = 1 6 r 2 δ α γ r α r γ J α B δ d 3 r = a γ δ p + a γ δ d , (42) a γ δ p ( ω ) = 1 j a 2 ω j a ω j a 2 ω 2 a a ^ γ j j m ^ δ a , (43) a γ δ d = a | a ^ γ δ d | a ,
and a straightforward calculation using the first identity in Equation (41) for different origins yields
a γ δ ( r ) = a γ δ ( r ) + 1 3 d β ( r β r β ) J γ B δ d 3 r 1 6 d α ( r γ r γ ) J α B δ d 3 r 1 6 d γ ( r α r α ) J α B δ d 3 r 1 6 ( d 2 δ α γ d α d γ ) J α B δ d 3 r .
The fifth term on the r.h.s. is expressed via the magnetic CDT (27), which is related to a density of electric dipole-magnetic dipole polarizability [24],
J α B β ( r , ω ) d 3 r = ω κ α β ,
and the other integrals are calculated allowing for Equations (23), (28) and (29). Thus one finds
(46) ( r β r β ) J γ B δ d 3 r = ω D δ , β γ r + ϵ β γ λ ξ λ δ (47) ( r γ r γ ) J α B δ d 3 r = ω D δ , α γ r + ϵ β γ λ ξ λ δ (48) ( r α r α ) J α B δ d 3 r = ω D δ , α α r ,
then the relationship connecting the anapole magnetizability a γ δ for the different origins becomes
a γ δ r = a γ δ r + ω 6 d α D δ , α γ r d γ D δ , α α r + 1 2 d β ϵ β γ λ ξ λ δ ω 6 d 2 δ α γ d α d γ κ α δ r .
Therefore, for the anapole magnetizability (41), only a few terms of those appearing in Equation (93) of Ref. [9] are recovered, for the reasons explained in the discussion following Equation (35).
For static anapole magnetizabilties, Equation (44) becomes [5]
a γ δ r = a γ δ r + 1 2 ϵ α β γ ξ δ α d β ,
thus the trace a γ γ is invariant of the origin for ω = 0 . For a freely tumbling chiral molecule in disordered media, the isotropic anapole moment induced in the direction of the static magnetic field is Δ a ^ α = ( 1 / 3 ) a γ γ B α , whereas the isotropic magnetic dipole moment induced by the curl of a static magnetic field B , that is, C = × B , is Δ m ^ α = ( 1 / 3 ) a γ γ C α .

4. Results and Discussion

A preliminary test of the theoretical scheme outlined in Section 2 and Section 3 has been conducted on the hydrogen peroxide molecule, by computing the anapole-electric dipole polarizability f α β , Equation (16), and anapole magnetizability a α β , Equation (41). The origin dependence of the former has been investigated via Equations (18) and (33), which account for all the terms reported in a previous study [9].
Calculations have been carried out for the R a enantiomer of H 2 O 2 shown in Figure 1, evaluating both f α β and a α β tensor components as functions of the dihedral angle ϕ = H-O-O-H, at the time dependent Hartree–Fock (TD-HF) level of theory, equivalent to random-phase approximation (RPA), adopting a fairly large basis set, that is, the uncontracted d-aug-cc-pVQZ [30,31,32] on hydrogen atoms and d-aug-cc-pVTZ [30,31,32] on oxygen atoms. Both basis sets have been downloaded from BSE [33,34]. In our implementation, outlined in detail in a previous reference [35], transition amplitudes S j and T j and corresponding transition energies have been obtained by means of TD-HF calculations. The complete procedure for computing frequency dependent anapole polarizabilities f α β and anapole magnetizabilities a α β has been coded within the freely available SYSMOIC program package [36].
The results obtained via numerical integration [21] of Equations (16) and (41) for two wavelengths, λ = 355 nm and λ = 589.3 nm, are displayed in the interval 0 ϕ 180 in Figure 2 and Figure 3 for the anapole polarizability and in Figure 4 and Figure 5 for the anapole magnetizability.
As can be observed in Figure 2 and Figure 3, the f y y component and the isotropic value f α α / 3 of the anapole-electric dipole polarizability are nearly independent of the dihedral angle ϕ H-O-O-H for both λ values, that is, insensitive to rotation about the O-O bond, which nearly coincides with the y axis. Thus, the anapole moment induced in the y direction by the time derivative of an electric field E y does not appreciably change on varying ϕ . The average anapole moment Δ a ^ = ( 1 / 3 ) f α α E ˙ is similarly unaffected by rotation.
Quite remarkably, the curves representing isotropic anapole magnetizability in Figure 4 and Figure 5, for the frequencies λ = 355 nm and λ = 589.3 nm, vanish in the proximity of a dihedral angle ϕ 100 . The presence of an accidental null value of a α α for a geometrically chiral conformation of hydrogen peroxide, also observed in previous studies of its static anapole magnetizability [37], is due to a different sign of diagonal components: a x x is positive and quite large for ϕ 100 , offsetting the algebraic sum of a y y and a z z .
A test of accuracy of computations is provided by results obtained via numerical integration, reported in Table 1. The quality of the basis sets adopted, and of its ability to give quite reliable representations of operators R ^ , P ^ and L ^ , is evident by comparison of electric dipole polarizabilities α α β in different gauges ( R , R ) , ( P , R ) and ( P , P ) , as well as electric dipole-magnetic dipole polarizabilities κ α β in ( R , L ) and ( P , L ) formalisms: computed tensor components are nearly coincident to three or more significant figures.
A similar conclusion is arrived at by inspection of Table 2, reporting electric dipole-electric quadrupole polarizability within dipole length ( R ) and dipole velocity ( P ) formalisms, Equations (36) and (A1).
A further proof of the near Hartree–Fock accuracy of the present calculations is obtained by comparing results from either side of Equation (23), fulfilling the off-diagonal hypervirial relations (22) quite nicely. The electric dipole electric anapole f α β tensor obtained by the integration of Equation (18) in the origin r = ( 1 , 2 , 3 ) is
f α β ( r ) = 2.3952 0.0103 0.5099 0.0559 2.5953 0.7955 0.7178 1.6330 1.1715 ,
virtually identical to that obtained via Equation (33),
f α β ( r ) = 2.3952 0.0104 0.5099 0.0560 2.5955 0.7955 0.7177 1.6330 1.1714 .
A broad discussion about the off-diagonal hypervirial relationship, concerning their fulfillment within incomplete basis sets, can be found in Ref. [28]. Quite small differences among ( R , R ) , ( R , P ) and ( P , P ) gauges should only be imputed to basis set quality and not to the numerical integration method.

5. Concluding Remarks and Outlook

The numerical results, arrived at in the present study within the theoretical scheme outlined in Section 3 and the discussion in Section 4, show that frequency dependent anapolar properties of chiral and achiral molecules responding to monochromatic radiation can be easily computed by simple relationships expressed in terms of induced electronic current densities.
Analogous procedures were found to be applicable to study the origin dependence of anapole polarizabilities and magnetizability.
The method proposed and implemented in the present work is an alternative to the conventional formulation based on perturbation theory. It appears to be practical to apply and is computationally more advantageous, since it avails itself of numerical integration, and is quite easy to perform whenever three-dimensional grids of values are available for the current density vector fields induced by the magnetic field, the time derivative of the electric field and of the electric field gradient associated with the impinging radiation. Visualizations of these current densities, which will be available in future investigations, are expected to yield very useful tools for the rationalization of mechanisms coming into play and giving rise to the anapolar response of electrons in a molecule. Effects due to electron correlation will also be investigated in future work.

Author Contributions

Conceptualization, P.L.; methodology, F.F.S. and P.L.; software, F.F.S.; validation, F.F.S. and P.L.; formal analysis, F.F.S. and P.L.; investigation, F.F.S. and P.L.; resources, F.F.S. and P.L.; data curation, F.F.S. and P.L.; writing—original draft preparation, P.L.; writing—review and editing, F.F.S. and P.L.; visualization, F.F.S.; supervision, P.L. All authors have read and agreed to the published version of the manuscript.

Funding

Financial support from MIUR is gratefully acknowledged.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank Professor G. Monaco for helpful discussions. This paper is dedicated to Professor Riccardo Zanasi on the occasion of his 70th birthday.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Electric dipole-electric quadrupole polarizability
α α , β γ = 1 j a 2 ω j a ω j a 2 ω 2 a μ ^ α j j μ ^ β γ a .
Electric quadrupole-electric quadrupole polarizability
α α β , γ δ = 1 j a 2 ω j a ω j a 2 ω 2 a μ ^ α β j j μ ^ γ δ a .
Magnetic dipole-electric quadrupole polarizability
D α , β γ = 1 j a 2 ω ω j a 2 ω 2 a m ^ α j j μ ^ β γ a .
Electric dipole-magnetic dipole polarizability
κ α β = 1 j a 2 ω ω j a 2 ω 2 a μ ^ α j j m ^ β a .
Magnetizability, diamagnetic and paramagnetic contributions,
(A5) ξ α β = ξ α β d + ξ α β p , ξ α β d = a ξ ^ α β d a , (A6) ξ α β p = 1 j a 2 ω j a ω j a 2 ω 2 a m ^ α j j m ^ β a .
Anapole-electric dipole polarizability
f α β ( ω ) = 1 j a 2 ω ω j a 2 ω 2 a a ^ α j j μ ^ β a .
Anapole-electric quadrupole polarizability
g α , β γ ( ω ) = 1 j a 2 ω ω j a 2 ω 2 a a ^ α j j μ ^ β γ a .
Anapole magnetizability, paramagnetic and diamagnetic contributions:
(A9) a α β p = 1 j a 2 ω j a ω j a 2 ω 2 a a ^ α j j m ^ β a , (A10) a α β d = a a ^ α β d a = e 2 12 m e ϵ α β γ a k = 1 n r 2 r γ k a , (A11) a α β = a α β p + a α β d .

References

  1. Pelloni, S.; Lazzeretti, P.; Monaco, G.; Zanasi, R. Magnetic-field induced electronic anapoles in small molecules. Rend. Lincei 2011, 22, 105–112. [Google Scholar] [CrossRef]
  2. Berger, R. Prediction of a Cyclic Helical Oligoacetylene Showing Anapolar Ring Currents in the Magnetic Field. Z. Naturforsch. 2012, 67, 1127–1131. [Google Scholar] [CrossRef]
  3. Tellgren, E.I.; Fliegl, H. Non-perturbative treatment of molecules in linear magnetic fields: Calculation of anapole susceptibilities. J. Chem. Phys. 2013, 139, 1–15. [Google Scholar] [CrossRef] [PubMed]
  4. Pagola, G.I.; Ferraro, M.B.; Provasi, P.F.; Pelloni, S.; Lazzeretti, P. Theoretical estimates of the anapole magnetizabilities of C4H4X2 cyclic molecules for X=O, S, Se, and Te. J. Chem. Phys. 2014, 141, 094305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zarycz, N.; Provasi, P.F.; Pagola, G.I.; Ferraro, M.B.; Pelloni, S.; Lazzeretti, P. Computational study of basis set and electron correlation effects on anapole magnetizabilities of chiral molecules. J. Comp. Chem. 2016, 37, 1552–1558. [Google Scholar] [CrossRef] [PubMed]
  6. Reiter, K.; Weigend, F.; Wirz, L.N.; Dimitrova, M.; Sundholm, D. Magnetically Induced Current Densities in Toroidal Carbon Nanotubes. J. Phys. Chem. C 2019, 123, 15354–15365. [Google Scholar] [CrossRef] [Green Version]
  7. Orozco-Ic, M.; Dimitrova, M.; Barroso, J.; Sundholm, D.; Merino, G. Magnetically Induced Ring-Current Strengths of Planar and Nonplanar Molecules: New Insights from the Pseudo-π Model. J. Phys. Chem. A 2021, 125, 5753–5764. [Google Scholar] [CrossRef]
  8. Pagola, G.I.; Ferraro, M.B.; Provasi, P.F.; Pelloni, S.; Lazzeretti, P. Physical achirality in geometrically chiral rotamers of hydrazine and boranylborane molecules. J. Comp. Chem. 2021, 42, 1772. [Google Scholar] [CrossRef] [PubMed]
  9. Lazzeretti, P. Static and optical anapole magnetizabilities and polarizabilities. J. Chem. Phys. 2020, 153, 074102. [Google Scholar] [CrossRef]
  10. Nanz, S. Toroidal Multipole Moments in Classical Electrodynamics; Springer: Wiesbaden, Germany, 2016. [Google Scholar] [CrossRef]
  11. Langhoff, P.W.; Epstein, S.T.; Karplus, M. Aspects of Time-Dependent Perturbation Theory. Rev. Mod. Phys. 1972, 44, 602–644. [Google Scholar] [CrossRef]
  12. Lazzeretti, P. The abstract GPT and GCPT groups of discrete C, P and T symmetries. J. Mol. Spectrosc. 2017, 337, 178–184. [Google Scholar] [CrossRef]
  13. Born, M.; Oppenheimer, R. Zur Quantentheorie der Molekeln. Ann. Phys. 1927, 389, 457–484. [Google Scholar] [CrossRef]
  14. McWeeny, R. Methods of Molecular Quantum Mechanics, 2nd ed.; Theoretical Chemistry, OCLC: 247760584; Academic Press: London, UK, 1992. [Google Scholar]
  15. Bloch, F. Zur Wirkung äußerer elektromagnetischer Felder auf kleine Systeme. In W. Heisenberg und die Physik Unserer Zeit; Bopp, F., Ed.; Friedr. Wieveg & Son: Braunschweig, Germany, 1961; pp. 93–102. [Google Scholar]
  16. Lazzeretti, P. Gauge invariance and origin independence of electronic charge density and current density induced by optical fields. J. Chem. Phys. 2018, 149, 154106. [Google Scholar] [CrossRef]
  17. Göppert-Mayer, M. Über Elementarakte mit zwei Quantensprüngen. Ann. Phys. 1931, 401, 273–294. [Google Scholar] [CrossRef]
  18. Cohen-Tannoudji, C.; Dupont-Roc, J.; Grynberg, G. Photons and Atoms: Introduction to Quantum Electrodynamics; Wiley: New York, NY, USA, 1989. [Google Scholar]
  19. Jinger, R.K.; Fliegl, H.; Bast, R.; Dimitrova, M.; Lehtola, S.; Sundholm, D. Spatial Contributions to Nuclear Magnetic Shieldings. J. Phys. Chem. A 2021, 125, 1778–1786. [Google Scholar] [CrossRef] [PubMed]
  20. Lehtola, S.; Dimitrova, M.; Fliegl, H.; Sundholm, D. Benchmarking Magnetizabilities with Recent Density Functionals. J. Chem. Theory Comput. 2021, 17, 1457–1468. [Google Scholar] [CrossRef] [PubMed]
  21. Becke, A.D. A multicenter numerical integration scheme for polyatomic molecules. J. Chem. Phys. 1988, 88, 2547–2553. [Google Scholar] [CrossRef]
  22. Lebedev, V. A quadrature formula for the sphere of 59th algebraic order of accuracy. Dokl. Akad. Nauk 1994, 338, 454–456. [Google Scholar]
  23. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables, 9th ed.; Dover Books on Mathematics; Dover Publ: New York, NY, USA, 2013. [Google Scholar]
  24. Lazzeretti, P. Frequency-dependent current density tensors as density functions of dynamic polarizabilities. J. Chem. Phys. 2019, 150, 184117. [Google Scholar] [CrossRef] [PubMed]
  25. Chen, J.C.Y. Off–Diagonal Hypervirial Theorem and Its Applications. J. Chem. Phys. 1964, 40, 615–621. [Google Scholar] [CrossRef]
  26. Chong, D.P.; Benston, M.L. Off-Diagonal Hypervirial Theorems as Constraints. J. Chem. Phys. 1968, 49, 1302–1306. [Google Scholar] [CrossRef]
  27. Epstein, S.T. The Variation Method in Quantum Chemistry; Number v. 33 in Physical Chemistry, a Series of Monographs; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  28. Summa, F.F.; Monaco, G.; Lazzeretti, P.; Zanasi, R. Assessment of the Performance of DFT Functionals in the Fulfillment of Off-Diagonal Hypervirial Relationships. Phys. Chem. Chem. Phys. 2021, 23, 15268–15274. [Google Scholar] [CrossRef]
  29. Moncrieff, D.; Wilson, S. On the accuracy of the algebraic approximation in molecular electronic structure calculations. III. Comparison of matrix Hartree-Fock and numerical Hartree-Fock calculations for the ground state of the nitrogen molecule. J. Phys. B At. Mol. Opt. Phys. 1993, 26, 1605. [Google Scholar] [CrossRef]
  30. Dunning, T.H. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  31. Woon, D.E.; Dunning, T.H. Gaussian basis sets for use in correlated molecular calculations. IV. Calculation of static electrical response properties. J. Chem. Phys. 1994, 100, 2975–2988. [Google Scholar] [CrossRef] [Green Version]
  32. Kendall, R.A.; Dunning, T.H.; Harrison, R.J. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 1992, 96, 6796–6806. [Google Scholar] [CrossRef] [Green Version]
  33. Pritchard, B.P.; Altarawy, D.; Didier, B.; Gibson, T.D.; Windus, T.L. New Basis Set Exchange: An Open, Up-to-Date Resource for the Molecular Sciences Community. J. Chem. Inf. Model. 2019, 59, 4814–4820. [Google Scholar] [CrossRef]
  34. Schuchardt, K.L.; Didier, B.T.; Elsethagen, T.; Sun, L.; Gurumoorthi, V.; Chase, J.; Li, J.; Windus, T.L. Basis Set Exchange: A Community Database for Computational Sciences. J. Chem. Inf. Model. 2007, 47, 1045–1052. [Google Scholar] [CrossRef] [Green Version]
  35. Summa, F.F.; Monaco, G.; Zanasi, R.; Pelloni, S.; Lazzeretti, P. Electronic Currents Induced by Optical Fields and Rotatory Power Density in Chiral Molecules. Molecules 2021, 26, 4195. [Google Scholar] [CrossRef]
  36. Monaco, G.; Summa, F.F.; Zanasi, R. Program Package for the Calculation of Origin-Independent Electron Current Density and Derived Magnetic Properties in Molecular Systems. J. Chem. Inf. Model. 2021, 61, 270–283. [Google Scholar] [CrossRef]
  37. Pelloni, S.; Provasi, P.F.; Pagola, G.I.; Ferraro, M.B.; Lazzeretti, P. Electric Dipole–Magnetic Dipole Polarizability and Anapole Magnetizability of Hydrogen Peroxide as Functions of the HOOH Dihedral Angle. J. Phys. Chem. A 2017, 121, 9369–9376. [Google Scholar] [CrossRef] [PubMed]
Figure 1. R α enantiomer of the hydrogen peroxide molecule.
Figure 1. R α enantiomer of the hydrogen peroxide molecule.
Chemistry 03 00073 g001
Figure 2. Tensor components of anapole polarizability f α β for the R α enantiomer of the hydrogen peroxide molecule computed at the TD-HF level by Equation (16) for λ = 355 nm, as a function of the dihedral angle.
Figure 2. Tensor components of anapole polarizability f α β for the R α enantiomer of the hydrogen peroxide molecule computed at the TD-HF level by Equation (16) for λ = 355 nm, as a function of the dihedral angle.
Chemistry 03 00073 g002
Figure 3. Tensor components of anapole polarizability f α β for the R α enantiomer of the hydrogen peroxide molecule at λ = 589.3 nm, as a function of the dihedral angle.
Figure 3. Tensor components of anapole polarizability f α β for the R α enantiomer of the hydrogen peroxide molecule at λ = 589.3 nm, as a function of the dihedral angle.
Chemistry 03 00073 g003
Figure 4. Tensor components of anapole magnetizability a α β for the R α enantiomer of the hydrogen peroxide molecule computed at the TD-HF level by Equation (41) for λ = 355 nm, as a function of the dihedral angle.
Figure 4. Tensor components of anapole magnetizability a α β for the R α enantiomer of the hydrogen peroxide molecule computed at the TD-HF level by Equation (41) for λ = 355 nm, as a function of the dihedral angle.
Chemistry 03 00073 g004
Figure 5. Tensor components of the anapole magnetizability a α β for the R α enantiomer of the hydrogen peroxide molecule at λ = 589.3 nm, as a function of the dihedral angle.
Figure 5. Tensor components of the anapole magnetizability a α β for the R α enantiomer of the hydrogen peroxide molecule at λ = 589.3 nm, as a function of the dihedral angle.
Chemistry 03 00073 g005
Table 1. Molecular tensors evaluated at λ = 589.3 nm, in atomic units.
Table 1. Molecular tensors evaluated at λ = 589.3 nm, in atomic units.
xx xy yx yy zz Av
f α β ( r ) −0.4788−0.19860.0441−0.1374−0.4130−0.3431
a α β ( r ) 0.09880.09920.0207−0.0136−0.1172−0.0107
α α β ( R , R ) ( r ) 11.99981.60391.603917.799511.056013.6184
α α β ( R , P ) ( r ) 11.98891.60141.602217.777911.047113.6046
α α β ( P , R ) ( r ) 11.98891.60221.601417.777911.047113.6046
α α β ( P , P ) ( r ) 11.97831.59971.599717.756611.038513.5911
κ α β ( R , L ) ( r ) −0.0873−0.0006−0.16010.01570.0668−0.0016
κ α β ( P , L ) ( r ) −0.0874−0.0006−0.15980.01570.0670−0.0016
Conversion factor for f α β : e 2 a 0 4 / = 1.908750473 × 10 45 F m 4 s 1 ; Conversion factor for a α β : e 2 a 0 3 / m e = 4.175756730 × 10 39 J T 2 m ; Conversion factor for α α β : e 2 a 0 2 / E h = 1.648777273 × 10 41 F m 2 ; Conversion factor for κ α β : e 2 a 0 3 / = 3.607015578 × 10 35 F m 3 s 1 .
Table 2. Mixed electric dipole-electric quadrupole polarizability within ( R ) and ( P ) formalisms, in atomic units.
Table 2. Mixed electric dipole-electric quadrupole polarizability within ( R ) and ( P ) formalisms, in atomic units.
RP
α z , x x 0.97710.9768
α z , y y 0.22700.2283
α z , z z 0.49200.4917
α z , x y 0.98260.9806
α x , x z 1.28971.2893
α x , y z 1.04181.0395
α y , x z 0.57110.5702
α y , y z 0.81970.8202
The conversion factor is e 2 a 0 3 / E h = 8.724953446 × 10 52 F m 3 .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Summa, F.F.; Lazzeretti, P. Electronic Currents and Anapolar Response Induced in Molecules by Monochromatic Light. Chemistry 2021, 3, 1022-1036. https://0-doi-org.brum.beds.ac.uk/10.3390/chemistry3030073

AMA Style

Summa FF, Lazzeretti P. Electronic Currents and Anapolar Response Induced in Molecules by Monochromatic Light. Chemistry. 2021; 3(3):1022-1036. https://0-doi-org.brum.beds.ac.uk/10.3390/chemistry3030073

Chicago/Turabian Style

Summa, Francesco Ferdinando, and Paolo Lazzeretti. 2021. "Electronic Currents and Anapolar Response Induced in Molecules by Monochromatic Light" Chemistry 3, no. 3: 1022-1036. https://0-doi-org.brum.beds.ac.uk/10.3390/chemistry3030073

Article Metrics

Back to TopTop