skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Heterogeneous diffusion in an harmonic potential: the role of the interpretation
Abstract Diffusion in heterogeneous energy and diffusivity landscapes is widespread in biological systems. However, solving the Langevin equation in such environments introduces ambiguity due to the interpretation parameter $$\alpha$$, which depends on the underlying physics and can take values in the range $$0<\alpha<1$$. The typical interpretations are It\^o ($$\alpha=0$$), Stratonovich ($$\alpha=1/2$$), and Hanggi-Klimontovich ($$\alpha=1$$). Here, we analyse the motion of a particle in an harmonic potential---modelled as an Ornstein-Uhlenbeck process---with diffusivity that varies in space. Our focus is on two-phase systems with a discontinuity in environmental properties at $x=0$. We derive the probability density of the particle position for the process, and consider two paradigmatic situations. In the first one, the damping coefficient remains constant, and fluctuation-dissipation relations are not satisfied. In the second one, these relations are enforced, leading to a position-dependent damping coefficient. In both cases, we provide solutions as a function of the interpretation parameter $$\alpha$$, with particular attention to the It\^o, Stratonovich, and Hanggi-Klimontovich interpretations, revealing fundamentally different behaviours, in particular with respect to an interface located at the potential minimum.  more » « less
Award ID(s):
2102832
PAR ID:
10596722
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
New Journal of Physics
ISSN:
1367-2630
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The distribution of natural frequencies of the Euler–Bernoulli beam subject to fully non-dissipative boundary conditions is investigated. The beam is clamped at the left end and equipped with a 4-parameter ($$\alpha ,\beta ,k_1,k_2$$) linear boundary feedback law at the right end. The $$2 \times 2$$ boundary feedback matrix relates the control input (a vector of velocity and its spatial derivative at the right end), to the output (a vector of shear and moment at the right end). The initial boundary value problem describing the dynamics of the beam has been reduced to the first order in time evolution equation in the state Hilbert space equipped with the energy norm. The dynamics generator has a purely discrete spectrum (the vibrational modes) denoted by $$\{\nu _n\}_{n\in \mathbb {Z}^{\prime}}$$. The role of the control parameters is examined and the following results have been proven: (i) when $$\beta \neq 0$$, the set of vibrational modes is asymptotically close to the vertical line on the complex $$\nu$$-plane given by the equation $$\Re \nu = \alpha + (1-k_1k_2)/\beta$$; (ii) when $$\beta = 0$$ and the parameter $$K = (1-k_1 k_2)/(k_1+k_2)$$ is such that $$\left |K\right |\neq 1$$ then the following relations are valid: $$\Re (\nu _n/n) = O\left (1\right )$$ and $$\Im (\nu _n/n^2) = O\left (1\right )$$ as $$\left |n\right |\to \infty$$; (iii) when $$\beta =0$$, $|K| = 1$, and $$\alpha = 0$$, then the following relations are valid: $$\Re (\nu _n/n^2) = O\left (1\right )$$ and $$\Im (\nu _n/n) = O\left (1\right )$$ as $$\left |n\right |\to \infty$$; (iv) when $$\beta =0$$, $|K| = 1$, and $$\alpha>0$$, then the following relations are valid: $$\Re (\nu _n/\ln \left |n\right |) = O\left (1\right )$$ and $$\Im (\nu _n/n^2) = O\left (1\right )$$ as $$\left |n\right |\to \infty$$. 
    more » « less
  2. null (Ed.)
    The distribution of natural frequencies of the Euler–Bernoulli beam resting on elastic foundation and subject to an axial force in the presence of several damping mechanisms is investigated. The damping mechanisms are: ( i ) an external or viscous damping with damping coefficient ( − a 0 ( x )), ( ii ) a damping proportional to the bending rate with the damping coefficient a 1 ( x ). The beam is clamped at the left end and equipped with a four-parameter (α, β, κ 1 , κ 2 ) linear boundary feedback law at the right end. The 2 × 2 boundary feedback matrix relates the control input (a vector of velocity and its spacial derivative at the right end) to the output (a vector of shear and moment at the right end). The initial boundary value problem describing the dynamics of the beam has been reduced to the first order in time evolution equation in the state Hilbert space of the system. The dynamics generator has a purely discrete spectrum (the vibrational modes). Explicit asymptotic formula for the eigenvalues as the number of an eigenvalue tends to infinity have been obtained. It is shown that the boundary control parameters and the distributed damping play different roles in the asymptotical formulas for the eigenvalues of the dynamics generator. Namely, the damping coefficient a 1 and the boundary controls κ 1 and κ 2 enter the leading asymptotical term explicitly, while damping coefficient a 0 appears in the lower order terms. 
    more » « less
  3. Abstract When compressed, certain lattices undergo phase transitions that may allow nuclei to gain sig- nificant kinetic energy. To explore the dynamics of this phenomenon, we develop a methodology to study Coulomb coupled N-body systems constrained to a sphere, as in the Thomson problem. We initialize N total Boron nuclei as point particles on the surface of the sphere, allowing them to equilibrate via Coulomb scattering with a viscous damping term. To simulate a phase transition, we remove Nrm particles, forcing the system to rearrange into a new equilibrium. With this model, we consider the Thomson problem as a dynamical system, providing a framework to explore how non-zero temperature affects structural imperfections in Thomson minima. We develop a scaling relation for the average peak kinetic energy attained by a single particle as a function of N and Nrm. For certain values of N , we find an order of magnitude energy gain when increasing Nrm from 1 to 6. The model may help to design a lattice that maximizes the energy output. 
    more » « less
  4. Abstract We theoretically and computationally investigate the cooling of antihydrogen, ¯H, using optical molasses cooling. This updates the results in Ref. [1] to the current capabilities of the ALPHA experiment. Through Monte Carlo simulation, we show that ¯Hs do not give the standard cooling even in an ideal optical molasses because of their small mass and large transition frequency. For optical molasses cooling in the ALPHA trap, the photons are constrained to travel in one direction only. It is only through the phase space mixing in the trap that cooling in all directions can be achieved. We explore the nontrivial role that laser intensity plays in the cooling. We also investigate the possibility for simultaneously cooling atoms in either of the trapped ground states. 
    more » « less
  5. A one-parameter family of coupled flows depending on a parameter $$\kappa>0$$ is introduced which reduces when $$\kappa=1$$ to the coupled flow of a metric $$\omega$$ with a $(1,1)$-form $$\alpha$$ due recently to Y. Li, Y. Yuan, and Y. Zhang. It is shown in particular that, for $$\kappa\not=1$$, estimates for derivatives of all orders would follow from $C^0$ estimates for $$\omega$$ and $$\alpha$$. Together with the monotonicity of suitably adapted energy functionals, this can be applied to establish the convergence of the flow in some situations, including on Riemann surfaces. Very little is known as yet about the monotonicity and convergence of flows in presence of couplings, and conditions such as $$\kappa\not=1$$ seem new and may be useful in the future. 
    more » « less