Source author record

Felix Hummel

Felix Hummel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Geometric blow-up of a dynamic Turing instability in the Swift-Hohenberg equation

We present a rigorous analysis of the slow passage through a Turing bifurcation in the Swift-Hohenberg equation using a novel approach based on geometric blow-up. We show that the formally derived multiple scales ansatz which is known from classical modulation theory can be adapted for use in the fast-slow setting, by reformulating it as a blow-up transformation. This leads to dynamically simpler modulation equations posed in the blown-up space, via a formal procedure which directly extends the established approach to the time-dependent setting. The modulation equations take the form of non-autonomous Ginzburg-Landau equations, which can be analysed within the blow-up. The asymptotics of solutions in weighted Sobelev spaces are given in two different cases: (i) A symmetric case featuring a delayed loss of stability, and (ii) A second case in which the symmetry is broken by a source term. In order to characterise the dynamics of the Swift-Hohenberg equation itself we derive rigorous estimates on the error of the dynamic modulation approximation. These estimates are obtained by bounding weak solutions to an evolution equation for the error which is also posed in the blown-up space. Using the error estimates obtained, we are able to infer the asymptotics of a large class of solutions to the dynamic Swift-Hohenberg equation. We provide rigorous asymptotics for solutions in both cases (i) and (ii). We also prove the existence of the delayed loss of stability in the symmetric case (i), and provide a lower bound for the delay time.

preprint2022arXiv

Interface to high-performance periodic coupled-cluster theory calculations with atom-centered, localized basis functions

Coupled cluster (CC) theory is often considered the gold standard of quantum-chemistry. For solids, however, the available software is scarce. We present CC-aims, which can interface ab initio codes with localized atomic orbitals and the CC for solids (CC4S) code by the group of A. Grüneis. CC4S features a continuously growing selection of wave function-based methods including perturbation and CC theory. The CC-aims interface was developed for the FHI-aims code (https://fhi-aims.org) but is implemented such that other codes may use it as a starting point for corresponding interfaces. As CC4S offers treatment of both molecular and periodic systems, the CC-aims interface is a valuable tool, where DFT is either too inaccurate or too unreliable, in theoretical chemistry and materials science alike.

preprint2021arXiv

Surface science using coupled cluster theory via local Wannier functions and in-RPA-embedding: the case of water on graphitic carbon nitride

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.

preprint2020arXiv

A periodic equation-of-motion coupled-cluster implementation applied to $F$-centers in alkaline earth oxides

We present an implementation of equation of motion coupled-cluster singles and doubles (EOM-CCSD) theory using periodic boundary conditions and a plane wave basis set. Our implementation of EOM-CCSD theory is applied to study $F$-centers in alkaline earth oxides employing a periodic supercell approach. The convergence of calculated electronic excitation energies for neutral color centers in MgO, CaO and SrO crystals with respect to orbital basis set and system size is explored. We discuss extrapolation techniques that approximate excitation energies in the complete basis set limit and reduce finite size errors. Our findings demonstrate that EOM-CCSD theory can predict optical absorption energies of $F$-centers in good agreement with experiment. Furthermore, we discuss calculated emission energies corresponding to the decay from triplet to singlet states, responsible for the photoluminescence properties. Our findings are compared to experimental and theoretical results available in literature.

preprint2020arXiv

Applying the Coupled-Cluster Ansatz to Solids and Surfaces in the Thermodynamic Limit

Modern electronic structure theories can predict and simulate a wealth of phenomena in surface science and solid-state physics. In order to allow for a direct comparison with experiment, such ab initio predictions have to be made in the thermodynamic limit, substantially increasing the computational cost of many-electron wave-function theories. Here, we present a method that achieves thermodynamic limit results for solids and surfaces using the "gold standard" coupled cluster ansatz of quantum chemistry with unprecedented efficiency. We study the energy difference between carbon diamond and graphite crystals, adsorption energies of water on h-BN, as well as the cohesive energy of the Ne solid, demonstrating the increased efficiency and accuracy of coupled cluster theory for solids and surfaces.

preprint2020arXiv

Sample Paths of White Noise in Spaces with Dominating Mixed Smoothness

The sample paths of white noise are proved to be elements of certain Besov spaces with dominating mixed smoothness. Unlike in isotropic spaces, here the regularity does not get worse with increasing space dimension. Consequently, white noise is actually much smoother than the known sharp regularity results in isotropic spaces suggest. An application of our techniques yields new results for the regularity of solutions of Poisson and heat equation on the half space with boundary noise. The main novelty is the flexible treatment of the interplay between the singularity at the boundary and the smoothness in tangential, normal and time direction.

preprint2020arXiv

Slow Manifolds for Infinite-Dimensional Evolution Equations

We extend classical finite-dimensional Fenichel theory in two directions to infinite dimensions. Under comparably weak assumptions we show that the solution of an infinite-dimensional fast-slow system is approximated well by the corresponding slow flow. After that we construct a two-parameter family of slow manifolds $S_{ε,ζ}$ under more restrictive assumptions on the linear part of the slow equation. The second parameter $ζ$ does not appear in the finite-dimensional setting and describes a certain splitting of the slow variable space in a fast decaying part and its complement. The finite-dimensional setting is contained as a special case in which $S_{ε,ζ}$ does not depend on $ζ$. Finally, we apply our new techniques to three examples of fast-slow systems of partial differential equations.

preprint2018arXiv

Dispersive mixed-order systems in $L^p$-Sobolev spaces and application to the thermoelastic plate equation

We study dispersive mixed-order systems of pseudodifferential operators in the setting of $L^p$-Sobolev spaces. Under the weak condition of quasi-hyperbolicity, these operators generate a semigroup in the space of tempered distributions. However, if the basic space is a tuple of $L^p$-Sobolev spaces, a strongly continuous semigroup is in many cases only generated if $p=2$ or $n=1$. The results are applied to the linear thermoelastic plate equation inertial term and with Fourier's or Maxwell-Cattaneo's law of heat conduction.

preprint2015arXiv

Adjacent Pairs Exchange correction to the Random Phase Approximation

The Random Phase Approximation (RPA) is a widely employed post Hartree-Fock or DFT method, capable of capturing van der Waal interactions and other dynamic correlation effects at relatively low costs of $\mathcal O(N^3)$ in time and $\mathcal O(N^2)$ in memory, if calculated from imaginary time propagators. However, since it neglects anti-symmetrization RPA is biased, overestimating the correlation energy and bond lengths in general. The Second Order Screened Exchange offers amelioration by anti-symmetrizing one Coulomb interaction but it comes at considerable costs of $\mathcal O(N^5)$ in time and $\mathcal O(N^4)$ in memory, since it has to be calculated from the direct ring Coupled Cluster Doubles (drCCD) amplitudes. We propose a diagrammatic method, exchanging adjacent pairs in the RPA diagrams - hence its name - offering similar accuracy but with memory requirement of only $\mathcal O(N^2)$. The correction is calculated from imaginary time propagators similar to efficient RPA implementations. It can be calculated in $\mathcal O(N^5)$ steps in momentum space or $\mathcal O(N^4)$ steps in real space. It improves on SOSEX in the uniform electron gas for low densities, where correlation effects are stronger, and first calculations of lattice constants in solids are in excellent agreement with experiment.

preprint2015arXiv

Hidden scale invariance of metals

Density functional theory (DFT) calculations of 58 liquid elements at their triple point show that most metals exhibit near proportionality between thermal fluctuations between virial and potential-energy in the isochoric ensemble. This demonstrates a general "hidden" scale invariance of metals making the dense part of the thermodynamic phase diagram effectively one dimensional with respect to structure and dynamics. DFT computed density scaling exponents, related to the Gr{ü}neisen parameter, are in good agreement with experimental values for 16 elements where reliable data were available. Hidden scale invariance is demonstrated in detail for magnesium by showing invariance of structure and dynamics. Computed melting curves of period three metals follow curves with invariance (isomorphs). The experimental structure factor of magnesium is predicted by assuming scale invariant inverse power-law (IPL) pair interactions. However, crystal packings of several transition metals (V, Cr, Mn, Fe, Nb, Mo, Ta, W and Hg), most post-transition metals (Ga, In, Sn, and Tl) and the metalloids Si and Ge cannot be explained by the IPL assumption. Thus, hidden scale invariance can be present even when the IPL-approximation is inadequate. The virial-energy correlation coefficient of iron and phosphorous is shown to increase at elevated pressures. Finally, we discuss how scale invariance explains the Gr{ü}neisen equation of state and a number of well-known empirical melting and freezing rules.

preprint2014arXiv

Computing the crystal growth rate by the interface pinning method

An essential parameter for crystal growth is the kinetic coefficient given by the proportionality between super-cooling and average growth velocity. Here we show that this coefficient can be computed in a single equilibrium simulation using the interface pinning method where two-phase configurations are stabilized by adding an spring-like bias field coupling to an order-parameter that discriminates between the two phases. Crystal growth is a Smoluchowski process and the crystal growth rate can therefore be computed from the terminal exponential relaxation of the order parameter. The approach is investigated in detail for the Lennard-Jones model. We find that the kinetic coefficient scales as the inverse square-root of temperature along the high temperature part of the melting line. The practical usability of the method is demonstrated by computing the kinetic coefficient of the elements Na, Mg, Al and Si from first principles. It is briefly discussed how a generalized version of the method is an alternative to forward flux sampling methods for computing rates along trajectories of rare events.

preprint2013arXiv

Computing Gibbs free energy differences by interface pinning

We propose an approach for computing the Gibbs free energy difference between phases of a material. The method is based on the determination of the average force acting on interfaces that separate the two phases of interest. This force, which depends on the Gibbs free energy difference between the phases, is computed by applying an external harmonic field that couples to a parameter which specifies the two phases. Validated first for the Lennard-Jones model, we demonstrate the flexibility, efficiency and practical applicability of this approach by computing the melting temperatures of sodium, magnesium, aluminum and silicon at ambient pressure using density functional theory. Excellent agreement with experiment is found for all four elements, except for silicon, for which the melting temperature is, in agreement with previous simulations, seriously underestimated.