-
Understanding Discrepancies of Wavefunction Theories for Large Molecules
Authors:
Tobias Schäfer,
Andreas Irmler,
Alejandro Gallo,
Andreas Grüneis
Abstract:
Quantum mechanical many-electron calculations can predict properties of atoms, molecules and even complex materials. The employed computational methods play a quintessential role in many scientifically and technologically relevant research fields. However, a question of paramount importance is whether approximations aimed at reducing the computational complexity for solving the many-electron Schrö…
▽ More
Quantum mechanical many-electron calculations can predict properties of atoms, molecules and even complex materials. The employed computational methods play a quintessential role in many scientifically and technologically relevant research fields. However, a question of paramount importance is whether approximations aimed at reducing the computational complexity for solving the many-electron Schrödinger equation, are accurate enough. Here, we investigate recently reported discrepancies of noncovalent interaction energies for large molecules predicted by two of the most widely-trusted many-electron theories: diffusion quantum Monte Carlo and coupled-cluster theory. We are able to unequivocally pin down the source of the puzzling discrepancies and present modifications to widely-used coupled-cluster methods needed for more accurate noncovalent interaction energies of large molecules on the hundred-atom scale. This enhances the reliability of predictions from quantum mechanical many-electron theories across a wide range of critical applications, including drug design, catalysis, and the innovation of new functional materials, such as those for renewable energy technologies.
△ Less
Submitted 4 July, 2024; v1 submitted 1 July, 2024;
originally announced July 2024.
-
Simulations of stochastic fluid dynamics near a critical point in the phase diagram
Authors:
Chandrodoy Chattopadhyay,
Josh Ott,
Thomas Schaefer,
Vladimir V. Skokov
Abstract:
We present simulations of stochastic fluid dynamics in the vicinity of a critical endpoint belonging to the universality class of the Ising model. This study is motivated by the challenge of modeling the dynamics of critical fluctuations near a conjectured critical endpoint in the phase diagram of Quantum Chromodynamics (QCD). We focus on the interaction of shear modes with a conserved scalar dens…
▽ More
We present simulations of stochastic fluid dynamics in the vicinity of a critical endpoint belonging to the universality class of the Ising model. This study is motivated by the challenge of modeling the dynamics of critical fluctuations near a conjectured critical endpoint in the phase diagram of Quantum Chromodynamics (QCD). We focus on the interaction of shear modes with a conserved scalar density, which is known as model H. We show that the observed dynamical scaling behavior depends on the correlation length and the shear viscosity of the fluid. As the correlation length is increased or the viscosity is decreased we observe a cross-over from the dynamical exponent of critical diffusion, $z\simeq 4$, to the expected scaling exponent of model H, $z\simeq 3$. We use our method to investigate time-dependent correlation function of non-Gaussian moments $M^n(t)$ of the order parameter. We find that the relaxation time depends in non-trivial manner on the power $n$.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Instantons, fluctuations and singularities in the supercritical stochastic nonlinear Schroedinger equation
Authors:
Sumeja Burekovic,
Tobias Schaefer,
Rainer Grauer
Abstract:
Recently, Josserand et al. proposed a stochastic nonlinear Schroedinger model for finite-time singularity-mediated turbulence [Phys. Rev. Fluids 5, 054607 (2020)]. Here, we use instanton calculus to quantify the effect of extreme fluctuations on the statistics of the energy dissipation rate. While the contribution of the instanton alone is insufficient, we obtain excellent agreement with direct si…
▽ More
Recently, Josserand et al. proposed a stochastic nonlinear Schroedinger model for finite-time singularity-mediated turbulence [Phys. Rev. Fluids 5, 054607 (2020)]. Here, we use instanton calculus to quantify the effect of extreme fluctuations on the statistics of the energy dissipation rate. While the contribution of the instanton alone is insufficient, we obtain excellent agreement with direct simulations when including Gaussian fluctuations and the corresponding zero mode. Fluctuations are crucial to obtain the correct scaling when quasi-singular events govern the turbulence statistics.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
Sampling the reciprocal Coulomb potential in finite and anisotropic cells
Authors:
Tobias Schäfer,
William Z. Van Benschoten,
James J. Shepherd,
Andreas Grüneis
Abstract:
We present a robust strategy to numerically sample the Coulomb potential in reciprocal space for periodic Born-von Karman cells of general shape. Our approach tackles two common issues of plane-wave based implementations of Coulomb integrals under periodic boundary conditions, the treatment of the singularity at the Brillouin-zone center, as well as quadrature errors, which can cause severe conver…
▽ More
We present a robust strategy to numerically sample the Coulomb potential in reciprocal space for periodic Born-von Karman cells of general shape. Our approach tackles two common issues of plane-wave based implementations of Coulomb integrals under periodic boundary conditions, the treatment of the singularity at the Brillouin-zone center, as well as quadrature errors, which can cause severe convergence problems in anisotropic cells, necessary for the calculation of low-dimensional systems. We apply our strategy to the Hartree-Fock (HF) and coupled cluster (CC) theory and discuss the consequences of different sampling strategies on the different theories. We show that sampling the Coulomb potential via the widely used probe-charge Ewald method is unsuitable for CC calculations in anisotropic cells. To demonstrate the applicability of our developed approach, we study two representative, low-dimensional use cases: the infinite carbon chain, for which we report the first periodic CCSD(T) potential energy surface, as well as a surface slab of lithium hydride, for which we demonstrate the impact of different sampling strategies for calculating surface energies. We find that our Coulomb sampling strategy serves as a vital solution, addressing the critical need for improved accuracy in plane-wave based CC calculations for low-dimensional systems.
△ Less
Submitted 23 October, 2023; v1 submitted 20 October, 2023;
originally announced October 2023.
-
Cerium Oxides without $U$: The Role of Many-Electron Correlation
Authors:
Tobias Schäfer,
Nathan Daelman,
Núria López
Abstract:
Electron transfer with changing occupation in the 4f subshell poses a considerable challenge for quantitative predictions in quantum chemistry. Using the example of cerium oxide, we identify the main deficiencies of common parameter-dependent one-electron approaches, such as density functional theory (DFT) with a Hubbard correction, or hybrid functionals. As a response, we present the first benchm…
▽ More
Electron transfer with changing occupation in the 4f subshell poses a considerable challenge for quantitative predictions in quantum chemistry. Using the example of cerium oxide, we identify the main deficiencies of common parameter-dependent one-electron approaches, such as density functional theory (DFT) with a Hubbard correction, or hybrid functionals. As a response, we present the first benchmark of ab initio many-electron theory for electron transfer energies and lattice parameters under periodic boundary conditions. We show that the direct random phase approximation clearly outperforms all DFT variations. From this foundation, we, then, systematically improve even further. Periodic second-order Møller-Plesset perturbation theory meanwhile manages to recover standard hybrid functional values. Using these approaches to eliminate parameter bias allows for highly accurate benchmarks of strongly correlated materials, the reliable assessment of various density functionals, and functional fitting via machine-learning.
△ Less
Submitted 6 December, 2021;
originally announced December 2021.
-
Surface science using coupled cluster theory via local Wannier functions and in-RPA-embedding: the case of water on graphitic carbon nitride
Authors:
Tobias Schäfer,
Alejandro Gallo,
Andreas Irmler,
Felix Hummel,
Andreas Grüneis
Abstract:
A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride employing an embedding approach is presented. The embedding approach involves an algorithm to obtain localized Wannier orbitals of various types expanded in a plane-wave basis and intrinsic atomic orbital projectors. The localized occupied orbitals are employed in combination with unoccupie…
▽ More
A first-principles study of the adsorption of a single water molecule on a layer of graphitic carbon nitride employing an embedding approach is presented. The embedding approach involves an algorithm to obtain localized Wannier orbitals of various types expanded in a plane-wave basis and intrinsic atomic orbital projectors. The localized occupied orbitals are employed in combination with unoccupied natural orbitals to perform many-electron perturbation theory calculations of local fragments. The fragments are comprised of a set of localized orbitals close to the adsorbed water molecule. Although the surface model contains more than 100 atoms in the simulation cell, the employed fragments are small enough to allow for calculations using high-level theories up to the coupled cluster ansatz with single, double and perturbative triple particle-hole excitation operators (CCSD(T)). To correct for the missing long-range correlation energy contributions to the adsorption energy, we embed CCSD(T) theory into the direct random phase approximation, yielding rapidly convergent adsorption energies with respect to the fragment size. Convergence of computed binding energies with respect to the virtual orbital basis set is achieved employing a number of recently developed techniques. Moreover, we discuss fragment size convergence for a range of approximate many-electron perturbation theories. The obtained benchmark results are compared to a number of density functional calculations.
△ Less
Submitted 14 October, 2021; v1 submitted 12 October, 2021;
originally announced October 2021.
-
GNOME and LBM Model Evaluation on Ocean Oil Spill Far-Field Impacts to Highly Sensitive Areas
Authors:
Zhanyang Zhang,
Tobias Schaefer,
Michael E. Kress
Abstract:
In case of an ocean oil spill, there are certain areas, e.g. shrimp farms, which are highly sensitive to small amounts of oil pollution while they are geographically far from the spill. We investigate the Lattice Boltzmann Method (LBM) and GNOME, a tool developed and used by NOAA, in terms of far-field impacts to sensitive areas. We present our simulation results of both models in limited scale (a…
▽ More
In case of an ocean oil spill, there are certain areas, e.g. shrimp farms, which are highly sensitive to small amounts of oil pollution while they are geographically far from the spill. We investigate the Lattice Boltzmann Method (LBM) and GNOME, a tool developed and used by NOAA, in terms of far-field impacts to sensitive areas. We present our simulation results of both models in limited scale (a sub area of Gulf of Mexico) under the same oil spill condition using real ocean current data from the Unified Wave Interface-Coupled Model (UWIN-CM). Our study shows that the kinetic theory based LBM model outperforms the stochastic particle based GNOME model in accuracy and computation time due to their fundamental difference in representation of oil pollution advection and diffusion mechanisms. We propose LBM as a viable alternative to the Lagrangian particle calculation component of GNOME in modeling the far-field impacts to highly sensitive areas.
△ Less
Submitted 11 May, 2021;
originally announced May 2021.
-
Effective Hamiltonians for the study of real metals using quantum chemical theories
Authors:
Tina N. Mihm,
Tobias Schäfer,
Sai Kumar Ramadugu,
Andreas Grüneis,
James J. Shepherd
Abstract:
Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are at the heart of computational materials science. In that respect the coupled cluster hierarchy of methods plays a central role in molecular quantum chemistry because of its systematic improvability and computational efficiency. In this hierarchy, coupled cluster singles and…
▽ More
Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are at the heart of computational materials science. In that respect the coupled cluster hierarchy of methods plays a central role in molecular quantum chemistry because of its systematic improvability and computational efficiency. In this hierarchy, coupled cluster singles and doubles (CCSD) is one of the most important steps in moving towards chemical accuracy and, in recent years, its scope has successfully been expanded to the study of insulating surfaces and solids. Here, we show that CCSD theory can also be applied to real metals. In so doing, we overcome the limitation of needing extremely large supercells to capture long range electronic correlation effects. An effective Hamiltonian can be found using the transition structure factor--a map of electronic excitations from the Hartree--Fock wavefunction--which has fewer finite size effects than conventional periodic boundary conditions. This not only paves the way of applying coupled cluster methods to real metals but also reduces the computational cost by two orders of magnitude compared to previous methods. Our applications to phases of lithium and silicon show a resounding success in reaching the thermodynamic limit, taking the first step towards a truly universal quantum chemical treatment of solids.
△ Less
Submitted 10 March, 2021;
originally announced March 2021.
-
Sharp Asymptotic Estimates for Expectations, Probabilities, and Mean First Passage Times in Stochastic Systems with Small Noise
Authors:
Tobias Grafke,
Tobias Schäfer,
Eric Vanden-Eijnden
Abstract:
Freidlin-Wentzell theory of large deviations can be used to compute the likelihood of extreme or rare events in stochastic dynamical systems via the solution of an optimization problem. The approach gives exponential estimates that often need to be refined via calculation of a prefactor. Here it is shown how to perform these computations in practice. Specifically, sharp asymptotic estimates are de…
▽ More
Freidlin-Wentzell theory of large deviations can be used to compute the likelihood of extreme or rare events in stochastic dynamical systems via the solution of an optimization problem. The approach gives exponential estimates that often need to be refined via calculation of a prefactor. Here it is shown how to perform these computations in practice. Specifically, sharp asymptotic estimates are derived for expectations, probabilities, and mean first passage times in a form that is geared towards numerical purposes: they require solving well-posed matrix Riccati equations involving the minimizer of the Freidlin-Wentzell action as input, either forward or backward in time with appropriate initial or final conditions tailored to the estimate at hand. The usefulness of our approach is illustrated on several examples. In particular, invariant measure probabilities and mean first passage times are calculated in models involving stochastic partial differential equations of reaction-advection-diffusion type.
△ Less
Submitted 15 September, 2021; v1 submitted 8 March, 2021;
originally announced March 2021.
-
Quantum limited sound attenuation in a dilute atomic Fermi gas
Authors:
Thomas Schaefer
Abstract:
A new experiment (P. B. Patel et al., Science 370 (2020) 6521, 1222) involving resonantly interacting atoms confined by laser beams sheds light on momentum and energy diffusion in quantum fluids.
A new experiment (P. B. Patel et al., Science 370 (2020) 6521, 1222) involving resonantly interacting atoms confined by laser beams sheds light on momentum and energy diffusion in quantum fluids.
△ Less
Submitted 21 December, 2020;
originally announced December 2020.
-
Local embedding of Coupled Cluster theory into the Random Phase Approximation using plane-waves
Authors:
Tobias Schäfer,
Florian Libisch,
Georg Kresse,
Andreas Grüneis
Abstract:
We present an embedding approach to treat local electron correlation effects in periodic environments. In a single, consistent framework, our plane-wave based scheme embeds a local high-level correlation calculation (here Coupled Cluster Theory, CC), employing localized orbitals, into a low-level correlation calculation (here the direct Random Phase Approximation, RPA). This choice allows for an a…
▽ More
We present an embedding approach to treat local electron correlation effects in periodic environments. In a single, consistent framework, our plane-wave based scheme embeds a local high-level correlation calculation (here Coupled Cluster Theory, CC), employing localized orbitals, into a low-level correlation calculation (here the direct Random Phase Approximation, RPA). This choice allows for an accurate and effcient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersion are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite, for instance, is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons.
△ Less
Submitted 11 December, 2020;
originally announced December 2020.
-
Transasymptotics and hydrodynamization of the Fokker-Planck equation for gluons
Authors:
A. Behtash,
S. Kamata,
M. Martinez,
T. Schaefer,
V. Skokov
Abstract:
We investigate the non-linear transport processes and hydrodynamization of a system of gluons undergoing longitudinal boost-invariant expansion. The dynamics is described within the framework of the Boltzmann equation in the small-angle approximation. The kinetic equations for a suitable set of moments of the one-particle distribution function are derived. By investigating the stability and asympt…
▽ More
We investigate the non-linear transport processes and hydrodynamization of a system of gluons undergoing longitudinal boost-invariant expansion. The dynamics is described within the framework of the Boltzmann equation in the small-angle approximation. The kinetic equations for a suitable set of moments of the one-particle distribution function are derived. By investigating the stability and asymptotic resurgent properties of this dynamical system, we demonstrate, that its solutions exhibit a rather different behavior for large (UV) and small (IR) effective Knudsen numbers. Close to the forward attractor in the IR regime the constitutive relations of each moment can be written as a multiparameter transseries. This resummation scheme allows us to extend the definition of a transport coefficient to the non-equilibrium regime naturally. Each transport coefficient is renormalized by the non-perturbative contributions of the non-hydrodynamic modes. The Knudsen number dependence of the transport coefficient is governed by the corresponding renormalization group flow equation. An interesting feature of the Yang-Mills plasma in this regime is that it exhibits transient non-Newtonian behavior while hydrodynamizing. In the UV regime the solution for the moments can be written as a power-law asymptotic series with a finite radius of convergence. We show that radius of convergence of the UV perturbative expansion grows linearly as a function of the shear viscosity to entropy density ratio. Finally, we compare the universal properties in the pullback and forward attracting regions to other kinetic models including the relaxation time approximation and the effective kinetic Arnold-Moore-Yaffe (AMY) theory.
△ Less
Submitted 15 March, 2021; v1 submitted 16 November, 2020;
originally announced November 2020.
-
Tracking the Footprints of Spin Fluctuations: A MultiMethod, MultiMessenger Study of the Two-Dimensional Hubbard Model
Authors:
Thomas Schäfer,
Nils Wentzell,
Fedor Šimkovic IV,
Yuan-Yao He,
Cornelia Hille,
Marcel Klett,
Christian J. Eckhardt,
Behnam Arzhang,
Viktor Harkov,
François-Marie Le Régent,
Alfred Kirsch,
Yan Wang,
Aaram J. Kim,
Evgeny Kozik,
Evgeny A. Stepanov,
Anna Kauch,
Sabine Andergassen,
Philipp Hansmann,
Daniel Rohe,
Yuri M. Vilk,
James P. F. LeBlanc,
Shiwei Zhang,
A. -M. S. Tremblay,
Michel Ferrero,
Olivier Parcollet
, et al. (1 additional authors not shown)
Abstract:
The Hubbard model represents the fundamental model for interacting quantum systems and electronic correlations. Using the two-dimensional half-filled Hubbard model at weak coupling as a testing ground, we perform a comparative study of a comprehensive set of state of the art quantum many-body methods. Upon cooling into its insulating antiferromagnetic ground-state, the model hosts a rich sequence…
▽ More
The Hubbard model represents the fundamental model for interacting quantum systems and electronic correlations. Using the two-dimensional half-filled Hubbard model at weak coupling as a testing ground, we perform a comparative study of a comprehensive set of state of the art quantum many-body methods. Upon cooling into its insulating antiferromagnetic ground-state, the model hosts a rich sequence of distinct physical regimes with crossovers between a high-temperature incoherent regime, an intermediate temperature metallic regime and a low-temperature insulating regime with a pseudogap created by antiferromagnetic fluctuations. We assess the ability of each method to properly address these physical regimes and crossovers through the computation of several observables probing both quasiparticle properties and magnetic correlations, with two numerically exact methods (diagrammatic and determinantal quantum Monte Carlo) serving as a benchmark. By combining computational results and analytical insights, we elucidate the nature and role of spin fluctuations in each of these regimes. Based on this analysis, we explain how quasiparticles can coexist with increasingly long-range antiferromagnetic correlations, and why dynamical mean-field theory is found to provide a remarkably accurate approximation of local quantities in the metallic regime. We also critically discuss whether imaginary time methods are able to capture the non-Fermi liquid singularities of this fully nested system.
△ Less
Submitted 29 March, 2021; v1 submitted 18 June, 2020;
originally announced June 2020.
-
High-resolution resonance-enhanced multiphoton photoelectron circular dichroism
Authors:
Alexander Kastner,
Greta Koumarianou,
Pavle Glodic,
Peter C. Samartzis,
Nicolas Ladda,
Simon T. Ranecky,
Tom Ring,
Vasudevan Sudheendran,
Constantin Witte,
Hendrike Braun,
Han-Gyeol Lee,
Arne Senftleben,
Robert Berger,
G. Barratt Park,
Tim Schäfer,
Thomas Baumert
Abstract:
Photoelectron circular dichroism (PECD) is a highly sensitive enantiospecific spectroscopy for studying chiral molecules in the gas phase using either single-photon ionization or multiphoton ionization. In the short pulse limit investigated with femtosecond lasers, resonance-enhanced multiphoton ionization (REMPI) is rather instantaneous and typically occurs simultaneously via more than one vibrat…
▽ More
Photoelectron circular dichroism (PECD) is a highly sensitive enantiospecific spectroscopy for studying chiral molecules in the gas phase using either single-photon ionization or multiphoton ionization. In the short pulse limit investigated with femtosecond lasers, resonance-enhanced multiphoton ionization (REMPI) is rather instantaneous and typically occurs simultaneously via more than one vibrational or electronic intermediate state due to limited frequency resolution. In contrast, vibrational resolution in the REMPI spectrum can be achieved using nanosecond lasers. In this work, we follow the high-resolution approach using a tunable narrow-band nanosecond laser to measure REMPI-PECD through distinct vibrational levels in the intermediate 3s and 3p Rydberg states of fenchone. We observe the PECD to be essentially independent of the vibrational level. This behaviour of the chiral sensitivity may pave the way for enantiomer specific molecular identification in multi-component mixtures: one can specifically excite a sharp, vibrationally resolved transition of a distinct molecule to distinguish different chiral species in mixtures.
△ Less
Submitted 31 January, 2020;
originally announced January 2020.
-
Gate-Tunable Graphene Hall Sensors with High Magnetic Field Sensitivity
Authors:
Brian T. Schaefer,
Lei Wang,
Alexander Jarjour,
Kenji Watanabe,
Takashi Taniguchi,
Paul L. McEuen,
Katja C. Nowack
Abstract:
Solid-state magnetic field sensors are important to both modern electronics and fundamental materials science. Many types of these sensors maintain high sensitivity only in a limited range of temperature and background magnetic field, but Hall-effect sensors are in principle able to operate over a broad range of these conditions. Here, we fabricate and characterize micrometer-scale graphene Hall s…
▽ More
Solid-state magnetic field sensors are important to both modern electronics and fundamental materials science. Many types of these sensors maintain high sensitivity only in a limited range of temperature and background magnetic field, but Hall-effect sensors are in principle able to operate over a broad range of these conditions. Here, we fabricate and characterize micrometer-scale graphene Hall sensors demonstrating high magnetic field sensitivity from liquid-helium to room temperature and in background magnetic field up to several Tesla. By tuning the charge carrier density with an electrostatic gate, we optimize the magnetic field sensitivity for different working conditions. From measurements of the Hall coefficient and the Hall voltage noise at 1 kHz, we estimate an optimum magnetic field sensitivity of 80 nT Hz$^{-1/2}$ at 4.2 K, 700 nT Hz$^{-1/2}$ at room temperature, and 3 $μ$T Hz$^{-1/2}$ in 3 T background magnetic field at 4.2 K. Our devices perform competitively with the best existing Hall sensor technologies at room temperature, outperform any Hall sensors reported in the literature at 4.2 K, and demonstrate high sensitivity for the first time in a few Tesla applied magnetic field.
△ Less
Submitted 29 December, 2019;
originally announced December 2019.
-
RPA natural orbitals and their application to post-Hartree-Fock electronic structure methods
Authors:
Benjamin Ramberger,
Zoran Sukurma,
Tobias Schäfer,
Georg Kresse
Abstract:
We present a method to approximate post-Hartree-Fock correlation energies by using approximate natural orbitals obtained by the random phase approximation (RPA). We demonstrate the method by applying it to the helium atom, the hydrogen and fluorine molecule, and to diamond as an example of a periodic system. For these benchmark systems, we show that RPA natural orbitals converge the MP2 correlatio…
▽ More
We present a method to approximate post-Hartree-Fock correlation energies by using approximate natural orbitals obtained by the random phase approximation (RPA). We demonstrate the method by applying it to the helium atom, the hydrogen and fluorine molecule, and to diamond as an example of a periodic system. For these benchmark systems, we show that RPA natural orbitals converge the MP2 correlation energy rapidly. Additionally, we calculated full configuration interaction energies for He and H$_2$, which are in excellent agreement with the literature and experimental values. We conclude that the proposed method may serve as a compromise to reach good approximations to correlation energies at moderate computational cost, and we expect the method to be especially useful for theoretical studies on surface chemistry by providing an efficient basis to correlated wave function based methods.
△ Less
Submitted 18 March, 2021; v1 submitted 16 September, 2019;
originally announced September 2019.
-
Towards a background-free neutrinoless double beta decay experiment based on a fluorescent bicolor sensor
Authors:
Iván Rivilla,
Borja Aparicio,
Juan M. Bueno,
David Casanova,
Claire Tonnelé,
Zoraida Freixa,
Pablo Herrero,
José I. Miranda,
Rosa M. Martínez-Ojeda,
Francesc Monrabal,
Beñat Olave,
Thomas Schäfer,
Pablo Artal,
David Nygren,
Fernando P. Cossío,
Juan J. Gómez-Cadenas
Abstract:
Searching for neutrinoless double beta decays ($β\beta0ν$) is the only practical way to establish if the neutrinos are their own antiparticles. Due to the smallness of neutrino masses, the lifetime of $β\beta0ν$ is expected to be at least ten orders of magnitude smaller than the noise associated with the natural radioactive chains. A positive identification of $β\beta0ν$ decays requires, ultimatel…
▽ More
Searching for neutrinoless double beta decays ($β\beta0ν$) is the only practical way to establish if the neutrinos are their own antiparticles. Due to the smallness of neutrino masses, the lifetime of $β\beta0ν$ is expected to be at least ten orders of magnitude smaller than the noise associated with the natural radioactive chains. A positive identification of $β\beta0ν$ decays requires, ultimately, finding a signal that cannot be mimicked by radioactive backgrounds. This signal could be the observation of the daughter atom in the decay, since no known background processes induce a Z+2 transformation. In particular, the $β\beta0ν$ decay of Xe-136 could be established by detecting the doubly ionised daughter atom, Ba$^{2+}$. Such a detection could be achieved via a sensor made of a monolayer of molecular indicators. The Ba$^{2+}$ would be captured by one of the molecules in the sensor, and the presence of the single chelated indicator would be subsequently revealed by a strong fluorescent response from repeated interrogation with a laser system. Here we describe a fluorescent bicolor indicator that binds strongly to Ba$^{2+}$ and shines very brightly, shifting its emission colour from green to blue when chelated in dry medium, thus allowing the unambiguous identification of single barium atoms in the sensor, and permitting a positive identification of the $β\beta0ν$ decay of Xe-136 in a gas chamber, that could led to a background-free experiment.
△ Less
Submitted 14 September, 2019; v1 submitted 6 September, 2019;
originally announced September 2019.
-
Poly-Sarcosine and Poly(ethylene-glycol) interactions with proteins investigated using molecular dynamics simulations
Authors:
Giovanni Settanni,
Timo Schäfer,
Christian Muhl,
Matthias Barz,
Friederike Schmid
Abstract:
Nanoparticles coated with hydrophilic polymers often show a reduction in unspecific interactions with the biological environment, which improves their biocompatibility. The molecular determinants of this reduction are not very well understood yet, and their knowledge may help improving nanoparticle design. Here we address, using molecular dynamics simulations, the interactions of human serum album…
▽ More
Nanoparticles coated with hydrophilic polymers often show a reduction in unspecific interactions with the biological environment, which improves their biocompatibility. The molecular determinants of this reduction are not very well understood yet, and their knowledge may help improving nanoparticle design. Here we address, using molecular dynamics simulations, the interactions of human serum albumin, the most abundant serum protein, with two promising hydrophilic polymers used for the coating of therapeutic nanoparticles, poly(ethylene-glycol) and poly-sarcosine. By simulating the protein immersed in a polymer-water mixture, we show that the two polymers have a very similar affinity for the protein surface, both in terms of the amount of polymer adsorbed and also in terms of the type of amino acids mainly involved in the interactions. We further analyze the kinetics of adsorption and how it affects the polymer conformations. Minor differences between the polymers are observed in the thickness of the adsorption layer, that are related to the different degree of flexibility of the two molecules. In comparison poly-alanine, an isomer of poly-sarcosine known to self-aggregate and induce protein aggregation, shows a significantly larger affinity for the protein surface than PEG and PSar, which we show to be related not to a different patterns of interactions with the protein surface, but to the different way the polymer interacts with water.
△ Less
Submitted 29 October, 2018;
originally announced October 2018.
-
Ab initio phase diagram of PbSe crystals calculated with the Random Phase Approximation
Authors:
Tobias Schäfer,
Zhaochuan Fan,
Michael Grünwald,
Georg Kresse
Abstract:
Understanding the phase behavior of semiconductor materials is important for applications in solid state physics and nanoscience. Accurate experimental data is often difficult to obtain due to strong kinetic effects. In this work, we calculate the temperature-pressure phase diagram of lead selenide (PbSe) using the random phase approximation (RPA), an accurate wavefunction based many-body techniqu…
▽ More
Understanding the phase behavior of semiconductor materials is important for applications in solid state physics and nanoscience. Accurate experimental data is often difficult to obtain due to strong kinetic effects. In this work, we calculate the temperature-pressure phase diagram of lead selenide (PbSe) using the random phase approximation (RPA), an accurate wavefunction based many-body technique. We consider three crystalline phases, the low pressure B1 phase (NaCl-type), the intermediate B33 phase (CrB-type), and the high pressure B2 phase (CsCl-type). The electronic contributions to the free energy (at T=0K) are calculated in the Born-Oppenheimer approximation using the RPA, whereas phononic contributions are computed in the quasi-harmonic approximation using DFT and the PBEsol functional. At room temperature, we find transition pressures of 4.6 +/- 0.3 GPa for the B1-B33 transition and 18.7 +/- 0.3 GPa for the B33-B2 transition, in good agreement with experiments. In contrast to the interpretation of recent experiments, we observe a negative Clapeyron slope for both transitions. Gibbs free energy differences between competing structures have small gradients close to coexistence, consistent with pronounced hysteresis observed in experiments. The phase diagram presented in this work can serve as a reference for future studies of PbSe and should prove useful in the development of accurate and efficient force fields.
△ Less
Submitted 30 July, 2018;
originally announced July 2018.
-
Laplace transformed MP2 for three dimensional periodic materials using stochastic orbitals in the plane wave basis and correlated sampling
Authors:
Tobias Schäfer,
Benjamin Ramberger,
Georg Kresse
Abstract:
We present an implementation and analysis of a stochastic high performance algorithm to calculate the correlation energy of three dimensional periodic systems in second-order Møller-Plesset perturbation theory (MP2). In particular we measure the scaling behavior of the sample variance and probe whether this stochastic approach is competitive if accuracies well below 1 meV per valence orbital are r…
▽ More
We present an implementation and analysis of a stochastic high performance algorithm to calculate the correlation energy of three dimensional periodic systems in second-order Møller-Plesset perturbation theory (MP2). In particular we measure the scaling behavior of the sample variance and probe whether this stochastic approach is competitive if accuracies well below 1 meV per valence orbital are required, as it is necessary for calculations of adsorption, binding, or surface energies. The algorithm is based on the Laplace transformed MP2 (LTMP2) formulation in the plane wave basis. The time-dependent Hartree-Fock orbitals, appearing in the LTMP2 formulation, are stochastically rotated in the occupied and unoccupied Hilbert space. This avoids a full summation over all combinations of occupied and unoccupied orbitals, as inspired by the work of D. Neuhauser, E. Rabani, and R. Baer in J. Chem. Theory Comput. 9, 24 (2013). Additionally, correlated sampling is introduced, accelerating the statistical convergence significantly.
△ Less
Submitted 9 February, 2018; v1 submitted 15 November, 2017;
originally announced November 2017.
-
Intermediate state dependence of the photoelectron circular dichroism of fenchone observed via femtosecond resonance-enhanced multi-photon ionization
Authors:
Alexander Kastner,
Tom Ring,
Bastian C. Krüger,
G. Barratt Park,
Tim Schäfer,
Arne Senftleben,
Thomas Baumert
Abstract:
The intermediate state dependence of photoelectron circular dichroism (PECD) in resonance-enhanced multi-photon ionization of fenchone in the gas phase is experimentally studied. By scanning the excitation wavelength from 359 to 431 nm we simultaneously excite up to three electronically distinct resonances. In the PECD experiment performed with a broadband femtosecond laser their respective contri…
▽ More
The intermediate state dependence of photoelectron circular dichroism (PECD) in resonance-enhanced multi-photon ionization of fenchone in the gas phase is experimentally studied. By scanning the excitation wavelength from 359 to 431 nm we simultaneously excite up to three electronically distinct resonances. In the PECD experiment performed with a broadband femtosecond laser their respective contributions to the photoelectron spectrum can be resolved. High-resolution spectroscopy allows us to identify two of the resonances as belonging to the B- and C-bands, which involve excitation to states with 3s and 3p Rydberg character, respectively. We observe a sign change in the PECD signal depending on which electronic state is used as an intermediate. Additionally, scanning the laser wavelength reveals a decrease of PECD magnitude with increasing photoelectron energy for the 3s state.
△ Less
Submitted 20 February, 2017;
originally announced February 2017.
-
Quartic scaling MP2 for solids: A highly parallelized algorithm in the plane wave basis
Authors:
Tobias Schäfer,
Benjamin Ramberger,
Georg Kresse
Abstract:
We present a low-complexity algorithm to calculate the correlation energy of periodic systems in second-order Møller-Plesset perturbation theory (MP2). In contrast to previous approximation-free MP2 codes, our implementation possesses a quartic scaling, $\mathcal O (N^4)$, with respect to the system size $N$ and offers an almost ideal parallelization efficiency. The general issue that the correlat…
▽ More
We present a low-complexity algorithm to calculate the correlation energy of periodic systems in second-order Møller-Plesset perturbation theory (MP2). In contrast to previous approximation-free MP2 codes, our implementation possesses a quartic scaling, $\mathcal O (N^4)$, with respect to the system size $N$ and offers an almost ideal parallelization efficiency. The general issue that the correlation energy converges slowly with the number of basis functions is eased by an internal basis set extrapolation. The key concept to reduce the scaling is to eliminate all summations over virtual orbitals which can be elegantly achieved in the Laplace transformed MP2 (LTMP2) formulation using plane wave basis sets and Fast Fourier transforms. Analogously, this approach could allow to calculate second order screened exchange (SOSEX) as well as particle-hole ladder diagrams with a similar low complexity. Hence, the presented method can be considered as a step towards systematically improved correlation energies.
△ Less
Submitted 15 March, 2017; v1 submitted 21 November, 2016;
originally announced November 2016.
-
Analytic Interatomic Forces in the Random Phase Approximation
Authors:
Benjamin Ramberger,
Tobias Schäfer,
Georg Kresse
Abstract:
We discuss that in the random phase approximation (RPA) the first derivative of the energy with respect to the Green's function is the self-energy in the GW approximation. This relationship allows us to derive compact equations for the RPA interatomic forces. We also show that position dependent overlap operators are elegantly incorporated in the present framework. The RPA force equations have bee…
▽ More
We discuss that in the random phase approximation (RPA) the first derivative of the energy with respect to the Green's function is the self-energy in the GW approximation. This relationship allows us to derive compact equations for the RPA interatomic forces. We also show that position dependent overlap operators are elegantly incorporated in the present framework. The RPA force equations have been implemented in the projector augmented wave formalism, and we present illustrative applications, including ab initio molecular dynamics simulations, the calculation of phonon dispersion relations for diamond and graphite, as well as structural relaxations for water on boron nitride. The present derivation establishes a concise framework for forces within perturbative approaches and is also applicable to more involved approximations for the correlation energy.
△ Less
Submitted 16 March, 2017; v1 submitted 2 November, 2016;
originally announced November 2016.
-
Longitudinal and transverse structure functions in high Reynolds-number magneto-hydrodynamic turbulence
Authors:
J. Friedrich,
H. Homann,
T. Schäfer,
R. Grauer
Abstract:
We investigate the scaling behavior of longitudinal and transverse structure functions in homogeneous and isotropic magneto-hydrodynamic (MHD) turbulence by means of an exact hierarchy of structure function equations as well as by direct numerical simulations of two- and three-dimensional MHD turbulence. In particular, rescaling relations between longitudinal and transverse structure functions are…
▽ More
We investigate the scaling behavior of longitudinal and transverse structure functions in homogeneous and isotropic magneto-hydrodynamic (MHD) turbulence by means of an exact hierarchy of structure function equations as well as by direct numerical simulations of two- and three-dimensional MHD turbulence. In particular, rescaling relations between longitudinal and transverse structure functions are derived and utilized in order to compare different scaling behavior in the inertial range. It is found that there are no substantial differences between longitudinal and transverse structure functions in MHD turbulence. This finding stands in contrast to the case of hydrodynamic turbulence which shows persistent differences even at high Reynolds numbers. We propose a physical picture that is based on an effective reduction of pressure contributions due to local regions of same magnitude and alignment of velocity and magnetic field fluctuations. Finally, our findings underline the importance of the pressure term for the actually observed scaling differences in hydrodynamic turbulence.
△ Less
Submitted 2 May, 2017; v1 submitted 19 September, 2016;
originally announced September 2016.
-
A generalized Theory of Diffusion based on Kinetic Theory
Authors:
Thomas Schaefer
Abstract:
We propose to use spin hydrodynamics, a two-fluid model of spin propagation, as a generalization of the diffusion equation. We show that in the dense limit spin hydrodynamics reduces to Fick's law and the diffusion equation. In the opposite limit spin hydrodynamics is equivalent to a collisionless Boltzmann treatment of spin propagation. Spin hydrodynamics avoids unphysical effects that arise when…
▽ More
We propose to use spin hydrodynamics, a two-fluid model of spin propagation, as a generalization of the diffusion equation. We show that in the dense limit spin hydrodynamics reduces to Fick's law and the diffusion equation. In the opposite limit spin hydrodynamics is equivalent to a collisionless Boltzmann treatment of spin propagation. Spin hydrodynamics avoids unphysical effects that arise when the diffusion equation is used to describe to a strongly interacting gas with a dilute corona. We apply spin hydrodynamics to the problem of spin diffusion in a trapped atomic gas. We find that the observed spin relaxation rate in the high temperature limit [Sommer et al., Nature 472, 201 (2011)] is consistent with the diffusion constant predicted by kinetic theory.
△ Less
Submitted 7 October, 2016; v1 submitted 17 August, 2016;
originally announced August 2016.
-
The instanton method and its numerical implementation in fluid mechanics
Authors:
Tobias Grafke,
Rainer Grauer,
Tobias Schäfer
Abstract:
A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-…
▽ More
A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-Wentzell action and known to be able to characterize rare events in such systems. While there is an impressive body of work concerning their analytical description, this review focuses on the question on how to compute these minimizers numerically. In a short introduction we present the relevant mathematical and physical background before we discuss the stochastic Burgers equation in detail. We present algorithms to compute instantons numerically by an efficient solution of the corresponding Euler-Lagrange equations. A second focus is the discussion of a recently developed numerical filtering technique that allows to extract instantons from direct numerical simulations. In the following we present modifications of the algorithms to make them efficient when applied to two- or three-dimensional fluid dynamical problems. We illustrate these ideas using the two-dimensional Burgers equation and the three-dimensional Navier-Stokes equations.
△ Less
Submitted 20 July, 2015; v1 submitted 29 June, 2015;
originally announced June 2015.
-
Relevance of instantons in Burgers turbulence
Authors:
Tobias Grafke,
Rainer Grauer,
Tobias Schäfer,
Eric Vanden-Eijnden
Abstract:
Instanton calculations are performed in the context of stationary Burgers turbulence to estimate the tails of the probability density function (PDF) of velocity gradients. These results are then compared to those obtained from massive direct numerical simulations (DNS) of the randomly forced Burgers equation. The instanton predictions are shown to agree with the DNS in a wide range of regimes, inc…
▽ More
Instanton calculations are performed in the context of stationary Burgers turbulence to estimate the tails of the probability density function (PDF) of velocity gradients. These results are then compared to those obtained from massive direct numerical simulations (DNS) of the randomly forced Burgers equation. The instanton predictions are shown to agree with the DNS in a wide range of regimes, including those that are far from the limiting cases previously considered in the literature. These results settle the controversy of the relevance of the instanton approach for the prediction of the velocity gradient PDF tail exponents. They also demonstrate the usefulness of the instanton formalism in Burgers turbulence, and suggest that this approach may be applicable in other contexts, such as 2D and 3D turbulence in compressible and incompressible flows.
△ Less
Submitted 30 November, 2014;
originally announced December 2014.
-
Averaging and spectral properties for the 2D advection-diffusion equation in the semi-classical limit for vanishing diffusivity
Authors:
Eugene Dedits,
Andrew C. Poje,
Tobias Schaefer,
Jesenko Vukadinovic
Abstract:
We consider the two-dimensional advection-diffusion equation on a bounded domain subject to either Dirichlet or von Neumann boundary conditions and study both time-independent and time-periodic cases involving Liouville integrable Hamiltonians that satisfy conditions conducive to applying the averaging principle. Transformation to action-angle coordinates permits averaging in time and angle, leadi…
▽ More
We consider the two-dimensional advection-diffusion equation on a bounded domain subject to either Dirichlet or von Neumann boundary conditions and study both time-independent and time-periodic cases involving Liouville integrable Hamiltonians that satisfy conditions conducive to applying the averaging principle. Transformation to action-angle coordinates permits averaging in time and angle, leading to an underlying eigenvalue equation that allows for separation of the angle and action coordinates. The result is a one-dimensional second-order equation involving an anti-symmetric imaginary potential. For radial flows on a disk or an annulus, we rigorously apply existing complex-plane WKBJ methods to study the spectral properties in the semi-classical limit for vanishing diffusivity. In this limit, the spectrum is found to be a complicated set consisting of lines related to Stokes graphs. Eigenvalues in the neighborhood of these graphs exhibit nonlinear scaling with respect to diffusivity leading to convection-enhanced rates of dissipation (relaxation, mixing) for initial data which are mean-free in the angle coordinate. These branches coexist with a diffusive branch of eigenvalues that scale linearly with diffusivity and contain the principal eigenvalue (no dissipation enhancement).
△ Less
Submitted 27 September, 2013;
originally announced September 2013.
-
Arclength parametrized Hamilton's equations for the calculation of instantons
Authors:
Tobias Grafke,
Rainer Grauer,
Tobias Schäfer,
Eric Vanden-Eijnden
Abstract:
A method is presented to compute minimizers (instantons) of action functionals using arclength parametrization of Hamilton's equations. This method can be interpreted as a local variant of the geometric minimum action method (gMAM) introduced to compute minimizers of the Freidlin-Wentzell action functional that arises in the context of large deviation theory for stochastic differential equations.…
▽ More
A method is presented to compute minimizers (instantons) of action functionals using arclength parametrization of Hamilton's equations. This method can be interpreted as a local variant of the geometric minimum action method (gMAM) introduced to compute minimizers of the Freidlin-Wentzell action functional that arises in the context of large deviation theory for stochastic differential equations. The method is particularly well-suited to calculate expectations dominated by noise-induced excursions from deterministically stable fixpoints. Its simplicity and computational efficiency are illustrated here using several examples: a finite-dimensional stochastic dynamical system (an Ornstein-Uhlenbeck model) and two models based on stochastic partial differential equations: the $φ^4$-model and the stochastically driven Burgers equation.
△ Less
Submitted 20 September, 2013;
originally announced September 2013.
-
Ellipsometry with polarisation analysis at cryogenic temperatures inside a vacuum chamber
Authors:
S. Bauer,
B. Grees,
D. Spitzer,
M. Beck,
R. Bottesch,
H. -W. Ortjohann,
B. Ostrick,
T. Schäfer,
H. H. Telle,
A. Wegmann,
M. Zbořil,
C. Weinheimer
Abstract:
In this paper we describe a new variant of null ellipsometry to determine thicknesses and optical properties of thin films on a substrate at cryogenic temperatures. In the PCSA arrangement of ellipsometry the polarizer and the compensator are placed before the substrate and the analyzer after it. Usually, the polarizer and the analyzer are rotated to find the intensity minimum searched for in null…
▽ More
In this paper we describe a new variant of null ellipsometry to determine thicknesses and optical properties of thin films on a substrate at cryogenic temperatures. In the PCSA arrangement of ellipsometry the polarizer and the compensator are placed before the substrate and the analyzer after it. Usually, the polarizer and the analyzer are rotated to find the intensity minimum searched for in null ellipsometry. In our variant we rotate the polarizer and the compensator instead, both being placed in the incoming beam before the substrate. Therefore the polarization analysis of the reflected beam can be realized by an analyzer at fixed orientation. We developed this method for investigations of thin cryogenic films inside a vacuum chamber, where the analyzer and detector had to be placed inside the cold shield at a temperature of T approx. 90K close to the substrate. All other optical components were installed at the incoming beam line outside the vacuum chamber, including all components which need to be rotated during the measurements. Our null ellipsometry variant has been tested with condensed krypton films on a highly oriented pyrolytic graphite substrate (HOPG) at a temperature of T approx. 25K. We show that it is possible to determine the indices of refraction of condensed krypton and of the HOPG substrate as well as thickness of krypton films with reasonable accuracy.
△ Less
Submitted 22 July, 2013;
originally announced July 2013.
-
Automated Image Analysis of Hodgkin lymphoma
Authors:
Alexander Schmitz,
Tim Schäfer,
Hendrik Schäfer,
Claudia Döring,
Jörg Ackermann,
Norbert Dichter,
Sylvia Hartmann,
Martin-Leo Hansmann,
Ina Koch
Abstract:
Hodgkin lymphoma is an unusual type of lymphoma, arising from malignant B-cells. Morphological and immunohistochemical features of malignant cells and their distribution differ from other cancer types. Based on systematic tissue image analysis, computer-aided exploration can provide new insights into Hodgkin lymphoma pathology. In this paper, we report results from an image analysis of CD30 immuno…
▽ More
Hodgkin lymphoma is an unusual type of lymphoma, arising from malignant B-cells. Morphological and immunohistochemical features of malignant cells and their distribution differ from other cancer types. Based on systematic tissue image analysis, computer-aided exploration can provide new insights into Hodgkin lymphoma pathology. In this paper, we report results from an image analysis of CD30 immunostained Hodgkin lymphoma tissue section images. To the best of our knowledge, this is the first systematic application of image analysis to a set of tissue sections of Hodgkin lymphoma. We have implemented an automatic procedure to handle and explore image data in Aperio's SVS format. We use pre-processing approaches on a down-scaled image to separate the image objects from the background. Then, we apply a supervised classification method to assign pixels to predefined classes. Our pre-processing method is able to separate the tissue content of images from the image background. We analyzed three immunohistologically defined groups, non-lymphoma and the two most common forms of Hodgkin lymphoma, nodular sclerosis and mixed cellularity type. We found that nodular sclerosis and non-lymphoma images exhibit different amounts of CD30 stain, whereas mixed cellularity type exhibits a large variance and overlaps with the other groups. The results can be seen as a first step to computationally identify tumor regions in the images. This allows us to focus on these regions when performing computationally expensive tasks like object detection in the high-resolution image.
△ Less
Submitted 14 September, 2012;
originally announced September 2012.
-
Hydrodynamic fluctuations and the minimum shear viscosity of the dilute Fermi gas at unitarity
Authors:
Clifford Chafin,
Thomas Schaefer
Abstract:
We study hydrodynamic fluctuations in a non-relativistic fluid. We show that in three dimensions fluctuations lead to a minimum in the shear viscosity to entropy density ratio $η/s$ as a function of the temperature. The minimum provides a bound on $η/s$ which is independent of the conjectured bound in string theory, $η/s \geq \hbar/(4πk_B)$, where $s$ is the entropy density. For the dilute Fermi g…
▽ More
We study hydrodynamic fluctuations in a non-relativistic fluid. We show that in three dimensions fluctuations lead to a minimum in the shear viscosity to entropy density ratio $η/s$ as a function of the temperature. The minimum provides a bound on $η/s$ which is independent of the conjectured bound in string theory, $η/s \geq \hbar/(4πk_B)$, where $s$ is the entropy density. For the dilute Fermi gas at unitarity we find $η/s\gsim 0.2\hbar$. This bound is not universal -- it depends on thermodynamic properties of the unitary Fermi gas, and on empirical information about the range of validity of hydrodynamics. We also find that the viscous relaxation time of a hydrodynamic mode with frequency $ω$ diverges as $1/\sqrtω$, and that the shear viscosity in two dimensions diverges as $\log(1/ ω)$.
△ Less
Submitted 18 February, 2013; v1 submitted 5 September, 2012;
originally announced September 2012.
-
Instanton filtering for the stochastic Burgers equation
Authors:
Tobias Grafke,
Rainer Grauer,
Tobias Schäfer
Abstract:
We address the question whether one can identify instantons in direct numerical simulations of the stochastically driven Burgers equation. For this purpose, we first solve the instanton equations using the Chernykh-Stepanov method [Phys. Rev. E 64, 026306 (2001)]. These results are then compared to direct numerical simulations by introducing a filtering technique to extract prescribed rare events…
▽ More
We address the question whether one can identify instantons in direct numerical simulations of the stochastically driven Burgers equation. For this purpose, we first solve the instanton equations using the Chernykh-Stepanov method [Phys. Rev. E 64, 026306 (2001)]. These results are then compared to direct numerical simulations by introducing a filtering technique to extract prescribed rare events from massive data sets of realizations. Using this approach we can extract the entire time history of the instanton evolution which allows us to identify the different phases predicted by the direct method of Chernykh and Stepanov with remarkable agreement.
△ Less
Submitted 5 September, 2012;
originally announced September 2012.
-
Conformal symmetry and non-relativistic second order fluid dynamics
Authors:
Jingyi Chao,
Thomas Schaefer
Abstract:
We study the constraints imposed by conformal symmetry on the equations of fluid dynamics at second order in gradients of the hydrodynamic variables. At zeroth order conformal symmetry implies a constraint on the equation of state, E=2/3 P, where E is the energy density and P is the pressure. At first order, conformal symmetry implies that the bulk viscosity must vanish. We show that at second ord…
▽ More
We study the constraints imposed by conformal symmetry on the equations of fluid dynamics at second order in gradients of the hydrodynamic variables. At zeroth order conformal symmetry implies a constraint on the equation of state, E=2/3 P, where E is the energy density and P is the pressure. At first order, conformal symmetry implies that the bulk viscosity must vanish. We show that at second order conformal invariance requires that two-derivative terms in the stress tensor must be traceless, and that it determines the relaxation of dissipative stresses to the Navier-Stokes form. We verify these results by solving the Boltzmann equation at second order in the gradient expansion. We find that only a subset of the terms allowed by conformal symmetry appear.
△ Less
Submitted 24 August, 2011;
originally announced August 2011.
-
Symmetrization of advection-diffusion operators
Authors:
Eugene Dedits,
Andrew C. Poje,
Tobias Schaefer,
Jesenko Vukadinovic
Abstract:
We present a new method to transform an expanded class of non-selfadjoint advection-diffusion operators into self-adjoint operators. The transform is based on a combination of a point transform and Lie transform in conjunction with an asymptotic expansion in terms of the diffusivity. We illustrate the method in the context of simple shear flow where the expansion is exact and all transformation st…
▽ More
We present a new method to transform an expanded class of non-selfadjoint advection-diffusion operators into self-adjoint operators. The transform is based on a combination of a point transform and Lie transform in conjunction with an asymptotic expansion in terms of the diffusivity. We illustrate the method in the context of simple shear flow where the expansion is exact and all transformation steps can be performed explicitly.
△ Less
Submitted 7 December, 2012; v1 submitted 19 February, 2010;
originally announced February 2010.
-
In search of the perfect fluid
Authors:
Thomas Schaefer
Abstract:
Shear viscosity measures the amount of internal friction in a simple fluid. In kinetic theory shear viscosity is related to momentum transport by quasi-particles, and the uncertainty relation implies that the ratio of shear viscosity eta to entropy density s is bounded by a constant multiplied by hbar/k_B, where hbar is Planck's constant and k_B is Boltzmann's constant. A specific bound has been…
▽ More
Shear viscosity measures the amount of internal friction in a simple fluid. In kinetic theory shear viscosity is related to momentum transport by quasi-particles, and the uncertainty relation implies that the ratio of shear viscosity eta to entropy density s is bounded by a constant multiplied by hbar/k_B, where hbar is Planck's constant and k_B is Boltzmann's constant. A specific bound has been proposed on the basis of string theory. In a large class of theories that can be studied using string theory methods the constant is 1/(4π). Experiments at RHIC indicate that eta/s of the quark gluon plasma is close to this prediction. We will refer to a fluid that saturates the string theory bound as a perfect fluid. In this contribution we summarizes the theoretical and experimental information on the fluidity of the main classes of strongly interacting quantum fluids.
△ Less
Submitted 29 June, 2009;
originally announced June 2009.
-
Averaged dynamics of time-periodic advection diffusion equations in the limit of small diffusivity
Authors:
Tobias Schaefer,
Andrew C. Poje,
Jesenko Vukadinovic
Abstract:
We study the effect of advection and small diffusion on passive tracers. The advecting velocity field is assumed to have mean zero and to possess time-periodic stream lines. Using a canonical transform to action-angle variables followed by a Lie-transform, we derive an averaged equation describing the effective motion of the tracers. An estimate for the time validity of the first-order approxima…
▽ More
We study the effect of advection and small diffusion on passive tracers. The advecting velocity field is assumed to have mean zero and to possess time-periodic stream lines. Using a canonical transform to action-angle variables followed by a Lie-transform, we derive an averaged equation describing the effective motion of the tracers. An estimate for the time validity of the first-order approximation is established. For particular cases of a regularized vortical flow we present explicit formulas for the coefficients of the averaged equation both at first and at second order. Numerical simulations indicate that the validity of the above first-order estimate extends to the second order.
△ Less
Submitted 10 June, 2008;
originally announced June 2008.