Surface-Enhanced Circular Dichroism of Oriented Chiral Molecules by

Dec 15, 2016 - We present a rigorous finite element method to calculate circular dichroism (CD) in various systems consisting of nanostructures and or...
3 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF REGINA

Article

Surface Enhanced Circular Dichroism of Oriented Chiral Molecules by Plasmonic Nanostructures Weixuan Zhang, Tong Wu, Rong-Yao Wang, and Xiangdong Zhang J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b09435 • Publication Date (Web): 15 Dec 2016 Downloaded from http://pubs.acs.org on December 21, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Surface Enhanced Circular Dichroism of Oriented Chiral Molecules by Plasmonic Nanostructures Weixuan Zhang, Tong Wu, Rongyao Wang and Xiangdong Zhang* Beijing Key Laboratory of Nanophotonics & Ultrafine Optoelectronic Systems, School of Physics, Beijing Institute of Technology, Beijing, 100081, China ABSTRACT: We present a rigorous finite element method to calculate circular dichroism (CD) in various systems consisting of nanostructures and oriented chiral molecules with electric quadrupole transitions. The interaction between oriented molecule materials, which are regarded as anisotropic chiral media, and metallic nanostructures has been investigated. Our results show that the plasmon-induced CD is sensitive to the orientations of the molecules. In many cases, the contribution of molecular electric quadrupole transitions to the total CD signal can play a key role. More interesting, we have demonstrated that both the quadrupole- and dipole-based CD signals can be improved greatly by matching the phases for the electromagnetic fields and their gradients at different regions around the nanostructures, which are occupied by the oriented chiral molecules. Different regions might produce CD of opposite sign. When integrating over regions with only one side of the proposed nanostructure, we find that the CD-peak may be nearly hundreds-fold over the case of integrating both sides. We believe that these findings would be helpful for realizing ultrasensitive probing of chiral information for oriented molecules by plasmon-based nanotechnology.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 32

I.INTRODUCTION Chirality, which refers to structures lacking any mirror symmetry planes, is a very intriguing property of molecules. Many biologically active molecules are chiral, which plays a pivotal role in biochemistry and the evolution of life itself.1-2 Detecting and characterizing chiral enantiomers of these biomolecules are of considerable importance for biomedical diagnostics and pathogen analyses.3 A common technique for chirality discrimination is CD spectroscopy describing the difference in molecular absorption of left- and right- handed circularly polarized photons.4-5 In general, the molecular CD signal is typically weak, thus, chiral analyses by such a spectroscopic technique have usually been restricted to analysis at a relatively high concentration.1-5 Recent investigations have shown that superchiral fields allow measuring the chiroptical properties of small amounts of molecules with high sensitivity.6-7 A weak molecular CD signal in the ultraviolet spectral region can be both enhanced and transferred to the visible-near-infrared region when chiral molecules are adsorbed at the surfaces of metallic nanoparticles or in the nanogaps (i.e., hot spots) of particle clusters.8-29 However, many related discussions assume that the molecules are randomly oriented, 6, 21-24 and the results have been obtained by averaging over the solid angles of the molecular directions. In fact, this rotational degree of freedom for molecules is absent in many cases, i.e. chiral molecules adopt geometries in which they have an axis with a well-defined orientation with respect to the surface of the nanostructures.9, 26-29 This leads to several orientation-selective signals for appropriately sculpted fields. Although there are also some theoretical researches concentrating on the plasmon-induced CD effect of oriented chiral molecules, they only focus on a fixed number of chiral molecules modeled as point dipoles and the molecular electric quadrupole transitions have also been neglected.17-20 As for the chiral molecular material with a large collections of molecules, it can be generally regarded as a

ACS Paragon Plus Environment

2

Page 3 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

homogeneous chiral medium being treated with full-field simulations.24, 30-31 In addition, up to now, all the studies have neglected the molecule electric quadrupolar contribution to the CD.17, 19, 20-22, 24

When the molecules are randomly oriented, the neglect is reasonable. Nevertheless, if the

molecules have defined orientations, the contribution of molecular electric quadrupole transitions to the total CD signal can play an important role.26, 27, 32-34 For example, the large dissymmetries observed for the experiment proposed by E. Hendry et al cannot be explained by dipolar chiral excitations.26 The large dissymmetry enhancements may result from the quadrupolar contribution to optical activity. Hence, it is not correct to ignore the molecular electric quadrupolar contribution in the oriented molecular systems. In this work, we present a rigorous finite element method to study electromagnetic interactions between nanostructures and anisotropic chiral molecular materials with the electric quadrupolar contribution. On the basis of such a method, we calculate the CD in the systems consisting of nanostructures and oriented molecules. We find that the molecular electric dipoleelectric quadrupole interaction can play primary roles in the plasmon-induced CD spectroscopy for the system with specific molecule orientations. Furthermore, we have also found that both the quadrupole- and dipole-based CD signals can be amplified tremendously by matching the phases for the electromagnetic fields and their gradients at different regions around the nanostructures, which are occupied by the oriented chiral molecules. Our calculations show that matched phases are as important as the intensities of the electromagnetic field at the molecular positions to enhance the CD signals of oriented molecular medium systems. This paper is organized as follows: Section 2A presents the effective medium theory of anisotropic chiral molecules based on the semi-classical theory of multipole oscillators induced by electromagnetic fields and the gradients of the electric field. Section 2B shows the method of

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 32

modeling wave propagation in anisotropic chiral medium. Section 3A is devoted to analyze the CD effects for the symmetric molecule distributions around nanostructures. Section 3B analyzes the CD effects of the asymmetric molecule distributions around nanostructures. Finally, a summary and the main conclusions of this research will be drawn.

II. Theory and method. A. Effective medium theory for the oriented chiral molecules with electric quadrupole transitions Although macroscopic chiral media are composed of microscopic chiral molecules, they are usually regarded as effective media when electromagnetic waves transport in them, which can be described by Maxwell’s equations with the following chiral constitutive relations: 35-36 D = ε 0ε% E + κ% B B = µ0 ( H − κ% T E ),

(1)

where the superscript ‘T’ indicates transposition; ε 0 and µ0 denote the vacuum permittivity and permeability, respectively. ε% is the relative permittivity dyadic and κ% is the magnetoelectric dyadic, a key parameter describing the chiral property of the system composed of oriented molecules. In the case of harmonic incident wave, κ% and ε% are given by the multipole theory purposed by Graham and Raab,37 and can be expressed as: (ε% )αβ = ε bδ αβ + (α% mac )αβ 1 (κ% )αβ = (η% mac )αβ + i ωΛ βγδ (ξ% mac )γδα . 2

(2)

ACS Paragon Plus Environment

4

Page 5 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Here Λ βαγ and δ αβ are the Levi-Civita symbol and the Kronecker delta operator; ε b is the relative permittivity of the embedded medium, ω is the angular frequency of the incident wave; (α% mac )αβ , (η% mac )αβ

and (ξ% mac )αβγ denote the corresponding matrix elements of the electric dipole-electric

dipole polarizability tensor, electric dipole-magnetic dipole polarizability tensor and electric dipole-electric quadrupole polarizability tensor of macroscopic chiral medium, respectively. In the limit of low density of molecules, the interactions among molecular dipoles can be ignored. Hence, the macroscopic polarizability tensors can be expressed by the polarizability tensors of the single molecule as: 21, 24 (α% mac )αβ = n0 (α% mic )αβ (η% mac )αβ = n0 (η% mic )αβ (ξ% mac )αβγ = n0 (ξ% mic )αβγ ,

(3)

where n0 is the molecular density; (α% mic )αβ , (η% mic )αβ and (ξ% mic )αβγ are the corresponding matrix elements of the electric dipole-electric dipole polarizability tensor, electric dipole-magnetic dipole polarizability tensor and electric dipole-electric quadrupole polarizability tensor for a single chiral molecule, respectively. These microscopic polarizability tensors can be derived from the following master equation for quantum states of the single molecule:

h

∂ρnm = i〈 n | [ρˆ , Hˆ 0 + VˆI ] | m〉 − (Γˆ ⋅ ρˆ )nm ∂t

(4)

where ρˆ is the density matrix, ρ nm (n (m)=1, 2) is the corresponding matrix element, Γˆ is the operator of the relaxation term which describes the phenomenological damping. The operator Hˆ describes the internal energy of the molecule.

VˆI

0

is the light-molecule interaction operator, which

can be written as:

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

) 1 ) ) VI = −µˆ.E − m.B − Q : ∇E 2

with

E = Re[ E0 e − iωt ]

and

B = Re[ B0 e − iωt ]

Page 6 of 32

(5)

being the incident electric and magnetic fields. µˆ , mˆ and Qˆ

are the molecular electric dipole, magnetic dipole and electric quadrupole operators (the primitive form38), respectively. By using the rotating-wave approximation in the linear regime, we can obtain the following solution from Eq. (4) as:

ρ 21 = Re[

With

1 VˆI21 = ( − µ 21 ⋅ E − m 21 ⋅ B − Q% 21 : ∇E ) 2

VˆI

21

(ω − ω0 ) h + iΓ 21

, where µ , 21

m21

e

− iωt

(6)

]

and

Q21

are the corresponding matrix

elements of molecular electric dipole, magnetic dipole and electric quadrupole moments. The parameter ω is the angular frequency of the molecular transition. Then, we can write the electric 0

dipole, magnetic dipole and electric quadrupole moments of the single molecule as: µ = µ12 ( ρ 21 + ρ12 ) m = m12 ρ 21 + m 21 ρ12

(7)

Q% = Q%12 ( ρ 21 + ρ12 ).

Inserting Eqs. (5) and (6) into Eq. (7), we obtain: 1 2

µ = α% mic E + η% mic B + ξ% mic (∇ E) m =-η% micT E Q% =ξ% micT E,

(8)

with

ACS Paragon Plus Environment

6

Page 7 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(α% mic )αβ = − ( (η% mic )αβ = − (

1 1 − )( µ12 ⊗ µ 21 )αβ hω − hω 0 + iΓ12 hω + hω 0 + i Γ 21 1 hω − hω 0 + iΓ12

(ξ% mic )αβγ = − (

1 hω − h ω0 + i Γ12

+

1 hω + hω 0 + iΓ 21



)( µ12 ⊗ m 21 )αβ

1 hω + hω 0 + iΓ 21

(9)

)( µ 21 )α (Q% 21 ) βγ

From Eqs. (2), (3) and (9), the general forms for the response tensors of anisotropic chiral medium with the electric quadrupole transitions can be obtained.

B. The method for modeling wave propagation in the anisotropic chiral medium A number of algorithms have been developed to calculate transport properties of the wave in the media.39-40 One of them is the finite element method (FEM) to solve the full Maxwell's equations. Over the last decades, numerical modeling algorithms are now being built into commercial or open-source software. COMSOL Multiphysics software,41 which is based on solution of second order differential equations for the electric field, is a numerical simulator based on the FEM. Recently, a rigorous finite element method that takes both the complex plasmonic structure and isotropic chiral medium into account has just been presented by using a modified COMSOL RF simulation environment.24 But an effective numerical approach to study the oriented molecular systems has not been elaborated. Here, we consider the molecular orientation-effects and the non-ignorable quadrupolar term in light-molecule interactions. The oriented molecule materials are regarded as anisotropic chiral media. The Maxwell’s equations for wave propagations in these chiral media can be expressed as: ⋅ ∇ × ∇ × E = ω 2 µ0 ε 0 ε% E + iωµ0κ% ( ∇ × E ) + iωµ 0 (∇ × κ% T ) ⋅ E − iωµ 0κ% T ( E∇ ∇ ). ×

(10)

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 32

⋅ ⋅ Here A% B% = Aij ei e j Bks ek es = Aij Bks ( ei × es )(e j ⋅ ek ) . If we assume that the magnetoelectric dyadic is ×

×

space-independent, the above equation can be reduced to: % 0 (∇ × E ) ∇ × ∇ × E = ω 2 µ 0 ε 0 ε% E + iωκµ − iωµ 0 [(κ% jy ∇ z E j − κ% jz ∇ y E j ) e x + (κ% jz ∇ x E j − κ% jx ∇ z E j ) e y + (κ% jx ∇ y E j − κ% jy ∇ x E j ) e z ].

(11)

In order to model wave propagations in the anisotropic chiral medium, we need to replace the wave equation for usual media by that of chiral media in Eq. 11.24, 30 Meanwhile, all parameters should be redefined through E according to the chiral constitutive relations described by Eq. 1. For example, the magnetic field vector should be expressed as:

H=

(∇ × E − i µ0ωκ% T E ) . −iωµ0

(12)

In addition, boundary conditions for Eq. 11 are not trivial, although they follow from the usual boundary conditions with: n × ( H 2 − H1 ) = 0 n × ( E2 − E1 ) = 0.

(13)

Combing Eq. 13 with the redefined magnetic field vector in Eq. 12, the surface current on the interface can be derived as:

n×(

B1

µ0



B2

µ0

) = n × (κ% T E ) = J b ,

(14)

where n is a vector normal to the boundary of the chiral medium (we only consider one kind of chiral medium in our system). What we need to emphasize is that the boundary conditions are indeed trivial (see in Eq. 13) for the bianisotropic media.42 There are other methods to simulate the chiral medium in Comsol. In fact, it is possible to leave the weak form of the wave equation

ACS Paragon Plus Environment

8

Page 9 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

in Comsol unchanged, and only to change the constitutive equations in “Equation View” of the “Wave Equation”. Comparing the present case with that of the isotropic chiral medium,24, 30 the response scalars have been replaced by tensors with the non-ignorable quadrupolar term in lightmolecule interactions. In addition, the chirality parameter κ% and dielectric tensors ε% for the oriented chiral molecular medium are derived from the first-principles calculations, shown in Section 2A. Both the ED-MD and ED-EQ interactions are considered during the process of solving the master equation for quantum states of the single oriented molecule. Consequently, the molecular electric quadrupolar contribution to the CD signal, which was always neglected in the past studies17, 19, 20-22, 24, can be calculated based on our method. The proposed method can be used to investigate the CD effect of the systems composed of more complex nanostructures and oriented chiral molecular medium of arbitrary shape with both ED-MD and ED-EQ interactions.

III. NUMERICAL RESULTS AND DISCUSSION. A. CD effect of symmetric molecule distributions around nanostructures. We first consider the oriented molecular systems without nanostructures. The schematics of the systems with different orientations of the molecular electric dipole moments are shown in the insets of Figure. 1. The chiral medium strips, which are arranged periodically with period p=320 nm, are infinite in the z direction, 0.75 nm and 10 nm along x- and y-axis (the size of the chiral medium is suited for the nanostructures discussed below). Two strips in each period are 80 nm apart. Throughout this work, the convenient parameters for the molecular electric dipole (ED) and magnetic dipole (MD) moments are

µ12 =| e | r12

and

m12

=

0.5i | e | r0 ω0 r12

molecular electric quadrupole (EQ) moments can be written as:

Q21,ij

. The form of the =

0.5 | e | rq 2 (i ≠ j )

and

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Q21,ii



| e | rq 2

Page 10 of 32

, where the subscripts i and j represent the coordinate components x, y and z,

respectively, and

|e|

is the absolute value of the electric charge. The form of the molecular EQ

moments will be changed with the rotation of the molecular ED moments. In this condition, minus sign is only used for the diagonal term of the EQ operators ( Q

21

,ii

) when the orientations of

molecular ED moments ( µ ) are along the i-axis. In the following calculations, we choose 12

parameters:19

r12 = 0.2 nm

,

r0 = 0.005 nm

molecular transition is taken as

,

rq = 0.025 nm

λ0 = 2π c / ω 0 = 300 nm

relative permittivity of the embedded medium is

and

Γ 21 =0.3 eV

; The wavelength of the

; The molecule density is

ε b =1 .

n0 =(2 nm)

−3

and the

Figure 1(a) and (b) show the numerical

results for the difference in molecular absorptivity (1-R-T, R and T are the reflectance and transmittance of the system) for left- and right- handed circularly polarized photons when the orientations of molecular ED moments are along x- and z-axis, respectively. Based on the numerical validation, we find that the influences of higher diffraction orders on the CD signal are negligible with the wavelength of the incident light larger than 200 nm in all of the following systems. We can see that if we neglect the molecular EQ transitions, the CD signal will decrease by increasing the angle ( θ , shown in Figure 1) between the molecular electric and magnetic dipole moments (black and green lines correspond to θ

= 0 deg

and θ

= 45 deg

) and equals to

zero with θ = 90 deg (blue line). Here, the molecular MD moments are rotated in x-y and x-z planes, respectively. As for the CD signals resulting from the molecular ED-EQ interaction (red lines), obviously, the molecular EQ transitions contribute in the same order as the MD on the whole CD for the oriented molecular systems.43 However, when the orientations of the molecular ED moments are parallel to the wave vector for the incident light, the CD signals nearly vanish as shown in Figure 1(c). On this occasion, the absorption of the molecular system becomes rather

ACS Paragon Plus Environment

10

Page 11 of 32

small, which results the corresponding CD signals can also be neglected. For comparison, we also plot the CD signal for the isotropic chiral medium with orange lines in Figure 1. It can be seen that the oriented molecular chiral medium may exhibit larger optical activity than the isotropic molecular system.

-8

θ=0deg only ED-MD θ=45deg only ED-MD θ=90deg only ED-MD Οnly ED-EQ

isotropy

∆Absorptivity

4.0x10

0.0 200

300

(a) 4.0x10

(b)

500

600 -8

∆Absorptivity

4.0x10

0.0 200

400

Wavelength(nm)

-8

∆Absorptivity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

300

400

500

Wavelength(nm)

600

0.0 200

(c)

300

400

500

Wavelength(nm)

600

Figure 1. (Color online) The difference in absorptivity of left- and right- handed circularly polarized photons excitation for the molecular systems that the chiral medium strips, arranged periodically with period p = 320 nm , are infinite in the z direction and 0.75 nm (10 nm) along x(y)- axis. Two strips in each period are

80 nm apart. The wave vector of the incident light is along the y-axis (identical with following systems). The orientations of ED moments are along (a) x-axis, (b) z-axis and (c) y-axis. The black, green and blue lines correspond to the CD signals, without EQ term, with the angles between the molecular ED and MD moments equaling to 0deg, 45deg and 90deg, respectively. The red lines represent the result of the EQ contribution. The orange lines show the solution of the isotropic chiral medium.

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 32

Some metallic nanostructures present a unique property, known as localized surface plasmon resonances (SPRs), under resonance excitations by external fields.44-45 The excitation of SPRs produces intense electric field enhancement in close proximity to the surfaces of metallic nanostructures, which promotes the light-matter interactions for the molecules located in the near field of nanostructures, leading to various plasmon-enhanced spectroscopies. In the following, we will investigate CD spectra of a few complex systems formed by the oriented chiral molecules and metallic nanostructures. As shown in the insets of Figure 2, the molecule-metal hybrids are arranged in an array of period p=320 nm. The chiral medium strips are infinite in the z direction, 0.75 nm and 10 nm along x- and y-axis, respectively. As for the gold nanoribbon (infinite along z-axis), the sizes in x and y directions are 80 nm and 10 nm. For the dielectric functions of gold, Johnson’s data were adopted.46 The orientations of the molecular ED moments are presented in the inset of Figure 2. Comparing with the case of bare chiral molecules (Figure 1), we found that the plasmon-induced CD signal appears in these complex systems when the local plasmonic resonance is excited (610 nm). As shown in Figure 2(a)-(c), we found that both the molecular ED-MD (black line) and ED-EQ (red line) interactions have a significant role on the plasmon-induced CD signal with different orientations of molecular ED moments. Meanwhile, when the molecules are achiral, the corresponding CD signals become 0 (blue line). Moreover, we need to emphasize that the above results do not depend on the orientations of the molecular MD moments. The averaged CD spectrum induced by the MD transitions (black line in Figure. 2(d)) is equal to the one when chiral media are considered to be isotropic (orange line in Figure. 2(d)). And the averaged CD signal resulting from the molecular EQ transitions is disappeared (red line in Figure. 2(d)), which supports the conclusion that EQ contribution to the differential absorption averages to zero in isotropic samples.6, 24

ACS Paragon Plus Environment

12

Page 13 of 32

-7

-7

1.50x10

1.50x10 80nm

10nm

0.75nm

y

x

ux

0.00 0deg only ED-MD only ED-EQ 90deg only ED-MD

200

(a)

400

600

800

1000

Wavelength(nm)

∆Absorbtivity

∆Absorbtivity

320nm

⊗ uz 0.00 0deg only ED-MD only ED-EQ 90deg only ED-MD

200

(b)

400

600

800

1000

Wavelength(nm) -7

0deg only ED-MD only ED-EQ 90deg only ED-MD

0.00

uy

-7

-1.50x10

(c)

200

400

600

800

Wavelength(nm)

1.50x10

∆Absorbtivity

∆Absorbtivity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

1000

r u Random 0.00 Averaged for ED-MD Isotropy Averaged for ED-EQ

200

(d)

400

600

800

1000

Wavelength(nm)

Figure 2. (Color online) The difference in the absorptivity of left- and right- handed circularly polarized photons excitation for the molecular systems formed by the oriented molecules and gold nanoribbons. The orientations of the molecular ED moments are along (a) x-axis (b) z-axis and (c) y-axis. The black and blue lines correspond to the CD signal, without EQ term, with the angle between the molecular ED and MD moments equaling to 0deg and 90deg, respectively. The red lines represent the result for the EQ contribution. (d) The comparison between the averaged CD signals, the ED-MD interaction (black) and ED-EQ interaction (red), and the CD effect of isotropic chiral medium system (orange).

As shown in the above results, the molecular ED-EQ interaction can play an equally important role as EQ-MD interaction on the plasmon-induced CD effects of oriented molecular systems. Next, we will demonstrate that the molecular EQ contribution can play a crucial role on

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 32

the plasmon-induced CD signal when chiral medium strips are put in the nanogaps. The corresponding schematics are presented in the insets of Figure 3. The dimensions of the gold nanoribbons are the same to the models as described in Figure 2. And the lengths along x- and yaxis of the chiral medium strips are 1.5 nm and 10 nm, respectively. The gap size (x direction) is 3 nm, and the period is also chosen to be 320 nm. It is shown clearly that the contributions of the molecular ED-MD interaction play a key role in the plasmon-induced circular dichroism spectra (the resonance is excited at 700 nm), and the contribution of the quadrupolar term almost disappears when the orientations of molecular ED moments are along x-axis, shown in Figure 3. However, if we alter the orientations of the molecular ED moments along z or y-axis, the situation is completely opposite. In this case, the molecular EQ transition plays a crucial role in the plasmon-induced CD spectra, this can be seen clearly in Figure 3(b) and (c). On the other hand, the averaged CD signals display the same behaviors as the above single nanoribbon system, shown in Figure 3(d), although the EQ transitions may play a more important role than the MD term on the CD spectra for the oriented molecules.

ACS Paragon Plus Environment

14

Page 15 of 32

-7

-7

1.0x10 ux 3nm

1.5nm

80nm y

320nm

x

0.0 0deg only ED-MD only ED-EQ 90deg only ED-MD

200

(a)

400

600

800

Wavelength(nm)

∆Absorbtivity

∆Absorbtivity

1.0x10

0.0

(b) 1.0x10

∆Absorbtivity

uy 0.0 0deg only ED-MD only ED-EQ 90deg only ED-MD

200

400

600

800

Wavelength(nm)

1000

-7

1.0x10

(c)

uz

0deg only ED-MD only ED-EQ 90deg only ED-MD

200

1000

-7



∆Absorbtivity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

400

600

800

Wavelength(nm)

1000

r u

Random

0.0

Averaged for ED-MD Isotropy Averaged for ED-EQ

200

(d)

400

600

800

Wavelength(nm)

1000

Figure 3. (Color online) The difference in the absorptivity of left- and right- handed circularly polarized photons excitation for the molecular systems formed by the oriented molecules and gap nanoribbons.

From the above results, we can see that the plasmon-induced CD is sensitive to the orientations of the molecules. These obvious oriented phenomena are mainly due to the remarkable anisotropy for the electromagnetic environment at the molecular positions produced by the metallic nanostructures, which has been shown in Figure 4. Figure 4(a)-(c) show the averaged intensities of the electric fields, magnetic fields and gradients of electric field within the volume Vblue (blue areas in the inset of Figure 4(a)), occupied by the chiral medium strips (show in Figure 2), as functions of the wavelength of the incident left-hand circular polarized light, respectively. The power of the incident light is taken as 1W/m in each period. For example,

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the averaged value for the x-polarized electric field can be expressed as: =

Page 16 of 32

∫ | E | dV / V x

blue

. In

Vblue

Figure 4(a) and (b), the black, red and blue lines represent the polarizations of the electromagnetic field along x, y and z directions, respectively. And for Figure 4(c), these lines illustrate the gradients of the electric field along x, y and z-axis. It is seen clearly that both the fields and their gradients exhibit the characteristics of anisotropy. As shown in Figure 4(a), we find that the averaged intensity of the electric field with the polarization direction along the xaxis is 10.8 and 9.8 times larger than that along y- and z-axis, respectively, when the local plasmonic resonance was excited (610 nm). The anisotropy of the electric fields at the positions of the chiral medium causes directional-selectivity for the ED moments during the process of light-molecules interaction. Similarly, the anisotropies for the magnetic fields (shown in Figure 4(b)) and gradients of electric field (shown in Figure 4(c)) cause directional-selectivity for the MD and EQ moments, respectively. Such a phenomenon related to the molecular orientations can be seen more clearly in the nanogap system. Figure 4(d)-(f) show the averaged intensities of the electric fields, magnetic fields and gradients of the electric field within the volume Vblue (blue areas in the inset of Figure 4(d)), occupied by the chiral medium strips (show in Figure 3), as functions of the wavelength of the incident left-hand circular polarized light, respectively. For comparative purposes, the power of the incident light is also chosen to be 1 W/m in each period. In such a case, the averaged intensity of the x-polarized electric field is 84.8 and 32.0 times larger than that polarized along y- and z-axis, respectively, when the local plasmonic resonance happens (700 nm). And comparing with the case of the single nanoribbon system, the averaged intensity of the x-polarized electric field is also amplified nearly 3-fold. From the above results, we can see that the value of x-polarized field is much larger than y- and z-polarized fields at molecular positions in the gap nanoribbon system. However, the magnetic fields and gradients of

ACS Paragon Plus Environment

16

Page 17 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

electric field have not been enhanced simultaneously compared to the single nanoribbon system, shown in Figure 4(d) and (e). The gradients of the electric field along x- and y-axis are only one third of the corresponding values of the single nanoribbon system. Consequently, we can speculate that the plasmonic CD signals induced by ED-MD and ED-EQ interactions are both dominated by the x-polarized electric fields in the nanogap system. As illustrated in Figure 3, the ED-MD interaction can induced large CD signals only when the orientations of the ED moments are along x-axis. Because only in this case, the coupling of the x-polarized electric field to EDMD transitions can occur. As for the case of ED-EQ interaction, the large CD signals can be induced when the orientations of the ED moments are along y- or z-axis with the coupling of the x-polarized electric field to ED-EQ transitions existing. These phenomena are obviously different from the cases of the single nanoribbon system that the CD effect produced by different orientations of molecules are all significant, shown in Figure 2. It can be explained by the fact that both the electric fields and gradients of electric field all have non-ignorable values in the single nanoribbon system, shown in Figure 4(a)-(c). Finally, we need to emphasize that the comparison between the results derived from the two nanostructures, shown in Figures 2 and 3, is not to distinguish the values of the CD signals. Our purpose is to disclose the physical origin of the phenomena that the dipole- and quadrupole-based CD signals have an important relationship with the molecular orientations.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

6

14

|Ex| |Ey| |Ez|

(a)

1000

Vblue

6

1.20x10

x

0

0.00 200

(d)

2

40

400 600 800 Wavelength[nm]

nanoribbon

(b)

1000

0 200

(c)

1x10

40

1000

400 600 800 Wavelength[nm]

120 80

|Ex| |Ey| |Ez| 400 600 800 Wavelength[nm]

200

|Hi|[A/m]

nanoribbon

|Hx| |Hy| |Hz|

0 200

(e)

14

2

0.00 200

80

|dE/di|[V/m ]

nanoribbon

|dE/di|[V/m ]

1x10

x

0

|Ei|[V/m]

120

Vblue

Vblue

|Hi|[A/m]

1.20x10

|Ei|[V/m]

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 32

|Hx| |Hy| |Hz|

400 600 800 Wavelength[nm]

1000

0 200

(f)

|dE/dx| |dE/dy| |dE/dz|

400

600

800

1000

Wavelength[nm] |dE/dx| |dE/dy| |dE/dz|

400

600

800

1000

Wavelength[nm]

Figure 4. (Color online) The averaged intensity for the (a) electric fields, (b) magnetic fields and (c) gradients of the electric field within the volume Vblue (blue areas in the inset of Figure 4(a)), occupied by the chiral medium strips, as a function of the wavelength of the periodic single-nanoribbons system. (d)-(f) plot the corresponding components for the periodic gap nanoribbons system.

B. CD effect of asymmetric molecule distributions around nanostructures. The above investigations only focus on the symmetric distributions of molecules around nanostructures. In fact, when the symmetry of the molecule distributions is broken, more prominent phenomenon can appear. In the following, we will investigate the CD of asymmetric distributions of the oriented molecular systems, which are presented in the insets of Figure 5. The chiral molecular media are only put at one side of the gold nanoribbons that the dimensions along x- and y-axis are 1.5 nm and 10 nm, respectively. As shown in Figure 5(a)-(f), the dipolebased CD spectra (black lines) are nearly identical with the corresponding symmetric systems (black lines in Figure 2(a)-(c)) when the orientations of the molecular ED and MD moments are

ACS Paragon Plus Environment

18

Page 19 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

paralleled. What may come as a surprise is that there are still non-ignorable CD signals (green lines) even the molecules are achiral, when the orientations of the molecular ED moments are along x- or y-axis, presented in Figure 5(a)-(d). Especially, when the molecular ED and MD moments are parallel to the y- and x-axis, the corresponding dipole-based CD peak (green lines in Figure 5(c)-(d)) is amplified tremendously comparing with the case that molecules were symmetrically distributed (black line in Figure 2(c)). However, when we change the directions of the molecular ED moments along the z-axis, shown in Figure 5(e) and (f), the ED-MD induced CD effects disappear if the molecules are achiral. In addition, the ED-EQ induced CD peaks (red lines) is nearly hundreds-fold over the corresponding symmetric system (red lines in Figure 2(b) and (c)), when the orientations of the molecular ED moments are along y- or z-axis, presented in Figure 5(c)-(f). On the other hand, in many case, the plasmon-induced CD signals are entirely different, when the molecules are positioned at different sides of the gold nanoribbon. Moreover, CD signals with asymmetric molecule distributions also depend on the orientations of the molecular MD moments, show in Figure 5. The above phenomena resulting from molecular asymmetric distributions can also exist in the nanogap system with the molecules clinging to the left or right side in the gap region, shown in the insets of Figure 6. From Figure 6(a)-(d), we can see that there are non-ignorable CD signals (green lines), even the molecules are achiral. However, when the directions of the molecular ED moments are along the z-axis, the corresponding CD signal will disappear with the molecules being achiral, shown in Figure 6(e) and (f). In this condition, the ED-EQ transition dominates the whole plasmon-induced CD signal, and the CD-peak is nearly hundreds-fold over the corresponding symmetric system (red line in Figure 3(b)).

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

-8

-8

8.0x10

8.0x10

∆Absorptivity

320nm 80nm

ux

1.5nm

10nm

y

x 0deg only ED-MD 90deg only ED-MD mz=0 90deg only ED-MD my=0 only ED-EQ

0.0 200

(a)

400

600

800

1000

ux 0.0 200

Wavelength(nm)

(b)

400

600

800

1000

Wavelength(nm)

-6

∆Absorptivity

1.0x10 0deg only ED-MD 90deg only ED-MD mz=0 90deg only ED-MD mx=0 only ED-EQ

0.0

uy 0.0

uy

-6

-1.0x10 200

(c)

400

600

800

1000

200

Wavelength(nm)

(d)

400

600

800

1000

Wavelength(nm)

-7

8.0x10

0.0

∆Absorptivity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 32

0deg only ED-MD 90deg only ED-MD mx=0 90deg only ED-MD my=0 only ED-EQ

⊗ uz 0.0

⊗ uz

-7

200

(e)

400

600

800

Wavelength(nm)

1000

-8.0x10 200

(f)

400

600

800

1000

Wavelength(nm)

Figure 5. (Color online) The difference of the absorptivity of left- and right- handed circularly polarized photons excitation with the asymmetric molecule distributions formed by the oriented molecules and gold nanoribbons. (a), (c) and (e) left side, (b), (d) and (f) right side exist molecules with the orientations of the molecular ED moments along x-axis, y-axis and z-axis, respectively. The black lines correspond to the CD signals, without EQ term, with the angle between the molecular ED and MD moments equaling to 0deg. The green and blue lines correspond to the CD signals of achiral molecules, which are not take the ED-EQ

ACS Paragon Plus Environment

20

Page 21 of 32

interaction into account, with different components of MD moments equaling to 0. The red lines represent the results for the EQ contributions.

-8

-8

5.0x10

∆Absorptivity

ux 3nm

1.5nm

80nm

-8

400

600

800

1000

-5.0x10 200

Wavelength(nm)

400

600

800

1000

Wavelength(nm)

(b) -7

8.0x10

-7

8.0x10

∆Absorptivity

ux

x

0deg only ED-MD 0.0 90deg only ED-MD mz=0 90deg only ED-MD my=0 only ED-EQ

-8

(a)

5.0x10

y

320nm

0.0

-5.0x10 200

uy 0deg only ED-MD 90deg only ED-MD mz=0 0.0 90deg only ED-MD mx=0 only ED-EQ

0.0

uy -7

-7

-8.0x10 200

(c)

400

600

800

1000

-8.0x10 200

400

-7

600

800

1000

Wavelength(nm)

(d)

Wavelength(nm)

-7

5.0x10

∆Absorptivity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

5.0x10

0deg only ED-MD 90deg only ED-MD my=0 90deg only ED-MD mx=0 only ED-EQ

0.0

⊗ uz 0.0

⊗ uz -7

-5.0x10 200

(e)

400

600

800

Wavelength(nm)

1000

-7

-5.0x10 200

(f)

400

600

800

1000

Wavelength(nm)

Figure 6. (Color online) The difference of the absorptivity of left- and right- handed circularly polarized photons excitation with the asymmetric molecule distributions formed by the oriented molecules and gap nanoribbons.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 32

To disclose the physical origin of the phenomena happened with asymmetric molecular distributions, we plot the differences of averaged phases for near-fields between different sides of the nanostructures as functions of the wavelength for the incident left-hand circular polarized light, shown in Figure 7. Figure 7(a)-(c) show differences of averaged phases of electric fields, magnetic fields and gradients of electric field between Vred (red areas in the inset of Figure 7(a)) and Vblue (blue areas in the inset of Figure 7(a)) areas around the single nanoribbon system as functions of the wavelength for the incident left-hand circular polarized light, respectively (the difference of averaged phase of x-polarized electric fields between Vred and Vblue areas can be expressed as: < ∆Arg ( Ex ) >=



Vred

Abs[ Arg ( Ex )]dV / Vred −



Abs[ Arg ( Ex )]dV / Vblue , with the variation range

Vblue

of the corresponding phases being −π → π ). In Figure 7(a) and (b), the black, red and blue lines represent the polarization of the electromagnetic fields along x, y and z directions, respectively. In Figure 7(c) these lines illustrate the electric field gradients along x, y and z-axis. We find that the difference of averaged phases between both sides of the single nanoribbon are nonzero for the y-polarized electric and magnetic fields and gradients of electric field along x-axis. Due to these non-uniform phase distributions for the electromagnetic components, which are coupled to the molecular ED-MD or ED-EQ transition, the oriented chiral molecular media located at different sides of the nanostructures may induced CD signals with the opposite sign, even mutual offset (green lines in Figure 5(c) and (d)). For example, comparing with the case that molecules are symmetrically distributed (Figure 2(b)), asymmetric molecular distributions (Figure 5(e) and (f)) may enlarge the CD-peak (red lines) induced by ED-EQ interaction over hundreds-fold. This can be explained by the condition that the gradients of the electric field along x-axis are coupled with molecular ED-EQ transitions when the orientation of the molecular electric dipole moment

ACS Paragon Plus Environment

22

Page 23 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

is along z-axis. This will result in the CD signals produced by the chiral molecular media, which lie in different sides of the single nanoribbon, possessing opposite sign (due to non-uniform sign of the averaged phase of the gradients of the electric field along x-axis in different sides of the single nanoribbon). Similar results also exist in the gap nanoribbon systems. Figure 7(d)-(f) show differences of averaged phases of electric fields, magnetic fields and gradients of electric field between Vred (red areas in the inset of Figure 7(d)) and Vblue (blue areas in the inset of Figure 7(d)) areas around the gap nanoribbon system as functions of the wavelength of the incident lefthand circular polarized light, respectively. We find that the nonzero components of differences of averaged phases are identical with the single nanoribbon system. Hence, similar phenomena happened with asymmetric molecular distributions for single nanoribbon and gap nanoribbon systems. Our calculations show that matching the phase for the electromagnetic components, which are coupling with molecular ED-MD or ED-EQ interaction, is significant for the CD enhancement of oriented chiral molecules system. Finally, we need to clarify that the phases discussed here are not the local phase difference between the electric and magnetic fields. What we concerned is the spatial distributions of the phases for the specific electromagnetic fields, which are coupled with ED-MD or ED-EQ interaction of the oriented chiral molecules to exhibit chiroptical effects.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

-2

Vblue

x

0

-4 200

(a)

nanoribbon 400 600 800 Wavelength[nm]

-4 1000 (b) 200 4

nanoribbon 0 nanoribbon 0

x

∆Arg(Ex) ∆Arg(Ey) ∆Arg(Ez)

200

400 600 800 Wavelength[nm]

0 -2

Vred Vblue

4

∆Arg(dE/di)

Vred

4

∆Arg(Hx) ∆Arg(Hy) ∆Arg(Hz)

1000

2

400 600 800 Wavelength[nm]

1000

∆Arg(Hx) ∆Arg(Hy) ∆Arg(Hz)

0 -2

(e)

-4 200

400 600 800 Wavelength[nm]

1000

2

∆Arg(dE/dx) ∆Arg(dE/dy) ∆Arg(dE/dz)

0 -2 -4 200

(c)

4 ∆Arg(dE/di)

0

2 ∆Arg(Hi)

∆Arg(Ei)

2

(d)

4

∆Arg(Ex) ∆Arg(Ey) ∆Arg(Ez)

∆Arg(Hi)

4

∆Arg(Ei)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 32

2

400 600 800 Wavelength[nm]

1000

∆Arg(dE/dx) ∆Arg(dE/dy) ∆Arg(dE/dz)

0 -2 -4 200

(f)

400 600 800 Wavelength[nm]

1000

Figure 7. (Color online) The difference of the averaged phases for the (a) electric fields, (b) magnetic fields and (c) gradients of the electric fields between Vred (red areas in the inset of Figure 7(a)) and Vblue (blue areas in the inset of Figure 7(a)) areas around the single nanoribbon system as a function of wavelength of the incident left-hand circular polarized light. (d)-(f) plot the corresponding components for the periodic gap nanoribbons system.

Conclusion. We have used a rigorous finite element method to study the CD effect of various systems with spatial symmetric and asymmetric molecular distributions around nanostructures. The orientation-sensitive CD effects are investigated with non-ignorable EQ term in light-molecule interaction, which was always omitted under the assumption of a large collections of molecules random oriented. Although the CD signal resulting from the molecular quadrupole moments in our system can play a key role in many cases, the relatively low enhancements of near field

ACS Paragon Plus Environment

24

Page 25 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

gradients proposed by our one-dimensional wires systems still cannot amplify the CD above the experimental detection limit (the same as the case of the isotropic chiral medium).24 By designing nanostructures with steeper field gradients, the much larger CD signal resulting from molecular quadrupole transition may be achieved. Furthermore, we have demonstrated that both the quadrupole- and dipole-based CD signals can be improved greatly by laying the chiral molecules asymmetrically around the nanostructures to match the phases for the electromagnetic components, which are coupled with molecular ED-MD or ED-EQ interaction, at the molecule positions. This method can lead to the uniform sign of the CD produced by the oriented chiral molecules, which are located at different spatial positions near the nanostructures. We believe that these findings would be helpful for realizing ultrasensitive probing of chiral information for oriented molecules by plasmon-based nanotechnology.

AUTHOR INFORMATION Corresponding Author * E-mail: [email protected] Notes The authors declare no competing financial interest. Acknowledgment This work was supported by the National Natural Science Foundation of China (Grant No. 11274042) and the National Key Basic Research Special Foundation of China under Grant 2013CB632704. Supporting Information

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 32

Details of analytical derivations for the anisotropic wave equation and the check of our numerical method are given. This material is available free of charge via the Internet at http://pubs.acs.org.

REFERENCES (1) Circular dichroism and the conformational analysis of biomolecules; Fasman, G. D., Ed.; Plenum: New York, 1996. (2) Wagniere, G. H. On chirality and the universal asymmetry: reflections on image and mirror image; VHCA with Wiley-VCH: Zurich, 2007; Busch, K. W. and Busch, M. A. (eds) Chiral Analysis (Elsevier, 2006).

(3) Nguyen, L. A.; He, H.; Pham-Huy, C. Chiral Drugs: an Overview. Int. J. Biomed. Sci.: IJBS 2006, 2, 85−100.

(4) Berova, N.; Polavarapu, P. L.; Nakanishi, K.; Woody, R. W. Comprehensive Chiroptical Spectroscopy; John Wiley and Sons Inc.: New York, 2012.

(5) Patterson, D.; Schnell, M.; Doyle, J. M. Enantiomer-specific detection of chiral molecules via microwave spectroscopy. Nature 2013, 497, 475-478. (6) Tang, Y.; Cohen, A. E. Optical Chirality and Its Interaction with Matter. Phys. Rev. Lett. 2010, 104, 163901. (7) Tang, Y.; Cohen, A. E. Enhanced Enantioselectivity in Excitation of Chiral Molecules by Superchiral Light. Science 2011, 332, 333−336.

ACS Paragon Plus Environment

26

Page 27 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(8) Maoz, B. M.; Chaikin, Y.; Tesler, A. B.; Elli, O. B.; Fan, Z.; Govorov, A. O.; Markovich, G. Amplification of Chiroptical Activity of Chiral Biomolecules by Surface Plasmons. Nano Lett. 2013, 13, 1203−1209.

(9) Lu, F.; Tian, Y.; Liu, M.; Su, D.; Zhang, H.; Govorov, A. O.; Gang, O. Discrete Nanocubes as Plasmonic Reporters of Molecular Chirality. Nano Lett. 2013, 13, 3145−3151. (10) Kuzyk, A.; Schreiber, R.; Fan, Z.; Pardatscher, G.; Roller, E.-M.; Hogele, A.; Simmel, F. C.; Govorov, A. O. ; Liedl, T. DNA-based self-assembly of chiral plasmonic nanostructures with tailored optical response. Nature (London) 2012, 483, 311-314. (11) Ma, W.; Kuang, H.; Xu, L.; Ding, L.; Xu, C.; Wang, L.; Kotov, N. A. Attomolar. DNA Detection with Chiral Nanorod Assemblies. Nat. Commun. 2013, 4, 2689. (12) Xu, Z.; Xu, L.; Liz Marzan, L. M.; Ma, W.; Kotov, N. A.; Wang, ́ L.; Kuang, H.; Xu, C. Sensitive Detection of Silver Ions Based on Chiroplasmonic Assemblies of Nanoparticles. Adv. Opt. Mater. 2013, 1, 626−630.

(13) Zhao, Y.; Xu, L.; Ma, W.; Wang, L.; Kuang, H.; Xu, C.; Kotov, N. A. Shell-Engineered Chiroplasmonic Assemblies of Nanoparticles for Zeptomolar DNA Detection. Nano Lett. 2014, 14, 3908−3913. (14) Xu, Z.; Xu, L.; Zhu, Y.; Ma, W.; Kuang, H.; Wang, L.; Xu, C. Chirality based sensor for bisphenol A Detection. Chem. Commun. 2012, 48, 5760-5762. (15) Li, Z.; Zhu, Z.; Liu, W.; Zhou, Y.; Han, B.; Gao, Y.; Tang, Z.; Reversible plasmonic circular dichroism of Au nanorod and DNA assemblies. J. Am. Chem. Soc. 2012, 134, 33223325.

ACS Paragon Plus Environment

27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 32

(16) Maoz, B. M.; Weegen, R.; van der Fan, Z.; Govorov, A. O.; Ellestad, G.; Berova, N.; Meijer, E. W.; Markovich, G. Plasmonic Chiroptical Response of Silver Nanoparticles Interacting with Chiral Supramolecular Assemblies. J. Am. Chem. Soc. 2012, 134, 17807−17813. (17) Govorov, A. O.; Fan, Z.; Hernandez, P.; Slocik, J. M.; Naik, R. R. Theory of Circular Dichroism of Nanomaterials Comprising Chiral Molecules and Nanocrystals: Plasmon Enhancement, Dipole Iteraction, and Dielectric Effects. Nano Lett. 2010, 10, 1374-1382. (18) Wu, T.; Ren, J.; Wang, R.; Zhang, X. Competition of Chiroptical Effect Caused by Nanostructure and Chiral Molecules. J. Phys. Chem. C 2014, 118, 20529−20537. (19) Zhang, H.; Govorov, A. O. Giant circular dichroism of a molecule in a region of strong plasmon resonances between two neighboring gold nanocrystals. Phys. Rev. B 2013, 87, 075410075417. (20) Govorov, A. O. Plasmon-Induced Circular Dichroism of a Chiral Molecule in the Vicinity of Metal Nanocrystals. Application to Various Geometries. J. Phys. Chem. C 2011, 115, 7914−7923. (21) Govorov, A. O.; Fan, Z. Theory of Chiral Plasmonic Nanostructures Comprising Metal Nanocrystals and Chiral Molecular Media. ChemPhysChem 2012, 13, 2551−2560. (22) Abdulrahman, N. A.; Fan, Z.; Tonooka, T.; Kelly, S. M.; Gadegaard, N.; Hendry, E.; Govorov, A. O.; Kadodwala, M. Induced Chirality Through Electromagnetic Coupling Between Chiral Molecular Layers and Plasmonic Nanostructures. Nano Lett. 2012, 12, 977−983.

ACS Paragon Plus Environment

28

Page 29 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(23) Garcia-Etxarri, A.; Dionne, J. A. Surface-Enhanced Circular Dichroism Spectroscopy Mediated by Nonchiral Nanoantennas. Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, 235409. (24) Nesterov, M. L.; Yin, X. H.; Schäferling, M.; Giessen, H.; Weiss, T. The Role of Plasmon-Generated Near Fields for Enhanced Circular Dichroism Spectroscopy. ACS. Photonics 2016, 3, 578-583. (25) Davis, T. J.; Gomez, D. E. Interaction of Localized Surface Plasmons with Chiral Molecules. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 90, 235424. (26) Hendry, E.; Carpy, T.; Johnston, J.; Popland, M.; Mikhaylovskiy, R. V.; Lapthorn, A. J.; Kelly, S. M.; Barron, L. D.; Gadegaard, N.; Kadodwala, M. Ultrasensitive detection and characterization of biomolecules using superchiral field. Nat. Nanotechnol. 2010, 5, 783-787. (27) Tullius, R.; Karimullah, A. S.; Rodier, M.; Fitzpatrick, B.; Gadegaard, N.; Barron, L. D.; Rotello, V. M.; Cooke, G.; Lapthorn, A.; Kadodwala, M. “Superchiral” spectroscopy: detection of protein higher order hierarchical structure with chiral plasmonic. J. Am. Chem. Soc. 2015, 137, 8380-8383. (28) Wang, R. Y.; Wang, P.; Liu, Y.; Zhao, W.; Zhai, D.; Hong, X.; Ji, Y.; Wu, X.; Wang, F.; Zhang, D.; Zhang, W.; Liu, R.; Zhang, X. Experimental observation of giant chiroptical amplification of small chiral molecules by gold nanosphere clusters. J. Phys. Chem. C 2014,118, 9690−9695. (29) Levi-Belenkova, T; Govorov, A. O.; Markovich, G. Orientation-sensitive peptide-induced plasmonic circular dichroism in silver nanocubes. J. Phys. Chem. C 2016,120, 12751−12756.

ACS Paragon Plus Environment

29

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 32

(30) Klimov, V. V.; Zabkov, I. V.; Pavlov, A. A.; Guzatov, D. V. Eigen Oscillations of a Chiral Sphere and Their Influence on Radiation of Chiral Molecules. Opt. Express 2014, 22, 18564−18578. (31) Schäferling, M.; Engheta, N.; Giessen, H.; Weiss, T. Reducing the complexity: Enantioselective chiral near-fields by diagonal slit and mirror configuration. ACS. Photonics 2016, 3, 1076-1084. (32) Efrima, S. Raman Optical Activity of Molecules Adsorbed on Metal Surfaces: Theory. J. Chem. Phys. 1985, 83, 1356−1362.

(33) Yang, N.; Cohen, A. E. Local geometry of electromagnetic fields and its role in molecular multipole transitions. J. Phys. Chem. B 2011,115, 5304−5311. (34) Wu, T.; Zhang, X. H.; Wang, R. Y.; Zhang, X. Strongly enhanced Raman optical activity in molecules by magnetic response of nanoparticles. J. Phys. Chem. C 2016, 118, 20529−20537. (35) Kong, J. A. Theory of Electromagnetic Waves; Wiley: New York, 1975. (36) Post, E. J. Formal Structure of Electromagnetics; North-Holland: Amsterdam, 1962. (37) Graham, E. B.; Raab, R. E. Covariant D and H fields for reflection from a magnetic anisotropic chiral medium. J. Opt. Soc. Am. A 1997, 14, 131−134. (38) Graham, E. B.; Pierrus, J.; Raab, R. E. Multipole moments and Maxwell’s equations. J. Phys. B 1992, 25, 4673−4684.

(39) Tai, C. T. Dyadic Green Functions in Electromagnetic Theory; IEEE Press: New York, 1993.

ACS Paragon Plus Environment

30

Page 31 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(40) Bohren, C. F.; Huffman, D. R. Absorption and Scattering of light by Small Particles; Wiley: Interscience, 1983. (41) www.comsol.com (42) Onishi, M.; Crabtree, K.; Chipman, R. A. Formulation of Rigorous Coupled-Wave Theory for Gratings in Bianisotropic Media. J. Opt. Soc. Am. A 2011, 28, 1747−1758. (43) Barron, L. D. Molecular Light Scattering and Optical Activity; 2nd Ed.; Cambridge University Press: New York, 2009. (44) Novotny, L.; Hecht, B. Principles of Nano-Optics; Cambridge University Press: U.K., 2006. (45) Barnes, W. L.; Dereux, A.; Ebbesen, T. W. Surface plasmon subwavelength optics. Nature 2003, 424, 824−830.

(46) Johnson, P. B.; Christy, R. W. Phys. Rev. B 1972, 6, 4370-4379.

ACS Paragon Plus Environment

31

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 32

TOC Graphic:

ACS Paragon Plus Environment

32