Bernd Schmidt

A bending-torsion theory for thin and ultrathin rods as a $Γ$ -limit of atomistic models

The purpose of this note is to establish two continuum theories for the bending and torsion of inextensible rods as $Γ$ -limits of 3D atomistic models. In our derivation we study simultaneous limits of vanishing rod thickness $h$ and interatomic distance $\varepsilon$ . First, we set up a novel theory for ultrathin rods composed of finitely many atomic fibres ( $\varepsilon\sim h$ ), which incorporates surface energy and new discrete terms in the limiting functional. This can be thought of as a contribution to the mechanical modelling of nanowires. Second, we treat the case where $\varepsilon\ll h$ and recover a nonlinear rod model $-$ the modern version of Kirchhoff's rod theory.

A Blake-Zisserman-Kirchhoff theory for plates with soft inclusions

We consider a two phase elastic thin film with soft inclusions subject to bending dominated deformations. The soft (void) phase may comprise asymptotically small droplets within the elastic matrix. We perform a dimension reduction analysis and obtain a novel `Blake-Zisserman-Kirchhoff' model that combines Kirchhoff's classical plate theory with Blake-Zisserman type surface energy contributions at cracks, folds and the boundary of voids.

A continuum model for brittle nanowires derived from an atomistic description by $Γ$ -convergence

Starting from a particle system with short-range interactions, we derive a continuum model for the bending, torsion, and brittle fracture of inextensible rods moving in three-dimensional space. As the number of particles tends to infinity, it is assumed that the rod's thickness is of the same order as the interatomic distance. Fracture energy in the $Γ$ -limit is expressed by an implicit cell formula, which covers different modes of fracture, including (complete) cracks, folds and torsional cracks. In special cases, the cell formula can be significantly simplified. Our approach applies e.g. to atomistic systems with Lennard-Jones-type potentials and is motivated by the research of ceramic nanowires.

Eigendamage: an Eigendeformation model for the variational approximation of cohesive fracture -- a one-dimensional case study

We study an approximation scheme for a variational theory of cohesive fracture in a one-dimensional setting. Here, the energy functional is approximated by a family of functionals depending on a small parameter $0 < \varepsilon \ll 1$ and on two fields: the elastic part of the displacement field and an eigendeformation field that describes the inelastic response of the material beyond the elastic regime. We measure the inelastic contributions of the latter in terms of a non-local energy functional. Our main result shows that, as $\varepsilon \to 0$ , the approximate functionals $Γ$ -converge to a cohesive zone model.

Two phase models for elastic membranes with soft inclusions

We derive an effective membrane theory in the thin film limit within a two phase material model for a specimen consisting of an elastic matrix and soft inclusions. The soft inclusions may lead to the formation of cracks within the elastic matrix and the corresponding limiting models are described by Griffith type fracture energy functionals. We also provide simplified proofs of relaxation results for bulk materials.

preprint2020arXiv

Geometric linearization of theories for incompressible elastic materials and applications

We derive geometrically linearized theories for incompressible materials from nonlinear elasticity theory in the small displacement regime. Our nonlinear stored energy densities may vary on the same (small) length scale as the typical displacements. This allows for applications to multiwell energies as, e.g., encountered in martensitic phases of shape memory alloys and models for nematic elastomers. Under natural assumptions on the asymptotic behavior of such densities we prove Gamma-convergence of the properly rescaled nonlinear energy functionals to the relaxation of an effective model. The resulting limiting theory is geometrically linearized in the sense that it acts on infinitesimal displacements rather than finite deformations, but will in general still have a limiting stored energy density that depends in a nonlinear way on the infinitesimal strains. Our results, in particular, establish a rigorous link of existing finite and infinitesimal theories for incompressible nematic elastomers.

preprint2019arXiv

Surface energy and boundary layers for a chain of atoms at low temperature

We analyze the surface energy and boundary layers for a chain of atoms at low temperature for an interaction potential of Lennard-Jones type. The pressure (stress) is assumed small but positive and bounded away from zero, while the temperature $β^{-1}$ goes to zero. Our main results are: (1) As $β\to \infty$ at fixed positive pressure $p>0$ , the Gibbs measures $μ_β$ and $ν_β$ for infinite chains and semi-infinite chains satisfy path large deviations principles. The rate functions are bulk and surface energy functionals $\overline{\mathcal{E}}_{\mathrm{bulk}}$ and $\overline{\mathcal{E}}_\mathrm{surf}$ . The minimizer of the surface functional corresponds to zero temperature boundary layers. (2) The surface correction to the Gibbs free energy converges to the zero temperature surface energy, characterized with the help of the minimum of $\overline{\mathcal{E}}_\mathrm{surf}$ . (3) The bulk Gibbs measure and Gibbs free energy can be approximated by their Gaussian counterparts. (4) Bounds on the decay of correlations are provided, some of them uniform in $β$ .

preprint2016arXiv

A Griffith-Euler-Bernoulli theory for thin brittle beams derived from nonlinear models in variational fracture mechanics

We study a planar thin brittle beam subject to elastic deformations and cracks described in terms of a nonlinear Griffith energy functional acting on $SBV$ deformations of the beam. In particular we consider the case in which elastic bulk contributions due to finite bending of the beam are comparable to the surface energy which is necessary to completely break the beam into several large pieces. In the limit of vanishing aspect ratio we rigorously derive an effective Griffith-Euler-Bernoulli functional which acts on piecewise $W^{2,2}$ regular curves representing the midline of the beam. The elastic part of this functional is the classical Euler-Bernoulli functional for thin beams in the bending dominated regime in terms of the curve's curvature. In addition there also emerges a fracture term proportional to the number of discontinuities of the curve and its first derivative.

preprint2016arXiv

Existence and Convergence of Solutions of the Boundary Value Problem in Atomistic and Continuum Nonlinear Elasticity Theory

We show existence of solutions for the equations of static atomistic nonlinear elasticity theory on a bounded domain with prescribed boundary values. We also show their convergence to the solutions of continuum nonlinear elasticity theory, with energy density given by the Cauchy-Born rule, as the interatomic distances tend to zero. These results hold for small data close to a stable lattice for general finite range interaction potentials. We also discuss the notion of stability in detail.

preprint2015arXiv

A quantitative geometric rigidity result in SBD

We present a quantitative geometric rigidity estimate for special functions of bounded deformation in a planar setting generalizing the result by Friesecke, James, Müller obtained in nonlinear elasticity theory and the piecewise rigidity result by Chambolle, Giacomini, Ponsiglione for brittle materials which do not store elastic energy. We show that for each deformation there is an associated triple consisting of a partition of the domain, a corresponding piecewise rigid motion being constant on each connected component of the cracked body and a displacement field measuring the distance of the deformation from the piecewise rigid motion. We also present a related estimate in the geometrically linear setting which can be interpreted as a `piecewise Korn-Poincaré inequality'.

preprint2014arXiv

An analysis of crystal cleavage in the passage from atomistic models to continuum theory

We study the behavior of atomistic models in general dimensions under uniaxial tension and investigate the system for critical fracture loads. We rigorously prove that in the discrete-to-continuum limit the minimal energy satisfies a particular cleavage law with quadratic response to small boundary displacements followed by a sharp constant cut-off beyond some critical value. Moreover, we show that the minimal energy is attained by homogeneous elastic configurations in the subcritical case and that beyond critical loading cleavage along specific crystallographic hyperplanes is energetically favorable. In particular, our results apply to mass spring models with full nearest and next-to-nearest pair interactions and provide the limiting minimal energy and minimal configurations.

preprint2014arXiv

On a discrete-to-continuum convergence result for a two dimensional brittle material in the small displacement regime

We consider a two-dimensional atomic mass spring system and show that in the small displacement regime the corresponding discrete energies can be related to a continuum Griffith energy functional in the sense of Gamma-convergence. We also analyze the continuum problem for a rectangular bar under tensile boundary conditions and find that depending on the boundary loading the minimizers are either homogeneous elastic deformations or configurations that are completely cracked generically along a crystallographic line. As applications we discuss cleavage properties of strained crystals and an effective continuum fracture energy for magnets.

Closure and commutability results for Gamma-limits and the geometric linearization and homogenization of multi-well energy functionals

Under a suitable notion of equivalence of integral densities we prove a $Γ$ -closure theorem for integral functionals: The limit of a sequence of $Γ$ -convergent families of such functionals is again a $Γ$ -convergent family. Its $Γ$ -limit is the limit of the $Γ$ -limits of the original problems. This result not only provides a common basic principle for a number of linearization and homogenization results in elasticity theory. It also allows for new applications as we exemplify by proving that geometric linearization and homogenization of multi-well energy functionals commute.

Dynamical arrest of ultracold lattice fermions

We theoretically investigate the thermodynamics of an interacting inhomogeneous two-component Fermi gas in an optical lattice. Motivated by a recent experiment by L. Hackermüller et al., Science, 327, 1621 (2010), we study the effect of the interplay between thermodynamics and strong correlations on the size of the fermionic cloud. We use dynamical mean-field theory to compute the cloud size, which in the experiment shows an anomalous expansion behavior upon increasing attractive interaction. We confirm this qualitative effect but, assuming adiabaticity, we find quantitative agreement only for weak interactions. For strong interactions we observe significant non-equilibrium effects which we attribute to a dynamical arrest of the particles due to increasing correlations.

Ground states of the 2D sticky disc model: fine properties and $N^{3/4}$ law for the deviation from the asymptotic Wulff shape

We investigate ground state configurations for a general finite number $N$ of particles of the Heitmann-Radin sticky disc pair potential model in two dimensions. Exact energy minimizers are shown to exhibit large microscopic fluctuations about the asymptotic Wulff shape which is a regular hexagon: There are arbitrarily large $N$ with ground state configurations deviating from the nearest regular hexagon by a number of $\sim N^{3/4}$ particles. We also prove that for any $N$ and any ground state configuration this deviation is bounded above by $\sim N^{3/4}$ . As a consequence we obtain an exact scaling law for the fluctuations about the asymptotic Wulff shape. In particular, our results give a sharp rate of convergence to the limiting Wulff shape.

On the infinite particle limit in Lagrangian dynamics and convergence of optimal transportation meshfree methods

We consider Lagrangian systems in the limit of infinitely many particles. It is shown that the corresponding discrete action functionals Gamma-converge to a continuum action functional acting on probability measures of particle trajectories. Also the convergence of stationary points of the action is established. Minimizers of the limiting functional and, more generally, limiting distributions of stationary points are investigated and shown to be concentrated on orbits of the Euler-Lagrange flow. We also consider time discretized systems. These results in particular provide a convergence analysis for optimal transportation meshfree methods for the approximation of particle flows by finite discrete Lagrangian dynamics.

Si nano-VINe: Electrical Percolation in Quantum-Confined nc-Si:SiO2 Systems

Nanostructures, especially of silicon, are of paramount importance for new generation solar cell applications. The key material requirements for solar cells are good electrical conductivity, confinement of excitons, and a tunable band-gap that allows for tandem arrangements to absorb a maximally large portion of the solar spectrum. However, to date, existing Si-based nanostructures have addressed these requirements individually. The difficulty is that these requirements are entangled in complex ways, so an effort to improve one can lead to a major trade-offs in the others. In this study, we report on a novel silicon nanostructure, where the key advance is that we have managed to satisfy all the abovementioned requirements simultaneously: We demonstrate that the excitons are confined, the bandgap is tunable, and the electrons can freely travel along the nc-Si network.

preprint2012arXiv

An atomistic-to-continuum analysis of crystal cleavage in a two-dimensional model problem

A two-dimensional atomic mass spring system is investigated for critical fracture loads and its crack path geometry. We rigorously prove that, in the discrete-to-continuum limit, the minimal energy of a crystal under uniaxial tension leads to a universal cleavage law and energy minimizers are either homogeneous elastic deformations or configurations that are completely cracked and do not store elastic energy. Beyond critical loading, the specimen generically cleaves along a unique optimal crystallographic hyperplane. For specific symmetric crystal orientations, however, cleavage might fail. In this case a complete characterization of possible limiting crack geometries is obtained.

preprint2012arXiv

On the passage from atomistic systems to nonlinear elasticity theory

We derive continuum limits of atomistic models in the realm of nonlinear elasticity theory rigorously as the interatomic distances tend to zero. In particular we obtain an integral functional acting on the deformation gradient in the continuum theory which depends on the underlying atomistic interaction potentials and the lattice geometry. The interaction potentials to which our theory applies are general finite range models on multilattices which in particular can also account for multi-pole interactions and bond-angle dependent contributions. Furthermore, we discuss the applicability of the Cauchy-Born rule. Our class of limiting energy densities consists of general quasiconvex functions and the class of linearized limiting energies consistent with the Cauchy-Born rule consists of general quadratic forms not restricted by the Cauchy relations.

preprint2011arXiv

Convergence Analysis of Meshfree Approximation Schemes

This work is concerned with the formulation of a general framework for the analysis of meshfree approximation schemes and with the convergence analysis of the Local Maximum-Entropy (LME) scheme as a particular example. We provide conditions for the convergence in Sobolev spaces of schemes that are n-consistent, in the sense of exactly reproducing polynomials of degree less or equal to n, and whose basis functions are of rapid decay. The convergence of the LME in a locally Sobolev space follows as a direct application of the general theory. The analysis shows that the convergence order is linear in h, a measure of the density of the point set. The analysis also shows how to parameterize the LME scheme for optimal convergence. Because of the convex approximation property of LME, its behavior near the boundary is singular and requires additional analysis. For the particular case of polyhedral domains we show that, away from a small singular part of the boundary, any Sobolev function can be approximated by means of the LME scheme. With the aid of a capacity argument, we further obtain approximation results with truncated LME basis functions in H1 and for spatial dimension d > 2.

preprint2011arXiv

Universal formulae for the limiting elastic energy of membrane networks

We provide universal formulae for the limiting stretching and bending energies of triangulated membrane networks endowed with nearest neighbor bond potentials and cosine-type dihedral angle potentials. The given formulae account for finite elasticity and solve some deficiencies of earlier results for Helfrich-type bending energies, due to shape-dependence and sensitivity to mesh distortion effects of the limiting elastic coefficients. We also provide the entire set of the elastic coefficients characterizing the limiting response of the examined networks, accounting for full bending-stretching coupling. We illustrate the effectiveness of the proposed formulae by way of example, on examining the special cases of cylindrical and spherical networks covered with equilateral triangles, and discussing possible strategies for the experimental characterization of selected elastic moduli.

A sharp version of Ehrenfest's theorem for general self-adjoint operators

We prove the Ehrenfest theorem of quantum mechanics under sharp assumptions on the operators involved.

Discretized vs. continuous models of p-wave interacting fermions in 1D

We present a general mapping between continuous and lattice models of Bose- and Fermi-gases in one dimension, interacting via local two-body interactions. For s-wave interacting bosons we arrive at the Bose-Hubbard model in the weakly interacting, low density regime. The dual problem of p-wave interacting fermions is mapped to the spin-1/2 XXZ model close to the critical point in the highly polarized regime. The mappings are shown to be optimal in the sense that they produce the least error possible for a given discretization length. As an application we examine the ground state of a interacting Fermi gas in a harmonic trap, calculating numerically real-space and momentum-space distributions as well as two-particle correlations. In the analytically known limits the convergence of the results of the lattice model to the continuous one is shown.

Fermionisation dynamics of a strongly interacting 1D Bose gas after an interaction quench

We study the dynamics of a one-dimensional Bose gas after a sudden change of the interaction strength from zero to a finite value using the numerical time-evolving block decimation (TEBD) algorithm. It is shown that despite the integrability of the system, local quantities such as the two-particle correlation $g^{(2)}(x,x)$ attain steady state values in a short characteristic time inversely proportional to the Tonks parameter $γ$ and the square of the density. The asymptotic values are very close to those of a finite temperature grand canonical ensemble with a local temperature corresponding to initial energy and density. Non-local density-density correlations on the other hand approach a steady state on a much larger time scale determined by the finite propagation velocity of oscillatory correlation waves.