Quasi-static Limit Research Articles

Optical bistability in two-dimensional nonlinear composites of coated cylinders with nonlinear core and graded shell

The intrinsic optical bistability of the nonlinear composite consisting of coated cylinders with nonlinear core and graded Drude shell is investigated with self-consistent mean-field approximation. We derive the nonlinear equation of near field for such graded composite in the quasistatic limit. We demonstrate that the bistability threshold and the bistable width are highly dependent on the core radius, the incident frequency, and the graded coefficient of the coated cylinder inclusion. It is found that the optical bistability appears only when the core radius and the incident frequency satisfy some specific conditions. Therefore, the optical bistability for nonlinear composite materials may be optimized by the suitable adjustment of the physical and geometrical parameters. Our results may be helpful for the design of the nonlinear device with appropriate bistability.

Computational nanoplasmonics in the quasistatic limit for biosensing applications.

The phenomenon of localized surface plasmon resonance (LSPR) provides high sensitivity in detecting biomolecules through shifts in resonance frequency when a target is present. Computational studies in this field have used the full Maxwell equations with simplified models of a sensor-analyte system, or they neglected the analyte altogether. In the long-wavelength limit, one can simplify the theory via an electrostatics approximation while adding geometrical detail in the sensor and analytes (at moderate computational cost). This work uses the latter approach, expanding the open-source PyGBe code to compute the extinction cross section of metallic nanoparticles in the presence of any target for sensing. The target molecule is represented by a surface mesh, based on its crystal structure. PyGBe is research software for continuum electrostatics, written in python with computationally expensive parts accelerated on GPU hardware, via PyCUDA. It is also accelerated algorithmically via a treecode that offers O(NlogN) computational complexity. These features allow PyGBe to handle problems with half a million boundary elements or more. In this work, we demonstrate the suitability of PyGBe, extended to compute LSPR response in the electrostatic limit, for biosensing applications. Using a model problem consisting of an isolated silver nanosphere in an electric field, our results show grid convergence as 1/N, and accurate computation of the extinction cross section as a function of wavelength (compared with an analytical solution). For a model of a sensor-analyte system, consisting of a spherical silver nanoparticle and a set of bovine serum albumin (BSA) proteins, our results again obtain grid convergence as 1/N (with respect to the Richardson extrapolated value). Computing the LSPR response as a function of wavelength in the presence of BSA proteins captures a redshift of 0.5 nm in the resonance frequency due to the presence of the analytes at 1-nm distance. The final result is a sensitivity study of the biosensor model, obtaining the shift in resonance frequency for various distances between the proteins and the nanoparticle. All results in this paper are fully reproducible, and we have deposited in archival data repositories all the materials needed to run the computations again and recreate the figures. PyGBe is open source under a permissive license and openly developed. Documentation is available at http://pygbe.github.io/pygbe/docs/.

Apparent violations of the second law in two-level systems.

We examine two-level systems which undergo an instantaneous change in the energy level differences due to the application of an external protocol. In particular the behavior of the work distribution as well as the probability of second-law violations as the thermodynamic limit, and separately the quasistatic limit, is approached are investigated.

On the Approximation of Quasistatic Evolutions for the Debonding of a Thin Film via Vanishing Inertia and Viscosity

In this paper we contribute to studying the issue of quasistatic limit in the context of Griffith's theory by investigating a one-dimensional debonding model. It describes the evolution of a thin film partially glued to a rigid substrate and subjected to a vertical loading. Taking friction into account and under suitable assumptions on the toughness of the glue, we prove that, in contrast to what happens in the undamped case, dynamic solutions converge to the quasistatic one when inertia and viscosity go to zero, except for a possible discontinuity at the initial time. We then characterise the size of the jump by means of an asymptotic analysis of the debonding front.

Assessment of hardening due to non-coherent precipitates in tungsten-rhenium alloys at the atomic scale

In metallurgical applications, precipitation strengthening is of great technological importance to engineer materials with the required strength. While precipitation hardening is essential for many applications involving operation at elevated temperatures, its subsequent embrittlement can be a showstopper for the overall performance of a component. In the nuclear industry, irradiation-induced/enhanced precipitation and the resulting embrittlement often limit the lifetime of components. In fusion applications, tungsten (W) based alloys are known to harden and embrittle as a result of irradiation-assisted transmutation to rhenium (Re) and its subsequent precipitation into non-coherent precipitates. Hence, a fundamental understanding of the interaction of dislocations with non-coherent precipitates is of great interest. In the present work, the interaction of dislocations with non-coherent Re-rich σ, χ and hcp phase precipitates embedded in a bcc W matrix is assessed. Large-scale atomistic simulations are performed to clarify the interaction mechanisms and derive the obstacle strength of the precipitates in the quasi-static limit. Thereby the impact of precipitate shape, size, interspacing and composition is assessed. Based on those results, an analytical model to predict precipitation hardening of σ, χ and hcp phase particles in bcc W is proposed and compared to available experimental data from mechanical tests on irradiated materials.

Effect of Pinning on the Yielding Transition of Amorphous Solids.

Using numerical simulations, we have studied the yielding response, in the athermal quasistatic limit, of a model amorphous material having inclusions in the form of randomly pinned particles. We show that, with increasing pinning concentration, the plastic activity becomes more spatially localized, resulting in smaller stress drops, and a corresponding increase in the magnitude of strain where yielding occurs. We demonstrate that, unlike the spatially heterogeneous and avalanche led yielding in the case of the unpinned glass, for the case of large pinning concentration, yielding takes place via a spatially homogeneous proliferation of localized events.

Equivalence principle on cosmological backgrounds in scalar–tensor theories

We study the extent up to which the equivalence principle is obeyed in models of modified gravity and dark energy involving a single scalar degree of freedom. We focus on the effective field theories of dark energy describing the late time acceleration in the presence of ordinary matter species. In their covariant form these coincide with the Horndeski theories on a cosmological background with a slowly varying Hubble rate and a time-dependent scalar field. We show that in the case of an exactly de Sitter universe both the weak and strong equivalence principles hold. This occurs due to the combination of the shift symmetry of the scalar and the time translational invariance of de Sitter space. When generalized to Friedmann–Robertson–Walker cosmologies (FRW) we show that the weak equivalence principle is obeyed for test particles and extended objects in the quasi-static subhorizon limit. We do this by studying the field created by an extended object moving in an FRW background as well as its equation of motion in an external gravitational field. We also estimate the corrections to the geodesic equation of the extended object due to the approximations made.

Dynamic enlargement of a hole in a sheet: Crater formation and propagation of cylindrical shock waves

Predicting the shape of a crater formed by high velocity impact is of interest in several fields. It can aid in design of more efficient protective structures, in forensic analysis of bullet holes, and in understanding the effects of meteorite impact in both space systems and in extreme geological events. In this paper we present, for the first time, a complete theoretical solution of the dynamic plane-stress problem. We consider the steady-state expansion of a cylindrical hole in a strain hardening elastoplastic sheet and find that a self-similar field emerges if the ‘specific cavitation energy’ is constant. It is shown that at the quasistatic limit this solution reduces to available classical solutions, while at high expansion velocities shock waves can appear. Investigation of the constitutive sensitivities of the expansion field is conducted and compared with available results for the spherical field which is commonly applied to predict resistance to high velocity penetration. It is shown that shock waves appear at significantly lower expansion velocities, in the plane-stress deformation pattern, for which material compressibility is found to have a negligible effect. This insensitivity can be taken advantage of in the future for design of light weight protective layers by incorporating porosity.

Power and efficiency of a thermal engine with a coherent bath.

We consider a quantum engine driven by repeated weak interactions with a heat bath of identical three-level atoms. This model was first introduced by Scully et al. [Science 299, 862 (2003)SCIEAS0036-807510.1126/science.1078955], who showed that coherence between the energy-degenerate ground states serves as a thermodynamic resource that allows operation of a thermal cycle with a coherence-dependent thermalization temperature. We consider a similar engine out of the quasistatic limit and find that the ground-state coherence also determines the rate of thermalization, therefore increasing the output power and the engine efficiency only when the thermalization temperature is reduced; revealing a more nuanced perspective of coherence as a resource. This allows us to optimize the output power by adjusting the coherence and relative stroke durations.

Slender-body theory for plasmonic resonance.

We develop a slender-body theory for plasmonic resonance of slender metallic nanoparticles, focusing on a general class of axisymmetric geometries with locally paraboloidal tips. We adopt a modal approach where one first solves the plasmonic eigenvalue problem, a geometric spectral problem which governs the surface-plasmon modes of the particle; then, the latter modes are used, in conjunction with spectral-decomposition, to analyse localized-surface-plasmon resonance in the quasi-static limit. We show that the permittivity eigenvalues of the axisymmetric modes are strongly singular in the slenderness parameter, implying widely tunable, high-quality-factor, resonances in the near-infrared regime. For that family of modes, we use matched asymptotics to derive an effective eigenvalue problem, a singular non-local Sturm-Liouville problem, where the lumped one-dimensional eigenfunctions represent axial voltage profiles (or charge line densities). We solve the effective eigenvalue problem in closed form for a prolate spheroid and numerically, by expanding the eigenfunctions in Legendre polynomials, for arbitrarily shaped particles. We apply the theory to plane-wave illumination in order to elucidate the excitation of multiple resonances in the case of non-spheroidal particles.

Stochastic work extraction in a colloidal heat engine in the presence of colored noise

From synthetic active devices such as self-propelling Janus colloids to micro-organisms like bacteria, micro-algae, living cells in tissues, active fluctuations are ubiquitous. Thermodynamics of small systems involving thermal as well as active fluctuations are of immense importance. They can be employed to extract thermodynamic work. Here we propose a simple model system that can produce thermodynamic work exploiting active fluctuations. We consider a Brownian particle, trapped by an externally controlled harmonic confinement that contracts and expands time-periodically by modulating its spring constant e.g. optical tweezers. The system produces work by being alternately connected to two baths one passive and other active. The active bath provides exponentially time-correlated noise to the particle, that breaks the fluctuation dissipation relation. The average efficiency of the system is calculated exactly in quasistatic limit. Nonquasistatic regime is explored by numerics. Comparing with its passive counterpart, we also show that the active micro heat engine can be more efficient depending on the chosen parameter space. We also believe that our model can be realised experimentally with the help of bacterial baths.

Effective drag of a rod in fluid-saturated granular beds

We measure the drag encountered by a vertically oriented rod moving across a sedimented granular bed immersed in a fluid under steady-state conditions. At low rod speeds, the presence of the fluid leads to a lower drag because of buoyancy, whereas a significantly higher drag is observed with increasing speeds. The drag as a function of the depth is observed to decrease from being quadratic at low speeds to appearing more linear at higher speeds. By scaling the drag with the average weight of the grains acting on the rod, we obtain the effective friction μ_{e} encountered over six orders of magnitude of speeds. While a constant μ_{e} is found when the grain size, rod depth, and fluid viscosity are varied at low speeds, a systematic increase is observed as the speed is increased. We analyze μ_{e} in terms of the inertial number I and viscous number J to understand the relative importance of inertia and viscous forces, respectively. For sufficiently high fluid viscosities, we find that the effect of varying the speed, depth, and viscosity can be described by the empirical function μ_{e}=μ_{o}+kJ^{n}, where μ_{o} is the effective friction measured in the quasistatic limit, and k and n are material constants. The drag is then analyzed in terms of the effective viscosity η_{e} and found to decrease systematically as a function of J. We further show that η_{e} as a function of J is directly proportional to the fluid viscosity and the μ_{e} encountered by the rod.

Vainshtein regime in scalar-tensor gravity: Constraints on degenerate higher-order scalar-tensor theories

We study the screening mechanism in the most general scalar-tensor theories that leave gravitational waves unaffected and are thus compatible with recent LIGO/Virgo observations. Using the effective field theory of dark energy approach, we consider the general action for perturbations beyond linear order, focussing on the quasi-static limit. When restricting to the subclass of theories that satisfy the gravitational wave constraints, the fully nonlinear effective Lagrangian contains only three independent parameters. One of these, $\beta_1$, is uniquely present in degenerate higher-order theories. We compute the two gravitational potentials for a spherically symmetric matter source and we find that for $\beta_1 \ge 0$ they decrease as the inverse of the distance, as in standard gravity, while the case $\beta_1 < 0$ is ruled out. For $\beta_1 > 0$, the two potentials differ and their gravitational constants are not the same on the inside and outside of the body. Generically, the bound on anomalous light bending in the Solar System constrains $\beta_1 \lesssim 10^{-5}$. Standard gravity can be recovered outside the body by tuning the parameters of the model, in which case $\beta_1 \lesssim 10^{-2}$ from the Hulse-Taylor pulsar.

Electrostatic T-matrix for a torus on bases of toroidal and spherical harmonics

Semi-analytic expressions for the static limit of the T-matrix for electromagnetic scattering are derived for a circular torus, expressed in bases of both toroidal and spherical harmonics. The scattering problem for an arbitrary static excitation is solved using toroidal harmonics and the extended boundary condition method to obtain analytic expressions for auxiliary Q and P-matrices, from which the T-matrix is given by their division. By applying the basis transformations between toroidal and spherical harmonics, the quasi-static limit of the T-matrix block for electric multipole coupling is obtained. For the toroidal geometry there are two similar T-matrices on a spherical basis, for computing the scattered field both near the origin and in the far field. Static limits of the optical cross-sections are computed, and analytic expressions for the limit of a thin ring are derived.

A variational approach to the quasistatic limit of viscous dynamic evolutions in finite dimension

In this paper we study the vanishing inertia and viscosity limit of a second order system set in an Euclidean space, driven by a possibly nonconvex time-dependent potential satisfying very general assumptions. By means of a variational approach, we show that the solutions of the singularly perturbed problem converge to a curve of stationary points of the energy and characterize the behavior of the limit evolution at jump times. At those times, the left and right limits of the evolution are connected by a finite number of heteroclinic solutions to the unscaled equation.

Application of shakedown theory in track substructure design

Tracks require high performance of earth structures against excessive residual settlement under repeated train loads. In the current design standards, the quality of soil materials and compaction control are strictly prescribed based on empirical knowledge and scale tests; however, construction costs become relatively high accordingly. A more rational design method based on shakedown theory is presented here. This method makes use of the lower-bound shakedown theorem and models the soils as Mohr–Coulomb materials. By calculating the analytical load-induced dynamic elastic stress field in a three-dimensional half-space and introducing a self-equilibrated residual stress field, the maximum magnitude of surface pressure on the soil foundation (i.e. shakedown limit) against long-term residual settlement is calculated through an optimisation programme. It is found the shakedown limit is dependent on the ratio of train velocity to shear wave velocity for a given earth structure. The maximum admissible train loads for the cases with different design velocities can be predicted by modifying the quasi-static shakedown limit with an attenuation factor.

Cosmology of surviving Horndeski theory: The road ahead

In the context of the effective field theory of dark energy (EFT) we perform agnostic explorations of Horndeski gravity. We choose two parametrizations for the free EFT functions, namely a power law and a dark energy density-like behaviour on a non trivial Chevallier-Polarski-Linder background. We restrict our analysis to those EFT functions which do not modify the speed of propagation of gravitational waves. Among those, we prove that one specific function cannot be constrained by data, since its contribution to the observables is below the cosmic variance, although we show it has a relevant role in defining the viable parameter space. We place constraints on the parameters of these models combining measurements from present day cosmological datasets and we prove that the next generation galaxy surveys can improve such constraints by one order of magnitude. We then verify the validity of the quasi-static limit within the sound horizon of the dark field, by looking at the phenomenological functions $\mu$ and $\Sigma$, associated respectively to clustering and lensing potentials. Furthermore, we notice up to $5\%$ deviations in $\mu, \Sigma$ with respect to General Relativity at scales smaller than the Compton one. For the chosen parametrizations and in the quasi-static limit, future constraints on $\mu$ and $\Sigma$ can reach the $1\%$ level and will allow us to discriminate between certain models at more than $3\sigma$, provided the present best-fit values remain.

Optical response of Luttinger semimetals in the normal and superconducting states

We investigate the optical response properties of three-dimensional Luttinger semimetals with the Fermi energy close to a quadratic band touching point. In particular, in order to address recent experiments on the spectroscopy of Pyrochlore Iridates and half-Heusler superconductors, we derive expressions for the optical conductivity in both the normal and general superconducting states in the linear response regime within the random phase approximation. The response functions can be decomposed into contributions from intraband and interband transitions, the latter comprising a genuine signature of the quadratic band touching point. We demonstrate the importance of interband transitions in the optical response in the normal state both in the homogeneous and quasi-static limit. Our analysis reveals a factorization property of the homogeneous conductivity in the spatially anisotropic case and the divergence of the conductivity for strong spatial anisotropy. In the quasi-static limit, the response is dominated by interband transitions and significantly different from systems with a single parabolic band. As an applications of the formalism in the superconducting state we compute the optical conductivity and superfluid density for the s-wave singlet superconducting case for both finite and vanishing chemical potential.

A Stochastic Bulk Damage Model Based on Mohr-Coulomb Failure Criterion for Dynamic Rock Fracture

We present a stochastic bulk damage model for rock fracture. The decomposition of strain or stress tensor to its negative and positive parts is often used to drive damage and evaluate the effective stress tensor. However, they typically fail to correctly model rock fracture in compression. We propose a damage force model based on the Mohr-Coulomb failure criterion and an effective stress relation that remedy this problem. An evolution equation specifies the rate at which damage tends to its quasi-static limit. The relaxation time of the model introduces an intrinsic length scale for dynamic fracture and addresses the mesh sensitivity problem of earlier damage models. The ordinary differential form of the damage equation makes this remedy quite simple and enables capturing the loading rate sensitivity of strain-stress response. The asynchronous Spacetime Discontinuous Galerkin (aSDG) method is used for macroscopic simulations. To study the effect of rock inhomogeneity, the Karhunen-Loeve method is used to realize random fields for rock cohesion. It is shown that inhomogeneity greatly differentiates fracture patterns from those of a homogeneous rock, including the location of zones with maximum damage. Moreover, as the correlation length of the random field decreases, fracture patterns resemble angled-cracks observed in compressive rock fracture.

Escape dynamics of a forced‐damped classical particle in an infinite‐range potential well

AbstractThe paper is devoted to analysis of the escape of periodically forced and damped particle from one‐dimensional potential well. The particle is initially at rest, and the forcing is switched on at a certain time instance. The present work is an extension of previous results, obtained in Hamiltonian setting, for much more realistic case with viscous damping. Assuming primary 1:1 resonance, one can consider the problem in terms of averaged transient dynamics. It turns out that, similar to the undamped case, the escape process can be reliably described in terms of topology of special trajectories on the resonant manifold. A theoretical prediction to the minimal force required for the escape as function of the excitation frequency for various damping coefficients is provided. In the explored frequency range, numeric simulations are in complete qualitative and reasonable quantitative agreement with the theoretical predictions except for small frequencies under 0.3. These discrepancies are related to quasistatic asymptotic limit of the considered model.