Analytical calculation of far infrared spectra of ice in terms of a molecular model

Wideband far infrared (FIR) spectra of complex permittivity e(p) of ice are calculated in terms of a simple analytical theory based on the method of dipolar autocorrelation functions. The molecular model represents a revision of the model recently presented for liquid water in Adv. Chem. Phys. 127 (2003) 65. A composite two-fractional model is proposed. The model is characterised by three phenomenological potential wells corresponding to the three FIR bands observed in ice. The first fraction comprises dipoles reorienting in a rather narrow and deep hat-like well; these dipoles generate the librational band centred at the frequency approximate to 880 cm(-1). The second fraction comprises elastically interacting particles; they generate two nearby bands placed around frequency 200 cm(-1). For description of one of these bands the harmonic oscillator (HO) model is used, in which translational oscillations of two charged molecules along the H-bond are considered. The other band is produced by the H-bond stretch, which governs hindered rotation of a rigid dipole. Such a motion and its dielectric response are described in terms of a new cut parabolic (CP) model applicable for any vibration amplitude. The composite hat-HO-CP model results in a smooth epsilon(nu) ice spectrum, which does not resemble the noise-like spectra of ice met in the known literature. The proposed theory satisfactorily agrees with the experimental ice spectrum measured at - 7 degrees C. The calculated longitudinal optic-transverse optic (LO-TO) splitting occurring at approximate to 250 cm(-1) qualitatively agrees with the measured data. (c) 2004 Elsevier B.V. All rights reserved.

Nonharmonic transverse vibration of the H-bonded molecules and the THz spectra in ice and water

A nonlinear equation of motion is found for the dimer comprising two charged H2O molecules. The THz dielectric response to nonharmonic vibration of a nonrigid dipole, forming the hydrogen bond (HB), is found in the direction transverse to this bond. An explicit expression is derived for the autocorrelator that governs the spectrum generated by transverse vibration (TV) of such a dipole. This expression is obtained by analytical solution of the truncated set of recurrence equations. The far infrared (FIR) spectra of ice at the temperature - 7 degrees C are calculated. The wideband, in the wavenumber (frequency) v range 0... 100.0 cm(-1), spectra are obtained for liquid water at room temperature and for supercooled water at -5.6 degrees C. All spectra are represented in terms of the complex permittivity epsilon(v) and the absorption coefficient alpha(v). The obtained analytical formula for epsilon comprises the term epsilon(perpendicular to) pertinent to the studied TV mechanism with three additional terms Delta epsilon(q), Delta epsilon(mu), and epsilon(or) arising, respectively, from: elastic harmonic vibration of charged molecules along the H-bond; elastic reorientation of HB permanent dipoles; and rather free libration of permanent dipoles in 'defects' of water/ice structure. The suggested TV-dielectric relaxation mechanism allows us: (a) to remove the THz 'deficit' of loss epsilon" inherent in previous theoretical studies; (b) to explain the THz loss and absorption spectra in supercooled (SC) water; and (c) to describe, in agreement with the experiment, the low- and high-frequency tails of the two bands of ice H2O located in the range 10...300 cm(-1). Specific THz dielectric properties of SC water are ascribed to association of water molecules, revealed in our study by transverse vibration of HB charged molecules. (C) 2006 Published by Elsevier B.V.

Three-dimensional photonic crystals based on macroporous silicon with modulated pore diameter

We report on the fabrication and optical characterization of a three-dimensional (3D) photonic crystal on the basis of macroporous silicon. The structure consists of a 2D array of air pores in silicon whose diameter is varied (modulated) periodically with depth. The bandstructure of the resulting 3D hexagonal photonic crystal is calculated and compared with transmission measurements. The described structure allows to adjust the dispersion relation along the pore axis almost independently from the dispersion relation in the plane perpendicular to the pore axis.

Optical bistability in nonlinear surface-plasmon polaritonic crystals

Nonlinear optical transmission through periodically nanostructured metal films (surface-plasmon polaritonic crystals) has been studied. The surface polaritonic crystals have been coated with a nonlinear polymer. The optical transmission of such nanostructures has been shown to depend on the control-light illumination conditions. The resonant transmission exhibits bistable behavior with the control-light intensity. The bistability is different at different resonant signal wavelengths and for different wavelengths of the control light. The effect is explained by the strong sensitivity of the surface-plasmon mode resonances at the signal wavelength to the surrounding dielectric environment and the electromagnetic field enhancement due to plasmonic excitations at the controlled light wavelengths.

Ferroelectricity and isotope effects in hydrogen-bonded KDP crystals

Based on an accurate first principles description of the energetics in H-bonded potassium-dihydrogen-phosphate crystals, we conduct a first study of nuclear quantum effects and of the changes brought about by deuteration. Tunneling is allowed only for clusters involving correlated protons and heavy ion displacements, the main effect of deuteration being a depletion of the proton probability density at the O-H-O bridge center, which in turn weakens its proton-mediated covalent bonding. The ensuing lattice expansion couples self-consistently with the proton off-centering, thus explaining both the giant isotope effect and its close connection with geometrical effects.

Interplay between proton ordering and ferroelectric polarization in H-bonded KDP-type crystals

The origin of ferroelectricity in KH2PO4 (KDP) is studied by first-principles electronic structure calculations. In the low-temperature phase, the collective off-centre ordering of the protons is accompanied by an electronic charge delocalization from the near and localization at the far oxygen within the O-H...O bonds. Electrostatic forces. then, push the K+ ions towards off-centre positions, and induce a macroscopic polarization. The analysis of the correlation between different geometrical and electronic quantities, in connection with experimental data. supports the idea that the role of tunnelling in isotopic effects is irrelevant. Instead, geometrical quantum effects appear to play a central role. (C) 2001 Elsevier Science B.V. All rights reserved.

First-principles study of ferroelectricity and isotope effects in H-bonded KH2PO4 crystals

By means of extensive first-principles calculations we studied the ferroelectric phase transition and the associated isotope effect in KH2PO4 (KDP). Our calculations revealed that the spontaneous polarization of the ferroelectric phase is due to electronic charge redistributions and ionic displacements which are a consequence of proton ordering, and not vice versa. The experimentally observed double-peaked proton distribution in the paraelectric phase cannot be explained by a dynamics of only protons. This requires, instead, collective displacements within clusters that include also the heavier ions. These tunneling clusters can explain the recent evidence of tunneling obtained from Compton scattering measurements. The sole effect of mass change upon deuteration is not sufficient to explain the huge isotope effect. Instead, we find that structural modifications deeply connected with the chemistry of the H bonds produce a feedback effect on tunneling that strongly enhances the phenomenon. The resulting influence of the geometric changes on the isotope effect agrees with experimental data from neutron scattering. Calculations under pressure allowed us to analyze the issue of universality in the disappearance of ferroelectricity upon compression. Compressing DKDP so that the distance between the two peaks in the deuteron distribution is the same as for protons in KDP, corresponds to a modification of the underlying double-well potential, which becomes 23 meV shallower. This energy difference is what is required to modify the O-O distance in such a way as to have the same distribution for protons and deuterons. At the high pressures required experimentally, the above feedback mechanism is crucial to explain the magnitude of the geometrical effect.

Local and semilocal density functional computations for crystals of 1-alkyl-3-methyl-imidazolium salts

The accuracy and reliability of popular density functional approximations for the compounds giving origin to room temperature ionic liquids have been assessed by computing the T=0 K crystal structure of several 1-alkyl-3-methyl-imidazolium salts. Two prototypical exchange-correlation approximations have been considered, i.e., the local density approximation (LDA) and one gradient corrected scheme [PBE-GGA, Phys. Rev. Lett. 77, 3865 (1996)]. Comparison with low-temperature x-ray diffraction data shows that the equilibrium volume predicted by either approximations is affected by large errors, nearly equal in magnitude (~10%), and of opposite sign. In both cases the error can be traced to a poor description of the intermolecular interactions, while the intramolecular structure is fairly well reproduced by LDA and PBE-GGA. The PBE-GGA optimization of atomic positions within the experimental unit cell provides results in good agreement with the x-ray structure. The correct system volume can also be restored by supplementing PBE-GGA with empirical dispersion terms reproducing the r-6 attractive tail of the van der Waals interactions.

Dispersing Light with Surface Plasmon Polaritonic Crystals

Spectral dispersion of light on a finite-size surface plasmon polaritonic (SPP) crystal has been studied. The angular wavelength separation of one or more orders of magnitude higher than in other state-of-the-art wavelength-splitting devices available to date has been demonstrated. The two-stage process is responsible for the dispersion value, which involves conversion of the incident light into SPP Bloch modes of a nanostructure followed by the SPP Bloch waves refraction at the SPP crystal boundary. The high spectral dispersion achievable in plasmonic devices may be useful for integrated high-resolution spectroscopy in nanophotonic, optical communication and lab-on-a-chip applications.