Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Shortcut to Adiabaticity in the Lipkin-Meshkov-Glick Model

We study transitionless quantum driving in an infinite-range many-body system described by the Lipkin-Meshkov-Glick model. Despite the correlation length being always infinite the closing of the gap at the critical point makes the driving Hamiltonian of increasing complexity also in this case. To this aim we develop a hybrid strategy combining a shortcut to adiabaticity and optimal control that allows us to achieve remarkably good performance in suppressing the defect production across the phase transition.

Macroscopic Quantum Resonators (MAQRO): 2015 update

Do the laws of quantum physics still hold for macroscopic objects - this is at the heart of Schrödinger’s cat paradox - or do gravitation or yet unknown effects set a limit for massive particles? What is the fundamental relation between quantum physics and gravity? Ground-based experiments addressing these questions may soon face limitations due to limited free-fall times and the quality of vacuum and microgravity. The proposed mission Macroscopic Quantum Resonators (MAQRO) may overcome these limitations and allow addressing such fundamental questions. MAQRO harnesses recent developments in quantum optomechanics, high-mass matter-wave interferometry as well as state-of-the-art space technology to push macroscopic quantum experiments towards their ultimate performance limits and to open new horizons for applying quantum technology in space. The main scientific goal is to probe the vastly unexplored ‘quantum-classical’ transition for increasingly massive objects, testing the predictions of quantum theory for objects in a size and mass regime unachievable in ground-based experiments. The hardware will largely be based on available space technology. Here, we present the MAQRO proposal submitted in response to the 4th Cosmic Vision call for a medium-sized mission (M4) in 2014 of the European Space Agency (ESA) with a possible launch in 2025, and we review the progress with respect to the original MAQRO proposal for the 3rd Cosmic Vision call for a medium-sized mission (M3) in 2010. In particular, the updated proposal overcomes several critical issues of the original proposal by relying on established experimental techniques from high-mass matter-wave interferometry and by introducing novel ideas for particle loading and manipulation. Moreover, the mission design was improved to better fulfill the stringent environmental requirements for macroscopic quantum experiments.

Nonequilibrium critical scaling from quantum thermodynamics

The emerging field of quantum thermodynamics is contributing important results and insights into archetypal many-body problems, including quantum phase transitions. Still, the question whether out-of-equilibrium quantities, such as fluctuations of work, exhibit critical scaling after a sudden quench in a closed system has remained elusive. Here, we take a novel approach to the problem by studying a quench across an impurity quantum critical point. By performing density matrix renormalization group computations on the two-impurity Kondo model, we are able to establish that the irreversible work produced in a quench exhibits finite-size scaling at quantum criticality. This scaling faithfully predicts the equilibrium critical exponents for the crossover length and the order parameter of the model, and, moreover, implies a new exponent for the rescaled irreversible work. By connecting the irreversible work to the two-impurity spin correlation function, our findings can be tested experimentally.

Electron energy-loss near-edge shape as a probe to investigate the stabilization of yttria-stabilized zirconia

The electron energy-loss near-edge structure (ELNES) at the O K edge has been studied in yttria-stabilized zirconia (YSZ). The electronic structure of YSZ for compositions between 3 and 15 mol % Y2O3 has been computed using a pseudopotential-based technique to calculate the local relaxations near the O vacancies. The results showed phase transition from the tetragonal to cubic YSZ at 10 mol % of Y2O3, reproducing experimental observations. Using the relaxed defect geometry, calculation of the ELNES was carried out using the full-potential linear muffin-tin orbital method. The results show very good agreement with the experimental O K-edge signal, demonstrating the power of using ELNES to probe the stabilization mechanism in doped metal oxides.

Electron impact ionization of He above the threshold

Near-threshold ionization of He has been studied by using a uniform semiclassical wavefunction for the two outgoing electrons in the final channel. The quantum mechanical transition amplitude for the direct and exchange scattering derived earlier by using the Kohn variational principle has been used to calculate the triple differential cross sections. Contributions from singlets and triplets are critically examined near the threshold for coplanar asymmetric geometry with equal energy sharing by the two outgoing electrons. It is found that in general the tripler contribution is much smaller compared to its singlet counterpart. However, at unequal scattering angles such as theta (1) = 60 degrees, theta (2) = 120 degrees the smaller peaks in the triplet contribution enhance both primary and secondary TDCS peaks. Significant improvements of the primary peak in the TDCS are obtained for the singlet results both in symmetric and asymmetric geometry indicating the need to treat the classical action variables without any approximation. Convergence of these cross sections are also achieved against the higher partial waves. Present results are compared with absolute and relative measurements of Rosel et al (1992 Phys. Rev. A 46 2539) and Selles et al (1987 J. Phys. B. At. Mel. Phys. 20 5195) respectively.

Density functional and Monte Carlo studies of sulfur. I. Structure and bonding in S-n rings and chains (n=2-18) (vol 118, pg 9257, 2003)

Density functional calculations have been performed for ring isomers of sulfur with up to 18 atoms, and for chains with up to ten atoms. There are many isomers of both types, and the calculations predict the existence of new forms. Larger rings and chains are very flexible, with numerous local energy minima. Apart from a small, but consistent overestimate in the bond lengths, the results reproduce experimental structures where known. Calculations are also performed on the energy surfaces of S8 rings, on the interaction between a pair of such rings, and the reaction between one S8 ring and the triplet diradical S8 chain. The results for potential energies, vibrational frequencies, and reaction mechanisms in sulfur rings and chains provide essential ingredients for Monte Carlo simulations of the liquid–liquid phase transition. The results of these simulations will be presented in Part II.

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.

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.

Understanding thickness effects in thin film capacitors

Pulsed Laser Deposition (PLD) was used to make Au/(Ba0.5Sr0.5)TiO3/(La0.5Sr0.5) CoO3/MgO thin film capacitor structures. Functional properties were studied with changing BST thickness from similar to1265 nm to similar to63 nm. The dielectric constant was found to decrease, and migration of T-m (the temperature at which the dielectric constant is maximum) to lower temperatures occurred as thickness was reduced. Curie-Weiss plots of the as-obtained dielectric data, indicated that the Curie temperature was also systemmatically progressively depressed. Further, fitting to expressions previously used to describe diffuse phase transitions suggested increased diffuseness in transformation behaviour as film thickness decreased. This paper discusses the care needed in interpreting the observations given above. We make particular distinction between the apparent Curie-temperature derived from Curie-Weiss plots of as-measured data, and the inherent Curie temperature determined after correction for the interfacial capacitance. We demonstrate that while the apparent Curie temperature decreases as thickness decreases, the inherent Curie temperature is thickness independent. Thickness-invariant phase transition behaviour is confirmed from analysis of polarisation loops, and from examination of the temperature dependence of the loss-tangent. We particularly note that correction of data for interfacial capacitance does not alter the position of T-m. We must therefore conclude that the position of T-m is not related simply to phase transformation behaviour in BST thin films.

Origin of Ferroelastic domains in Free-Standing Single Crystal Ferroelectric Films

The origin of the unusual 90 degrees ferroelectric/ferroelastic domains, consistently observed in recent studies on mesoscale and nanoscale free-standing single crystals of BaTiO3 [Schilling , Phys. Rev. B 74, 024115 (2006); Schilling , Nano Lett. 7, 3787 (2007)], has been considered. A model has been developed which postulates that the domains form as a response to elastic stress induced by a surface layer which does not undergo the paraelectric-ferroelectric cubic-tetragonal phase transition. This model was found to accurately account for the changes in domain periodicity as a function of size that had been observed experimentally. The physical origin of the surface layer might readily be associated with patterning damage, seen in experiment; however, when all evidence of physical damage is removed from the BaTiO3 surfaces by thermal annealing, the domain configuration remains practically unchanged. This suggests a more intrinsic origin, such as the increased importance of surface tension at small dimensions. The effect of surface tension is also shown to be proportional to the difference in hardness between the surface and the interior of the ferroelectric. The present model for surface-tension induced twinning should also be relevant for finely grained or core-shell structured ceramics.

The Influence of Point Defects and Inhomogeneous Strain on the Functional Behaviour of Thin film Ferroelectrics

An attempt has been made to unequivocally identify the influence that inhomogeneous strain fields, surrounding point defects, have on the functional properties of thin film ferroelectrics. Single crystal thin film lamellae of BaTiO3 have been integrated into capacitor structures, and the functional differences between those annealed in oxygen and those annealed in nitrogen have been mapped. Key features, such as the change in the paraelectric-ferroelectric phase transition from first to second order were noted and found to be consistent with mean field modeling predictions for the effects of inhomogeneous strain. Switching characteristics appeared to be unaffected, suggesting that point defects have a low efficacy in domain wall pinning.

Domains in Ferroelectric Nanodots

Almost free-standing single crystal mesoscale and nanoscale dots of ferroelectric BaTiO3 have been made by direct focused ion beam patterning of bulk single crystal material. The domain structures which appear in these single crystal dots, after cooling through the Curie temperature, were observed to form into quadrants, with each quadrant consisting of fine 90° stripe domains. The reason that these rather complex domain configurations form is uncertain, but we consider and discuss three possibilities for their genesis: first, that the quadrant features initially form to facilitate field-closure, but then develop 90° shape compensating stripe domains in order to accommodate disclination stresses; second, that they are the result of the impingement of domain packets which nucleate at the sidewalls of the dots forming “Forsbergh” patterns (essentially the result of phase transition kinetics); and third, that 90° domains form to conserve the shape of the nanodot as it is cooled through the Curie temperature but arrange into quadrant packets in order to minimize the energy associated with uncompensated surface charges (thus representing an equilibrium state). While the third model is the preferred one, we note that the second and third models are not mutually exclusive.

Temperature-Driven Mixing-Demixing Behavior of Binary Mixtures of the Ionic Liquid Choline Bis(trifluoromethylsulfonyl)imide and Water

The ionic liquid (2-hydroxyethylammonium)trimethylammonium) bis(trifluoromethylsulfonyl)imide (choline bistriflimide) was obtained as a supercooled liquid at room temperature (melting point = 30 degrees C). Crystals of choline bistriflimide suitable for structure determination were grown from the melt in situ on the X-ray diffractometer. The choline cation adopts a folded conformation, whereas the bistriflimide anion exhibits a transoid conformation. The choline cation and the bistriflimide anion are held together by hydrogen bonds between the hydroxyl proton and a sulfonyl oxygen atom. This hydrogen bonding is of importance for the temperature-dependent solubility proper-ties of the ionic liquid. Choline bistriflimide is not miscible with water at room temperature, but forms one phase with water at temperatures above 72 degrees C (equals upper critical solution temperature). H-1 NMR studies show that the hydrogen bonds between the choline cation and the bistriflimide anion are substantially weakened above this temperature. The thermophysical properties of water-choline bistriflimide binary mixtures were furthermore studied by a photopyroelectric technique and by adiabatic scanning calorimetry (ASC). By photothermal analysis, besides highly accurate values for the thermal conductivity and effusivity of choline bistriflimide at 30 degrees C, the detailed temperature dependence of both the thermal conductivity and effusivity of the upper and lower part of a critical water-choline bistriflimide mixture in the neighborhood of the mixing-demixing phase transition could be determined with high resolution and accuracy. Together with high resolution ASC data for the heat capacity, experimental values were obtained for the critical exponents alpha and beta, and for the critical amplitude ratio G(+)/G(-). These three values were found to be consistent with theoretical expectations for a three dimensional Ising-type of critical behavior of binary liquid mixtures.

Quantitative mapping of aqueous microfluidic temperature with sub-degree resolution using fluorescence lifetime imaging microscopy

The use of a water-soluble, thermo-responsive polymer as a highly sensitive fluorescence-lifetime probe of microfluidic temperature is demonstrated. The fluorescence lifetime of poly(N-isopropylacrylamide) labelled with a benzofurazan fluorophore is shown to have a steep dependence on temperature around the polymer phase transition and the photophysical origin of this response is established. The use of this unusual fluorescent probe in conjunction with fluorescence lifetime imaging microscopy (FLIM) enables the spatial variation of temperature in a microfluidic device to be mapped, on the micron scale, with a resolution of less than 0.1 degrees C. This represents an increase in temperature resolution of an order of magnitude over that achieved previously by FLIM of temperature-sensitive dyes