Theory of strongly correlated electron systems. I. Intersite coulomb interaction and the approximation of renormalized fermions in total energy calculations

The diagrammatic strong-coupling perturbation theory (SCPT) for correlated electron systems is developed for intersite Coulomb interaction and for a nonorthogonal basis set. The construction is based on iterations of exact closed equations for many - electron Green functions (GFs) for Hubbard operators in terms of functional derivatives with respect to external sources. The graphs, which do not contain the contributions from the fluctuations of the local population numbers of the ion states, play a special role: a one-to-one correspondence is found between the subset of such graphs for the many - electron GFs and the complete set of Feynman graphs of weak-coupling perturbation theory (WCPT) for single-electron GFs. This fact is used for formulation of the approximation of renormalized Fermions (ARF) in which the many-electron quasi-particles behave analogously to normal Fermions. Then, by analyzing: (a) Sham's equation, which connects the self-energy and the exchange- correlation potential in density functional theory (DFT); and (b) the Galitskii and Migdal expressions for the total energy, written within WCPT and within ARF SCPT, a way we suggest a method to improve the description of the systems with correlated electrons within the local density approximation (LDA) to DFT. The formulation, in terms of renormalized Fermions LIDA (RF LDA), is obtained by introducing the spectral weights of the many electron GFs into the definitions of the charge density, the overlap matrices, effective mixing and hopping matrix elements, into existing electronic structure codes, whereas the weights themselves have to be found from an additional set of equations. Compared with LDA+U and self-interaction correction (SIC) methods, RF LDA has the advantage of taking into account the transfer of spectral weights, and, when formulated in terms of GFs, also allows for consideration of excitations and nonzero temperature. Going beyond the ARF SCPT, as well as RF LIDA, and taking into account the fluctuations of ion population numbers would require writing completely new codes for ab initio calculations. The application of RF LDA for ab initio band structure calculations for rare earth metals is presented in part 11 of this study (this issue). (c) 2005 Wiley Periodicals, Inc.

Theory of strongly correlated electron systems. II. Including correlation effects into electronic structure calculations

We have previously shown that a division of the f-shell into two subsystems gives a better understanding of the cohesive properties as well the general behavior of lanthanide systems. In this article, we present numerical computations, using the suggested method. We show that the picture is consistent with most experimental data, e.g., the equilibrium volume and electronic structure in general. Compared with standard energy band calculations and calculations based on the self-interaction correction and LIDA + U, the f-(non-f)-mixing interaction is decreased by spectral weights of the many-body states of the f-ion. (c) 2005 Wiley Periodicals, Inc.

Complex numbers and symmetries in quantum mechanics, and a nonlinear superposition principle for Wigner functions

Complex numbers appear in the Hilbert space formulation of quantum mechanics, but not in the formulation in phase space. Quantum symmetries are described by complex, unitary or antiunitary operators defining ray representations in Hilbert space, whereas in phase space they are described by real, true representations. Equivalence of the formulations requires that the former representations can be obtained from the latter and vice versa. Examples are given. Equivalence of the two formulations also requires that complex superpositions of state vectors can be described in the phase space formulation, and it is shown that this leads to a nonlinear superposition principle for orthogonal, pure-state Wigner functions. It is concluded that the use of complex numbers in quantum mechanics can be regarded as a computational device to simplify calculations, as in all other applications of mathematics to physical phenomena.

First-principles quantum dynamics in interacting Bose gases II: stochastic gauges

First principles simulations of the quantum dynamics of interacting Bose gases using the stochastic gauge representation are analysed. In a companion paper, we showed how the positive-P representation can be applied to these problems using stochastic differential equations. That method, however, is limited by increased sampling error as time evolves. Here, we show how the sampling error can be greatly reduced and the simulation time significantly extended using stochastic gauges. In particular, local stochastic gauges (a subset) are investigated. Improvements are confirmed in numerical calculations of single-, double- and multi-mode systems in the weak-mode coupling regime. Convergence issues are investigated, including the recognition of two modes by which stochastic equations produced by phase-space methods in general can diverge: movable singularities and a noise-weight relationship. The example calculated here displays wave-like behaviour in spatial correlation functions propagating in a uniform 1D gas after a sudden change in the coupling constant. This could in principle be tested experimentally using Feshbach resonance methods.

First-principles quantum simulations of dissociation of molecular condensates: Atom correlations in momentum space

We investigate the quantum many-body dynamics of dissociation of a Bose-Einstein condensate of molecular dimers into pairs of constituent bosonic atoms and analyze the resulting atom-atom correlations. The quantum fields of both the molecules and atoms are simulated from first principles in three dimensions using the positive-P representation method. This allows us to provide an exact treatment of the molecular field depletion and s-wave scattering interactions between the particles, as well as to extend the analysis to nonuniform systems. In the simplest uniform case, we find that the major source of atom-atom decorrelation is atom-atom recombination which produces molecules outside the initially occupied condensate mode. The unwanted molecules are formed from dissociated atom pairs with nonopposite momenta. The net effect of this process-which becomes increasingly significant for dissociation durations corresponding to more than about 40% conversion-is to reduce the atom-atom correlations. In addition, for nonuniform systems we find that mode mixing due to inhomogeneity can result in further degradation of the correlation signal. We characterize the correlation strength via the degree of squeezing of particle number-difference fluctuations in a certain momentum-space volume and show that the correlation strength can be increased if the signals are binned into larger counting volumes.

Studies on Wilson loops, correlators and localization in supersymmetric quantum field theories

In this thesis we study at perturbative level correlation functions of Wilson loops (and local operators) and their relations to localization, integrability and other quantities of interest as the cusp anomalous dimension and the Bremsstrahlung function. First of all we consider a general class of 1/8 BPS Wilson loops and chiral primaries in N=4 Super Yang-Mills theory. We perform explicit two-loop computations, for some particular but still rather general configuration, that confirm the elegant results expected from localization procedure. We find notably full consistency with the multi-matrix model averages, obtained from 2D Yang-Mills theory on the sphere, when interacting diagrams do not cancel and contribute non-trivially to the final answer. We also discuss the near BPS expansion of the generalized cusp anomalous dimension with L units of R-charge. Integrability provides an exact solution, obtained by solving a general TBA equation in the appropriate limit: we propose here an alternative method based on supersymmetric localization. The basic idea is to relate the computation to the vacuum expectation value of certain 1/8 BPS Wilson loops with local operator insertions along the contour. Also these observables localize on a two-dimensional gauge theory on S^2, opening the possibility of exact calculations. As a test of our proposal, we reproduce the leading Luscher correction at weak coupling to the generalized cusp anomalous dimension. This result is also checked against a genuine Feynman diagram approach in N=4 super Yang-Mills theory. Finally we study the cusp anomalous dimension in N=6 ABJ(M) theory, identifying a scaling limit in which the ladder diagrams dominate. The resummation is encoded into a Bethe-Salpeter equation that is mapped to a Schroedinger problem, exactly solvable due to the surprising supersymmetry of the effective Hamiltonian. In the ABJ case the solution implies the diagonalization of the U(N) and U(M) building blocks, suggesting the existence of two independent cusp anomalous dimensions and an unexpected exponentation structure for the related Wilson loops.

A micro-opto-electro-mechanical system (MOEMS) for microstructure manipulation

Microstructure manipulation is a fundamental process to the study of biology and medicine, as well as to advance micro- and nano-system applications. Manipulation of microstructures has been achieved through various microgripper devices developed recently, which lead to advances in micromachine assembly, and single cell manipulation, among others. Only two kinds of integrated feedback have been demonstrated so far, force sensing and optical binary feedback. As a result, the physical, mechanical, optical, and chemical information about the microstructure under study must be extracted from macroscopic instrumentation, such as confocal fluorescence microscopy and Raman spectroscopy. In this research work, novel Micro-Opto-Electro-Mechanical-System (MOEMS) microgrippers are presented. These devices utilize flexible optical waveguides as gripping arms, which provide the physical means for grasping a microobject, while simultaneously enabling light to be delivered and collected. This unique capability allows extensive optical characterization of the structure being held such as transmission, reflection, or fluorescence. The microgrippers require external actuation which was accomplished by two methods: initially with a micrometer screw, and later with a piezoelectric actuator. Thanks to a novel actuation mechanism, the "fishbone", the gripping facets remain parallel within 1 degree. The design, simulation, fabrication, and characterization are systematically presented. The devices mechanical operation was verified by means of 3D finite element analysis simulations. Also, the optical performance and losses were simulated by the 3D-to-2D effective index (finite difference time domain FDTD) method as well as 3D Beam Propagation Method (3D-BPM). The microgrippers were designed to manipulate structures from submicron dimensions up to approximately 100 μm. The devices were implemented in SU-8 due to its suitable optical and mechanical properties. This work demonstrates two practical applications: the manipulation of single SKOV-3 human ovarian carcinoma cells, and the detection and identification of microparts tagged with a fluorescent "barcode" implemented with quantum dots. The novel devices presented open up new possibilities in the field of micromanipulation at the microscale, scalable to the nano-domain.

Ab initio quantum chemical studies on neutral-radical reactions of ethynyl (C2H) and cyano (CN) with unsaturated hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA's Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth's primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan's atmosphere, the ISM, and cold celestial bodies.<.

A Micro-Opto-Electro-Mechanical System (MOEMS) for Microstructure Manipulation

Microstructure manipulation is a fundamental process to the study of biology and medicine, as well as to advance micro- and nano-system applications. Manipulation of microstructures has been achieved through various microgripper devices developed recently, which lead to advances in micromachine assembly, and single cell manipulation, among others. Only two kinds of integrated feedback have been demonstrated so far, force sensing and optical binary feedback. As a result, the physical, mechanical, optical, and chemical information about the microstructure under study must be extracted from macroscopic instrumentation, such as confocal fluorescence microscopy and Raman spectroscopy. In this research work, novel Micro-Opto-Electro-Mechanical-System (MOEMS) microgrippers are presented. These devices utilize flexible optical waveguides as gripping arms, which provide the physical means for grasping a microobject, while simultaneously enabling light to be delivered and collected. This unique capability allows extensive optical characterization of the structure being held such as transmission, reflection, or fluorescence. The microgrippers require external actuation which was accomplished by two methods: initially with a micrometer screw, and later with a piezoelectric actuator. Thanks to a novel actuation mechanism, the “fishbone”, the gripping facets remain parallel within 1 degree. The design, simulation, fabrication, and characterization are systematically presented. The devices mechanical operation was verified by means of 3D finite element analysis simulations. Also, the optical performance and losses were simulated by the 3D-to-2D effective index (finite difference time domain FDTD) method as well as 3D Beam Propagation Method (3D-BPM). The microgrippers were designed to manipulate structures from submicron dimensions up to approximately 100 µm. The devices were implemented in SU-8 due to its suitable optical and mechanical properties. This work demonstrates two practical applications: the manipulation of single SKOV-3 human ovarian carcinoma cells, and the detection and identification of microparts tagged with a fluorescent “barcode” implemented with quantum dots. The novel devices presented open up new possibilities in the field of micromanipulation at the microscale, scalable to the nano-domain.

Quantum-limited estimation of continuous spontaneous localization

We apply the formalism of quantum estimation theory to extract information about potential collapse mechanisms of the continuous spontaneous localisation (CSL) form.
In order to estimate the strength with which the field responsible for the CSL mechanism couples to massive systems, we consider the optomechanical interaction
between a mechanical resonator and a cavity field. Our estimation strategy passes through the probing of either the state of the oscillator or that of the electromagnetic field that drives its motion. In particular, we concentrate on all-optical measurements, such as homodyne and heterodyne measurements.
We also compare the performances of such strategies with those of a spin-assisted optomechanical system, where the estimation of the CSL parameter is performed
through time-gated spin-like measurements.

Approximate ab initio calculations and the method of molecular fragments

A two stage approach to performing ab initio calculations on medium and large sized molecules is described. The first step is to perform SCF calculations on small molecules or molecular fragments using the OPIT Program. This employs a small basis set of spherical and p-type Gaussian functions. The Gaussian functions can be identified very closely with atomic cores, bond pairs, lone pairs, etc. The position and exponent of any of the Gaussian functions can be varied by OPIT to produce a small but fully optimised basis set. The second stage is the molecular fragments method. As an example of this, Gaussian exponents and distances are taken from an OPIT calculation on ethylene and used unchanged in a single SCF calculation on benzene. Approximate ab initio calculations of this type give much useful information and are often preferable to semi-empirical approaches, since the nature of the approximations involved is much better defined.

Tomographic readout of an opto-mechanical interferometer

The quantum state of light changes its nature when being reflected off a mechanical oscillator due to the latter's susceptibility to radiation pressure. As a result, a coherent state can transform into a squeezed state and can get entangled with the motion of the oscillator. Full information of the state of light can only be gathered by a tomographic measurement. Here we demonstrate a tomographic interferometer readout by measuring arbitrary quadratures of the light field exiting a Michelson-Sagnac interferometer that contains a thermally excited high-quality silicon nitride membrane. A readout noise of 1.9 x 10(-16) mHz(-1/2) around the membrane's fundamental oscillation mode at 133 kHz has been achieved, going below the peak value of the standard quantum limit by a factor of 8.2 (9 dB). The readout noise was entirely dominated by shot noise in a rather broad frequency range around the mechanical resonance.

Finite difference calculations of permeability in large domains in a wide porosity range.

Determining effective hydraulic, thermal, mechanical and electrical properties of porous materials by means of classical physical experiments is often time-consuming and expensive. Thus, accurate numerical calculations of material properties are of increasing interest in geophysical, manufacturing, bio-mechanical and environmental applications, among other fields. Characteristic material properties (e.g. intrinsic permeability, thermal conductivity and elastic moduli) depend on morphological details on the porescale such as shape and size of pores and pore throats or cracks. To obtain reliable predictions of these properties it is necessary to perform numerical analyses of sufficiently large unit cells. Such representative volume elements require optimized numerical simulation techniques. Current state-of-the-art simulation tools to calculate effective permeabilities of porous materials are based on various methods, e.g. lattice Boltzmann, finite volumes or explicit jump Stokes methods. All approaches still have limitations in the maximum size of the simulation domain. In response to these deficits of the well-established methods we propose an efficient and reliable numerical method which allows to calculate intrinsic permeabilities directly from voxel-based data obtained from 3D imaging techniques like X-ray microtomography. We present a modelling framework based on a parallel finite differences solver, allowing the calculation of large domains with relative low computing requirements (i.e. desktop computers). The presented method is validated in a diverse selection of materials, obtaining accurate results for a large range of porosities, wider than the ranges previously reported. Ongoing work includes the estimation of other effective properties of porous media.

Ab Initio Quantum Chemical Studies on Neutral-Radical Reactions of Ethynyl (C2H) and Cyano (CN) with Unsaturated Hydrocarbons

An Ab Initio/RRKM study of the reaction mechanism and product branching ratios of neutral-radical ethynyl (C2H) and cyano (CN) radical species with unsaturated hydrocarbons is performed. The reactions studied apply to cold conditions such as planetary atmospheres including Titan, the Interstellar Medium (ISM), icy bodies and molecular clouds. The reactions of C2H and CN additions to gaseous unsaturated hydrocarbons are an active area of study. NASA’s Cassini/Huygens mission found a high concentration of C2H and CN from photolysis of ethyne (C2H2) and hydrogen cyanide (HCN), respectively, in the organic haze layers of the atmosphere of Titan. The reactions involved in the atmospheric chemistry of Titan lead to a vast array of larger, more complex intermediates and products and may also serve as a chemical model of Earth’s primordial atmospheric conditions. The C2H and CN additions are rapid and exothermic, and often occur barrierlessly to various carbon sites of unsaturated hydrocarbons. The reaction mechanism is proposed on the basis of the resulting potential energy surface (PES) that includes all the possible intermediates and transition states that can occur, and all the products that lie on the surface. The B3LYP/6-311g(d,p) level of theory is employed to determine optimized electronic structures, moments of inertia, vibrational frequencies, and zero-point energy. They are followed by single point higher-level CCSD(T)/cc-vtz calculations, including extrapolations to complete basis sets (CBS) of the reactants and products. A microcanonical RRKM study predicts single-collision (zero-pressure limit) rate constants of all reaction paths on the potential energy surface, which is then used to compute the branching ratios of the products that result. These theoretical calculations are conducted either jointly or in parallel to experimental work to elucidate the chemical composition of Titan’s atmosphere, the ISM, and cold celestial bodies.

Quantum Monte Carlo study of effective masses in a polarized unitary Fermi gas

Ultracold gases provide an ideal platform for quantum simulations of many-body systems. Here we are interested in a particular system which has been the focus of most experimental and theoretical works on ultracold fermionic gases: the unitary Fermi gas. In this work we study with Quantum Monte Carlo simulations a two-component gas of fermionic atoms at zero temperature in the unitary regime. Specifically, we are interested in studying how the effective masses for the quasi-particles of the two components of the Fermi liquid evolve as the polarization is progressively reduced from full to lower values. A recent theoretical work, based on alternative diagrammatic methods, has indeed suggested that such effective masses should diverge at a critical polarization. To independently verify such predictions, we perform Variational Monte Carlo (VMC) calculations of the energy based on Jastrow-Slater wavefunctions after adding or subtracting a particle with a given momentum to a full Fermi sphere. In this way, we determine the quasi-particle dispersions, from which we extract the effective masses for different polarizations. The resulting effective masses turn out to be quite close to the non-interacting values, even though some evidence of an increase for the effective mass of the minority component appears close to the predicted value for the critical polarization. Preliminary results obtained for the majority component with the Fixed-node Diffusion Monte Carlo (DMC) method seem to indicate that DMC could lead to an increase of the effective masses in comparison with the VMC results. Finally, we point out further improvements of the trial wave-function and boundary conditions that would be necessary in future simulations to draw definite conclusions on the effective masses of the polarized unitary Fermi gas.