Strong second harmonic generation in SiC, ZnO, GaN two-dimensional hexagonal crystals from first-principles many-body calculations

The second harmonic generation (SHG) intensity spectrum of SiC, ZnO, GaN two-dimensional hexagonal crystals is calculated by using a real-time first-principles approach based on Green's function theory [Attaccalite et al., Phys. Rev. B: Condens. Matter Mater. Phys. 2013 88, 235113]. This approach allows one to go beyond the independent particle description used in standard first-principles nonlinear optics calculations by including quasiparticle corrections (by means of the GW approximation), crystal local field effects and excitonic effects. Our results show that the SHG spectra obtained using the latter approach differ significantly from their independent particle counterparts. In particular they show strong excitonic resonances at which the SHG intensity is about two times stronger than within the independent particle approximation. All the systems studied (whose stabilities have been predicted theoretically) are transparent and at the same time exhibit a remarkable SHG intensity in the range of frequencies at which Ti:sapphire and Nd:YAG lasers operate; thus they can be of interest for nanoscale nonlinear frequency conversion devices. Specifically the SHG intensity at 800 nm (1.55 eV) ranges from about 40-80 pm V(-1) in ZnO and GaN to 0.6 nm V(-1) in SiC. The latter value in particular is 1 order of magnitude larger than values in standard nonlinear crystals.

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.

Single-ionization of helium at Ti:sapphire wavelengths: rates and scaling laws

We present a numerical and theoretical study of intense-field single-electron ionization of helium at 390 nm and 780 nm. Accurate ionization rates (over an intensity range of (0.175-34) X10^14 W/ cm^2 at 390 nm, and (0.275 - 14.4) X 10^14 W /cm^2 at 780 nm) are obtained from full-dimensionality integrations of the time-dependent helium-laser Schroedinger equation. We show that the power law of lowest order perturbation theory, modified with a ponderomotive-shifted ionization potential, is capable of modelling the ionization rates over an intensity range that extends up to two orders of magnitude higher than that applicable to perturbation theory alone. Writing the modified perturbation theory in terms of scaled wavelength and intensity variables, we obtain to first approximation a single ionization law for both the 390 nm and 780 nm cases. To model the data in the high intensity limit as well as in the low, a new function is introduced for the rate. This function has, in part, a resemblance to that derived from tunnelling theory but, importantly, retains the correct frequency-dependence and scaling behaviour derived from the perturbative-like models at lower intensities. Comparison with the predictions of classical ADK tunnelling theory confirms that ADK performs poorly in the frequency and intensity domain treated here.

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.

Existence of multibreathers in one-dimensional Debye crystals

The existence of highly localized multisite oscillatory structures (discrete multibreathers) in a nonlinear Klein-Gordon chain which is characterized by an inverse dispersion law is proven and their linear stability is investigated. The results are applied in the description of vertical (transverse, off-plane) dust grain motion in dusty plasma crystals, by taking into account the lattice discreteness and the sheath electric and/or magnetic field nonlinearity. Explicit values from experimental plasma discharge experiments are considered. The possibility for the occurrence of multibreathers associated with vertical charged dust grain motion in strongly coupled dusty plasmas (dust crystals) is thus established. From a fundamental point of view, this study aims at providing a rigorous investigation of the existence of intrinsic localized modes in Debye crystals and/or dusty plasma crystals and, in fact, suggesting those lattices as model systems for the study of fundamental crystal properties.

Nonlinear dynamics and solitons in Debye bi-crystals

The nonlinear dynamics of longitudinal dust lattice waves propagating in a dusty plasma bi-crystal is investigated. A “diatomic”-like one-dimensional dust lattice configuration is considered, consisting of two distinct dust grain species with different charges and masses. Two different frequency dispersion modes are obtained in the linear limit, namely, an optical and an acoustic wave dispersion branch. Nonlinear solitary wave solutions are shown to exist in both branches, by considering the continuum limit for lattice excitations in different nonlinear potential regimes. For this purpose, a generalized Boussinesq and an extended Korteweg de Vries equation is derived, for the acoustic mode excitations, and their exact soliton solutions are provided and compared. For the optic mode, a nonlinear Schrödinger-type equation is obtained, which is shown to possess bright- (dark-) type envelope soliton solutions in the long (short, respectively) wavelength range. Optic-type longitudinal wavepackets are shown to be generally unstable in the continuum limit, though this is shown not to be the rule in the general (discrete) case.

Nanoscale ferroelectrics machined from single crystals

This paper summarises some of the most recent work that has been done on nanoscale ferroelectrics as a result of a joint collaborative research effort involving groups in Queen's University Belfast, the University of Cambridge and the University of St. Andrews. Attempts have been made to observe fundamental effects of reduced size, and increasing morphological complexity, on ferroelectric behaviour by studying the functional response and domain characteristics in nanoscale single crystal material, whose size and morphology have been defined by Focused Ion Beam (FIB) patterning. This approach to nanoshape fabrication has allowed the following broad statements to be made: (i) in single crystal BaTiO3 sheets, permittivity and phase transition behaviour is not altered from that of bulk material down to a thickness of similar to 75 nm; (ii) in single crystal BaTiO3 sheets and nanowires changes in observed domain morphologies are consistent with large scale continuum modeling.

Thin film capacitor cut from single crystals using focused ion beam milling

The focused ion beam microscope (FIB) has been used to fabricate thin parallel-sided ferroelectric capacitors from single crystals of BaTiO3 and SrTiO3. A series of nano-sized capacitors ranging in thickness from similar to660 nm to similar to300 nm were made. Cross-sectional high resolution transmission electron microscopy (HRTEM) revealed that during capacitor fabrication, the FIB rendered around 20 nm of dielectric at the electrode-dielectric interface amorphous, associated with local gallium impregnation. Such a region would act electrically in series with the single crystal and would presumably have a considerable negative influence on the dielectric properties. However, thermal annealing prior to gold electrodes deposition was found to fully recover the single crystal capacitors and homogenise the gallium profile. The dielectric testing of the STO ultra-thin single crystal capacitors was performed yielding a room temperature dielectric constant of similar to300, as is the case in bulk. Therefore, there was no evidence of a collapse in dielectric constant associated with thin film dimensions.