-
Short-time Fokker-Planck propagator beyond the Gaussian approximation
Authors:
Julian Kappler
Abstract:
We present a perturbation approach to calculate the short-time propagator, or transition density, of the one-dimensional Fokker-Planck equation, to in principle arbitrary order in the time increment. Our approach preserves probability exactly and allows us to evaluate expectation values of analytical observables to in principle arbitrary accuracy; to showcase this, we derive perturbation expansion…
▽ More
We present a perturbation approach to calculate the short-time propagator, or transition density, of the one-dimensional Fokker-Planck equation, to in principle arbitrary order in the time increment. Our approach preserves probability exactly and allows us to evaluate expectation values of analytical observables to in principle arbitrary accuracy; to showcase this, we derive perturbation expansions for the moments of the spatial increment, the finite-time Kramers-Moyal coefficients, and the mean medium entropy production rate. For an explicit multiplicative-noise system with available analytical solution, we validate all our perturbative results. Throughout, we compare our perturbative results to those obtained from the widely used Gaussian approximation of the short-time propagator; we demonstrate that this Gaussian propagator leads to errors that can be many orders of magnitude larger than those resulting from our perturbation approach. Potential applications of our results include parametrizing diffusive stochastic dynamics from observed time series, and sampling path spaces of stochastic trajectories numerically.
△ Less
Submitted 28 May, 2024;
originally announced May 2024.
-
Resolution dependence of most probable pathways with state-dependent diffusivity
Authors:
Alice L. Thorneywork,
Jannes Gladrow,
Ulrich F. Keyser,
Michael E. Cates,
Ronojoy Adhikari,
Julian Kappler
Abstract:
Recent experiments have probed the relative likelihoods of trajectories in stochastic systems by observing survival probabilities within a tube of radius $R$ in spacetime. We measure such probabilities here for a colloidal particle in a corrugated channel, corresponding to a bistable potential with state-dependent diffusivity. In contrast to previous findings for state-independent noise, we find t…
▽ More
Recent experiments have probed the relative likelihoods of trajectories in stochastic systems by observing survival probabilities within a tube of radius $R$ in spacetime. We measure such probabilities here for a colloidal particle in a corrugated channel, corresponding to a bistable potential with state-dependent diffusivity. In contrast to previous findings for state-independent noise, we find that the most probable pathway changes qualitatively as the tube radius $R$ is altered. We explain this by computing the survival probabilities predicted by overdamped Langevin dynamics. At high enough resolution (small enough $R$), survival probabilities depend solely on diffusivity variations, independent of deterministic forces; finite $R$ corrections yield a generalization of the Onsager-Machlup action. As corollary, ratios of survival probabilities are singular as $R \to 0$, but become regular, and described by the classical Onsager-Machlup action, only in the special case of state-independent noise.
△ Less
Submitted 2 February, 2024;
originally announced February 2024.
-
Barrier-crossing times for different non-Markovian friction in well and barrier $-$ A numerical study
Authors:
Florian N. Brünig,
Roland R. Netz,
Julian Kappler
Abstract:
We introduce a generalized Langevin model system for different non-Markovian effects in the well and barrier regions of a potential, and use it to numerically study the dependence of the barrier-crossing time. In the appropriate limits, our model interpolates between the theoretical barrier-crossing-time predictions by Grote and Hynes (GH), as well as by Pollak et al., which for a single barrier m…
▽ More
We introduce a generalized Langevin model system for different non-Markovian effects in the well and barrier regions of a potential, and use it to numerically study the dependence of the barrier-crossing time. In the appropriate limits, our model interpolates between the theoretical barrier-crossing-time predictions by Grote and Hynes (GH), as well as by Pollak et al., which for a single barrier memory time can differ by several orders of magnitude. Our model furthermore allows to test an analytic rate theory for space-inhomogeneous memory, which disagrees with our numerical results in the long well-memory regime. In this regime, we find that short barrier memory decreases the barrier-crossing time as compared to long barrier memory. This is in contrast with the short well-memory regime, where both our numerical results and GH theory predict an acceleration of the barrier crossing time with increasing barrier memory time. Both effects, the `Markovian-barrier acceleration' and GH `non-Markovian-barrier acceleration' can be understood from a committor analysis. Our model combines finite relaxation times of orthogonal degrees of freedom with a space-inhomogeneous coupling to such degrees, and represents a step towards more realistic modeling of physical reaction coordinates.
△ Less
Submitted 11 July, 2022;
originally announced July 2022.
-
Diffusivity dependence of the transition path ensemble
Authors:
Lukas Kikuchi,
Ronojoy Adhikari,
Julian Kappler
Abstract:
Transition pathways of stochastic dynamical systems are typically approximated by instantons. Here we show, using a dynamical system containing two competing pathways, that at low-to-intermediate temperatures, instantons can fail to capture the most likely transition pathways. We construct an approximation which includes fluctuations around the instanton and, by comparing with the results of an ac…
▽ More
Transition pathways of stochastic dynamical systems are typically approximated by instantons. Here we show, using a dynamical system containing two competing pathways, that at low-to-intermediate temperatures, instantons can fail to capture the most likely transition pathways. We construct an approximation which includes fluctuations around the instanton and, by comparing with the results of an accurate and efficient path-space Monte Carlo sampling method, find this approximation to hold for a wide range of temperatures. Our work delimits the applicability of large deviation theory and provides methods to probe these limits numerically.
△ Less
Submitted 4 August, 2022; v1 submitted 24 March, 2022;
originally announced March 2022.
-
Time-dependent friction effects on vibrational infrared frequencies and line shapes of liquid water
Authors:
Florian N. Brünig,
Otto Geburtig,
Alexander von Canal,
Julian Kappler,
Roland R. Netz
Abstract:
From ab initio simulations of liquid water, the time-dependent friction functions and time-averaged non-linear effective bond potentials for the OH stretch and HOH bend vibrations are extracted. The obtained friction exhibits adiabatic contributions at and below the vibrational time scales, but also much slower non-adiabatic contributions, reflecting homogeneous and inhomogeneous line broadening m…
▽ More
From ab initio simulations of liquid water, the time-dependent friction functions and time-averaged non-linear effective bond potentials for the OH stretch and HOH bend vibrations are extracted. The obtained friction exhibits adiabatic contributions at and below the vibrational time scales, but also much slower non-adiabatic contributions, reflecting homogeneous and inhomogeneous line broadening mechanisms, respectively. Compared to the gas phase, hydration softens both stretch and bend potentials, which by itself would lead to a red-shift of the corresponding vibrational bands. In contrast, non-adiabatic friction contributions cause a spectral blue shift. For the stretch mode, the potential effect dominates and thus a significant red shift when going from gas to the liquid phase results. For the bend mode, potential and non-adiabatic friction effects are of comparable magnitude, so that a slight blue shift results, in agreement with well-known but puzzling experimental findings. The observed line broadening is shown to be roughly equally caused by adiabatic and non-adiabatic friction contributions for both, the stretch and bend modes in liquid water. Thus, the understanding of infrared vibrational frequencies and line shapes is considerably advanced by the quantitative analysis of the time-dependent friction that acts on vibrational modes in liquids.
△ Less
Submitted 5 January, 2022;
originally announced January 2022.
-
Efficient Bayesian inference of fully stochastic epidemiological models with applications to COVID-19
Authors:
Yuting I. Li,
Günther Turk,
Paul B. Rohrbach,
Patrick Pietzonka,
Julian Kappler,
Rajesh Singh,
Jakub Dolezal,
Timothy Ekeh,
Lukas Kikuchi,
Joseph D. Peterson,
Austen Bolitho,
Hideki Kobayashi,
Michael E. Cates,
R. Adhikari,
Robert L. Jack
Abstract:
Epidemiological forecasts are beset by uncertainties about the underlying epidemiological processes, and the surveillance process through which data are acquired. We present a Bayesian inference methodology that quantifies these uncertainties, for epidemics that are modelled by (possibly) non-stationary, continuous-time, Markov population processes. The efficiency of the method derives from a func…
▽ More
Epidemiological forecasts are beset by uncertainties about the underlying epidemiological processes, and the surveillance process through which data are acquired. We present a Bayesian inference methodology that quantifies these uncertainties, for epidemics that are modelled by (possibly) non-stationary, continuous-time, Markov population processes. The efficiency of the method derives from a functional central limit theorem approximation of the likelihood, valid for large populations. We demonstrate the methodology by analysing the early stages of the COVID-19 pandemic in the UK, based on age-structured data for the number of deaths. This includes maximum a posteriori estimates, MCMC sampling of the posterior, computation of the model evidence, and the determination of parameter sensitivities via the Fisher information matrix. Our methodology is implemented in PyRoss, an open-source platform for analysis of epidemiological compartment models.
△ Less
Submitted 13 August, 2021; v1 submitted 21 October, 2020;
originally announced October 2020.
-
Sojourn probabilities in tubes and pathwise irreversibility for Itô processes
Authors:
Julian Kappler,
Michael E. Cates,
Ronojoy Adhikari
Abstract:
The sojourn probability of an Itô diffusion process, i.e. its probability to remain in the tubular neighborhood of a smooth path, is a central quantity in the study of path probabilities. For $N$-dimensional Itô processes with state-dependent full-rank diffusion tensor, we derive a general expression for the sojourn probability in tubes whose radii are small but finite, and fixed by the metric of…
▽ More
The sojourn probability of an Itô diffusion process, i.e. its probability to remain in the tubular neighborhood of a smooth path, is a central quantity in the study of path probabilities. For $N$-dimensional Itô processes with state-dependent full-rank diffusion tensor, we derive a general expression for the sojourn probability in tubes whose radii are small but finite, and fixed by the metric of the ambient Euclidean space. The central quantity in our study is the exit rate at which trajectories leave the tube for the first time. This has an interpretation as a Lagrangian and can be measured directly in experiment, unlike previously defined sojourn probabilities which depend on prior knowledge of the state-dependent diffusivity. We find that while in the limit of vanishing tube radius the ratio of sojourn probabilities for a pair of distinct paths is in general divergent, the same for a path and its time-reversal is always convergent and finite. This provides a pathwise definition of irreversibility for Itô processes that is agnostic to the state-dependence of the diffusivity. For one-dimensional systems we derive an explicit expression for our Lagrangian in terms of the drift and diffusivity, and find that our result differs from previously reported multiplicative-noise Lagrangians. We confirm our result by comparing to numerical simulations, and relate our theory to the Stratonovich Lagrangian for multiplicative noise. For one-dimensional systems, we discuss under which conditions the vanishing-radius limiting ratio of sojourn probabilities for a pair of forward and backward paths recovers the established pathwise entropy production. Finally, we demonstrate for our one-dimensional example system that the most probable tube for a barrier crossing depends sensitively on the tube radius, and hence on the tolerated amount of fluctuations around the smooth reference path.
△ Less
Submitted 2 February, 2024; v1 submitted 8 September, 2020;
originally announced September 2020.
-
Multiscale Topology Characterises Dynamic Tumour Vascular Networks
Authors:
Bernadette J. Stolz,
Jakob Kaeppler,
Bostjan Markelc,
Franziska Mech,
Florian Lipsmeier,
Ruth J. Muschel,
Helen M. Byrne,
Heather A. Harrington
Abstract:
Advances in imaging techniques enable high resolution 3D visualisation of vascular networks over time and reveal abnormal structural features such as twists and loops, and their quantification is an active area of research. Here we showcase how topological data analysis (TDA), the mathematical field that studies `shape' of data, can characterise the geometric, spatial and temporal organisation of…
▽ More
Advances in imaging techniques enable high resolution 3D visualisation of vascular networks over time and reveal abnormal structural features such as twists and loops, and their quantification is an active area of research. Here we showcase how topological data analysis (TDA), the mathematical field that studies `shape' of data, can characterise the geometric, spatial and temporal organisation of vascular networks. We propose two topological lenses to study vasculature, which capture inherent multi-scale features and vessel connectivity, and surpass the single scale analysis of existing methods. We analyse images collected using intravital and ultramicroscopy modalities and quantify spatio-temporal variation of twists, loops, and avascular regions (voids) in 3D vascular networks. This topological approach validates and quantifies known qualitative trends such as dynamic changes in tortuosity and loops in response to antibodies that modulate vessel sprouting; furthermore, it quantifies the effect of radiotherapy on vessel architecture.
△ Less
Submitted 26 April, 2022; v1 submitted 19 August, 2020;
originally announced August 2020.
-
Measurement of irreversibility and entropy production via the tubular ensemble
Authors:
Julian Kappler,
Ronojoy Adhikari
Abstract:
The appealing theoretical measure of irreversibility in a stochastic process, as the ratio of the probabilities of a trajectory and its time reversal, cannot be accessed directly in experiment since the probability of a single trajectory is zero. We regularize this definition by considering, instead, the limiting ratio of probabilities for trajectories to remain in the tubular neighborhood of a sm…
▽ More
The appealing theoretical measure of irreversibility in a stochastic process, as the ratio of the probabilities of a trajectory and its time reversal, cannot be accessed directly in experiment since the probability of a single trajectory is zero. We regularize this definition by considering, instead, the limiting ratio of probabilities for trajectories to remain in the tubular neighborhood of a smooth path and its time reversal. The resulting pathwise medium entropy production agrees with the formal expression from stochastic thermodynamics, and can be obtained from measurable tube probabilities. Estimating the latter from numerically sampled trajectories for Langevin dynamics yields excellent agreement with theory. By combining our measurement of pathwise entropy production with a Markov Chain Monte Carlo algorithm, we infer the entropy-production distribution for a transition path ensemble directly from short recorded trajectories. Our work enables the measurement of irreversibility along individual paths, and path ensembles, in a model-free manner.
△ Less
Submitted 13 April, 2022; v1 submitted 22 July, 2020;
originally announced July 2020.
-
Experimental measurement of relative path probabilities and stochastic actions
Authors:
Jannes Gladrow,
Ulrich F. Keyser,
R. Adhikari,
Julian Kappler
Abstract:
For diffusive stochastic dynamics, the probability to observe any individual trajectory is vanishingly small, making it unclear how to experimentally validate theoretical results for ratios of path probabilities. We provide the missing link between theory and experiment, by establishing a protocol to extract ratios of path probabilities from measured time series. For experiments on a single colloi…
▽ More
For diffusive stochastic dynamics, the probability to observe any individual trajectory is vanishingly small, making it unclear how to experimentally validate theoretical results for ratios of path probabilities. We provide the missing link between theory and experiment, by establishing a protocol to extract ratios of path probabilities from measured time series. For experiments on a single colloidal particle in a microchannel, we extract both ratios of path probabilities, and the most probable path for a barrier crossing, and find excellent agreement with independently calculated predictions based on the Onsager-Machlup stochastic action. Our experimental results at room temperature are found to be inconsistent with the low-noise Freidlin-Wentzell stochastic action, and we discuss under which circumstances the latter is expected to describe the most probable path. Furthermore, while the experimentally accessible ratio of path probabilities is uniquely determined, the formal path-integral action is known to depend on the time-discretization scheme used for deriving it; we reconcile these two seemingly contradictory facts by careful analysis of the time-slicing derivation of the path integral. Our experimental protocol enables us to probe probability distributions on path space, and allows us to relate theoretical single-trajectory results to measurement.
△ Less
Submitted 31 May, 2021; v1 submitted 30 June, 2020;
originally announced June 2020.
-
Stochastic action for tubes: Connecting path probabilities to measurement
Authors:
Julian Kappler,
Ronojoy Adhikari
Abstract:
The trajectories of diffusion processes are continuous but non-differentiable, and each occurs with vanishing probability. This introduces a gap between theory, where path probabilities are used in many contexts, and experiment, where only events with non-zero probability are measurable. Here we bridge this gap by considering the probability of diffusive trajectories to remain within a tube of sma…
▽ More
The trajectories of diffusion processes are continuous but non-differentiable, and each occurs with vanishing probability. This introduces a gap between theory, where path probabilities are used in many contexts, and experiment, where only events with non-zero probability are measurable. Here we bridge this gap by considering the probability of diffusive trajectories to remain within a tube of small but finite radius around a smooth path. This probability can be measured in experiment, via the rate at which trajectories exit the tube for the first time, thereby establishing a link between path probabilities and physical observables. Considering $N$-dimensional overdamped Langevin dynamics, we show that the tube probability can be obtained theoretically from the solution of the Fokker-Planck equation. Expressing the resulting exit rate as a functional of the path and ordering it as a power series in the tube radius, we identify the zeroth-order term as the Onsager-Machlup stochastic action, thereby elevating it from a mathematical construct to a physical observable. The higher-order terms reveal, for the first time, the form of the finite-radius contributions which account for fluctuations around the path. To demonstrate the experimental relevance of this action functional for tubes, we numerically sample trajectories of Brownian motion in a double-well potential, compute their exit rate, and show an excellent agreement with our analytical results. Our work shows that smooth tubes are surrogates for non-differentiable diffusive trajectories, and provide a direct way of comparing theoretical results on single trajectories, such as path-wise definitions of irreversibility, to measurement.
△ Less
Submitted 30 June, 2020;
originally announced June 2020.
-
Inference, prediction and optimization of non-pharmaceutical interventions using compartment models: the PyRoss library
Authors:
R. Adhikari,
Austen Bolitho,
Fernando Caballero,
Michael E. Cates,
Jakub Dolezal,
Timothy Ekeh,
Jules Guioth,
Robert L. Jack,
Julian Kappler,
Lukas Kikuchi,
Hideki Kobayashi,
Yuting I. Li,
Joseph D. Peterson,
Patrick Pietzonka,
Benjamin Remez,
Paul B. Rohrbach,
Rajesh Singh,
Günther Turk
Abstract:
PyRoss is an open-source Python library that offers an integrated platform for inference, prediction and optimisation of NPIs in age- and contact-structured epidemiological compartment models. This report outlines the rationale and functionality of the PyRoss library, with various illustrations and examples focusing on well-mixed, age-structured populations. The PyRoss library supports arbitrary s…
▽ More
PyRoss is an open-source Python library that offers an integrated platform for inference, prediction and optimisation of NPIs in age- and contact-structured epidemiological compartment models. This report outlines the rationale and functionality of the PyRoss library, with various illustrations and examples focusing on well-mixed, age-structured populations. The PyRoss library supports arbitrary structured models formulated stochastically (as master equations) or deterministically (as ODEs) and allows mid-run transitioning from one to the other. By supporting additional compartmental subdivision ad libitum, PyRoss can emulate time-since-infection models and allows medical stages such as hospitalization or quarantine to be modelled and forecast. The PyRoss library enables fitting to epidemiological data, as available, using Bayesian parameter inference, so that competing models can be weighed by their evidence. PyRoss allows fully Bayesian forecasts of the impact of idealized NPIs by convolving uncertainties arising from epidemiological data, model choice, parameters, and intrinsic stochasticity. Algorithms to optimize time-dependent NPI scenarios against user-defined cost functions are included. PyRoss's current age-structured compartment framework for well-mixed populations will in future reports be extended to include compartments structured by location, occupation, use of travel networks and other attributes relevant to assessing disease spread and the impact of NPIs. We argue that such compartment models, by allowing social data of arbitrary granularity to be combined with Bayesian parameter estimation for poorly-known disease variables, could enable more powerful and robust prediction than other approaches to detailed epidemic modelling. We invite others to use the PyRoss library for research to address today's COVID-19 crisis, and to plan for future pandemics.
△ Less
Submitted 19 May, 2020;
originally announced May 2020.
-
Molecular interpretation of the non-Newtonian viscoelastic behavior of liquid water at high frequencies
Authors:
Julius C. F. Schulz,
Alexander Schlaich,
Matthias Heyden,
Roland R. Netz,
Julian Kappler
Abstract:
Using classical as well as ab-initio molecular dynamics simulations, we calculate the frequency-dependent shear viscosity of pure water and water-glycerol mixtures. In agreement with recent experiments, we find deviations from Newtonian-fluid behavior in the THz regime. Based on an extension of the Maxwell model, we introduce a viscoelastic model to describe the observed viscosity spectrum of pure…
▽ More
Using classical as well as ab-initio molecular dynamics simulations, we calculate the frequency-dependent shear viscosity of pure water and water-glycerol mixtures. In agreement with recent experiments, we find deviations from Newtonian-fluid behavior in the THz regime. Based on an extension of the Maxwell model, we introduce a viscoelastic model to describe the observed viscosity spectrum of pure water. We find four relaxation modes in the spectrum which we attribute to i) hydrogen-bond network topology changes, ii) hydrogen-bond stretch vibrations of water pairs, iii) collective vibrations of water molecule triplets, and iv) librational excitations of individual water molecules. Our model quantitatively describes the viscoelastic response of liquid water on short timescales, where the hydrodynamic description via a Newtonian-fluid model breaks down.
△ Less
Submitted 18 March, 2020;
originally announced March 2020.
-
Non-Markovian barrier crossing with two-time-scale memory is dominated by the faster memory component
Authors:
Julian Kappler,
Victor B. Hinrichsen,
Roland R. Netz
Abstract:
We investigate non-Markovian barrier-crossing kinetics of a massive particle in one dimension in the presence of a memory function that is the sum of two exponentials with different memory times $τ_1$ and $τ_2$. Our Langevin simulations for the special case where both exponentials contribute equally to the total friction show that the barrier crossing time becomes independent of the longer memory…
▽ More
We investigate non-Markovian barrier-crossing kinetics of a massive particle in one dimension in the presence of a memory function that is the sum of two exponentials with different memory times $τ_1$ and $τ_2$. Our Langevin simulations for the special case where both exponentials contribute equally to the total friction show that the barrier crossing time becomes independent of the longer memory time if at least one of the two memory times is larger than the intrinsic diffusion time. When we associate memory effects with coupled degrees of freedom that are orthogonal to a one-dimensional reaction coordinate, this counterintuitive result shows that the faster orthogonal degrees of freedom dominate barrier-crossing kinetics in the non-Markovian limit and that the slower orthogonal degrees become negligible, quite contrary to the standard time-scale separation assumption and with important consequences for the proper setup of coarse-graining procedures in the non-Markovian case. By asymptotic matching and symmetry arguments, we construct a crossover formula for the barrier crossing time that is valid for general multi-exponential memory kernels. This formula can be used to estimate barrier-crossing times for general memory functions for high friction, i.e. in the overdamped regime, as well as for low friction, i.e. in the inertial regime. Typical examples where our results are important include protein folding in the high-friction limit and chemical reactions such as proton-transfer reactions in the low-friction limit.
△ Less
Submitted 30 July, 2019;
originally announced July 2019.
-
Extracting 3D Vascular Structures from Microscopy Images using Convolutional Recurrent Networks
Authors:
Russell Bates,
Benjamin Irving,
Bostjan Markelc,
Jakob Kaeppler,
Ruth Muschel,
Vicente Grau,
Julia A. Schnabel
Abstract:
Vasculature is known to be of key biological significance, especially in the study of cancer. As such, considerable effort has been focused on the automated measurement and analysis of vasculature in medical and pre-clinical images. In tumors in particular, the vascular networks may be extremely irregular and the appearance of the individual vessels may not conform to classical descriptions of vas…
▽ More
Vasculature is known to be of key biological significance, especially in the study of cancer. As such, considerable effort has been focused on the automated measurement and analysis of vasculature in medical and pre-clinical images. In tumors in particular, the vascular networks may be extremely irregular and the appearance of the individual vessels may not conform to classical descriptions of vascular appearance. Typically, vessels are extracted by either a segmentation and thinning pipeline, or by direct tracking. Neither of these methods are well suited to microscopy images of tumor vasculature. In order to address this we propose a method to directly extract a medial representation of the vessels using Convolutional Neural Networks. We then show that these two-dimensional centerlines can be meaningfully extended into 3D in anisotropic and complex microscopy images using the recently popularized Convolutional Long Short-Term Memory units (ConvLSTM). We demonstrate the effectiveness of this hybrid convolutional-recurrent architecture over both 2D and 3D convolutional comparators.
△ Less
Submitted 26 May, 2017;
originally announced May 2017.
-
Nonlinear fractional waves at elastic interfaces
Authors:
Julian Kappler,
Shamit Shrivastava,
Matthias F. Schneider,
Roland R. Netz
Abstract:
We derive the nonlinear fractional surface wave equation that governs compression waves at an interface that is coupled to a viscous bulk medium. The fractional character of the differential equation comes from the fact that the effective thickness of the bulk layer that is coupled to the interface is frequency dependent. The nonlinearity arises from the nonlinear dependence of the interface compr…
▽ More
We derive the nonlinear fractional surface wave equation that governs compression waves at an interface that is coupled to a viscous bulk medium. The fractional character of the differential equation comes from the fact that the effective thickness of the bulk layer that is coupled to the interface is frequency dependent. The nonlinearity arises from the nonlinear dependence of the interface compressibility on the local compression, which is obtained from experimental measurements and reflects a phase transition at the interface. Numerical solutions of our nonlinear fractional theory reproduce several experimental key features of surface waves in phospholipid monolayers at the air-water interface without freely adjustable fitting parameters. In particular, the propagation length of the surface wave abruptly increases at a threshold excitation amplitude. The wave velocity is found to be of the order of 40 cm/s both in experiments and theory and slightly increases as a function of the excitation amplitude. Nonlinear acoustic switching effects in membranes are thus shown to arise purely based on intrinsic membrane properties, namely the presence of compressibility nonlinearities that accompany phase transitions at the interface.
△ Less
Submitted 28 February, 2017;
originally announced February 2017.
-
Comment on "Correlation between Bulk Thermodynamic Measurements and the Low-Temperature-Resistance Plateau in SmB6"
Authors:
Kai Chen,
Jean-Paul Kappler
Abstract:
Low-temperature-resistivity plateau observed in $\rm SmB_6$ single crystal,which is due to surface, not bulk, conduction has been confirmed from electrical transport measurements. Recently, the correlation between bulk thermodynamic measurements and the low-temperature-resistance plateau in $\rm SmB_6$ have been investigated and a change in Sm valence at the surface has been obtained from x-ray ab…
▽ More
Low-temperature-resistivity plateau observed in $\rm SmB_6$ single crystal,which is due to surface, not bulk, conduction has been confirmed from electrical transport measurements. Recently, the correlation between bulk thermodynamic measurements and the low-temperature-resistance plateau in $\rm SmB_6$ have been investigated and a change in Sm valence at the surface has been obtained from x-ray absorption spectroscopy and x-ray magnetic circular dichroism. Here we show that the statement of the report are not supported by the results from x-ray absorption spectroscopy and x-ray magnetic circular dichroism.
△ Less
Submitted 12 March, 2015;
originally announced March 2015.
-
From itinerant to local transformation and critical point in Ni-rich Ce_2(Ni_{1-y}Pd_y)_2Sn
Authors:
J. G. Sereni,
G. Schmerber,
M. Gomez Berisso,
J. P. Kappler
Abstract:
Structural, magnetic (M) and thermal (C_m) studies on Ce_2(Ni_{1-y}Pd_y)_2Sn alloys are presented within the 0<y<0.55 range of concentration, showing evidences for itinerant to local electronic transformation. At variance with RKKY type interactions between localized moments μ_{eff}, the substitution of Ni by isoelectronic Pd leads the antiferromagnetic transition to decrease from T_N~3.8K to ~1.2…
▽ More
Structural, magnetic (M) and thermal (C_m) studies on Ce_2(Ni_{1-y}Pd_y)_2Sn alloys are presented within the 0<y<0.55 range of concentration, showing evidences for itinerant to local electronic transformation. At variance with RKKY type interactions between localized moments μ_{eff}, the substitution of Ni by isoelectronic Pd leads the antiferromagnetic transition to decrease from T_N~3.8K to ~1.2K between y=0 and 0.48, while M(H) measured at H=5Tesla and 1.8K rises from 0.12 up to 0.75μ_B/Ce-at. Furthermore, the C_m(T_N) jump increases with concentration whereas |θ_P| decreases. The magnetic entropy S_m(T) grows moderately with temperature for y=0 due to a significant contribution of excited levels at low energy, while at y=0.5 it shows a incipient plateau around S_m=Rln2. All these features reflect the progressive ground state transformation of from itinerant to a local character. Another peculiarity of this system is the nearly constant value of C_m(T_N) that ends in an "entropy bottleneck" as T_N decreases. Consequently, the system shows a critical point at y_{cr}~0.48 with signs of ferromagnetic behavior above H_{cr}~0.3T. A splitting of the C_m(T_N) maximum, tuned by field and concentration, indicates a competition between two magnetic phases, with respective peaks at T_N~1.2K and T_I~1.45K.
△ Less
Submitted 20 June, 2014;
originally announced June 2014.
-
Direct observation of a highly spin-polarized organic spinterface at room temperature
Authors:
F. Djeghloul,
F. Ibrahim,
M. Cantoni,
M. Bowen,
L. Joly,
S. Boukari,
P. Ohresser,
F. Bertran,
P. Lefèvre,
P. Thakur,
F. Scheurer,
T. Miyamachi,
R. Mattana,
P. Seneor,
A. Jaafar,
C. Rinaldi,
S. Javaid,
J. Arabski,
J. -P. Kappler,
W. Wulfhekel,
N. B. Brookes,
R. Bertacco,
A. Taleb-Ibrahimi,
M. Alouani,
E. Beaurepaire
, et al. (1 additional authors not shown)
Abstract:
The design of large-scale electronic circuits that are entirely spintronics-driven requires a current source that is highly spin-polarised at and beyond room temperature, cheap to build, efficient at the nanoscale and straightforward to integrate with semiconductors. Yet despite research within several subfields spanning nearly two decades, this key building block is still lacking. We experimental…
▽ More
The design of large-scale electronic circuits that are entirely spintronics-driven requires a current source that is highly spin-polarised at and beyond room temperature, cheap to build, efficient at the nanoscale and straightforward to integrate with semiconductors. Yet despite research within several subfields spanning nearly two decades, this key building block is still lacking. We experimentally and theoretically show how the interface between Co and phthalocyanine molecules constitutes a promising candidate. Spin-polarised direct and inverse photoemission experiments reveal a high degree of spin polarisation at room temperature at this interface. We measured a magnetic moment on the molecules's nitrogen pi orbitals, which substantiates an ab-initio theoretical description of highly spin-polarised charge conduction across the interface due to differing spinterface formation mechanims in each spin channel. We propose, through this example, a recipe to engineer simple organic-inorganic interfaces with remarkable spintronic properties that can endure well above room temperature.
△ Less
Submitted 2 October, 2012; v1 submitted 6 September, 2012;
originally announced September 2012.
-
Tri-critical point and suppression of the Shastry-Sutherland phase in Ce$_{2}$Pd$_{2}$Sn by Ni doping
Authors:
J. G. Sereni,
G. Schmerber,
M. Gómez Berisso,
B. Chevalier,
J. P. Kappler
Abstract:
Structural, magnetization and heat capacity measurements were performed on Ce$_2$(Pd$_{1-x}$Ni$_x$)$_2$Sn ($0 \leq x \leq 0.25$) alloys, covering the full range of the Mo$_2$FeB$_2$ structure stability. In this system, the two transitions observed in Ce$_2$Pd$_2$Sn (at $T_N=4.8$\,K and $T_C=2.1$\,K respectively) converge into a tri-critical point at $T_{cr}\approx 3.4$\,K for $x\approx 0.3$, where…
▽ More
Structural, magnetization and heat capacity measurements were performed on Ce$_2$(Pd$_{1-x}$Ni$_x$)$_2$Sn ($0 \leq x \leq 0.25$) alloys, covering the full range of the Mo$_2$FeB$_2$ structure stability. In this system, the two transitions observed in Ce$_2$Pd$_2$Sn (at $T_N=4.8$\,K and $T_C=2.1$\,K respectively) converge into a tri-critical point at $T_{cr}\approx 3.4$\,K for $x\approx 0.3$, where the intermediate antiferromagnetic AF phase is suppressed. The $T_N(x)$ phase boundary decrease is due to an incipient Kondo screening of the Ce-4f moments and local atomic disorder in the alloy. Both mechanisms affect the formation of Ce-magnetic dimers on which the Shastry-Sutherland lattice (SSL) builds up. On the contrary, the $T_C(x)$ transition to the ferromagnetic ground state increases as a consequence of the weakening of the AF-SSL phase. Applied magnetic field also suppresses the AF phase like in the stoichiometric compound.
△ Less
Submitted 3 October, 2011;
originally announced October 2011.
-
Structural, magnetic and thermal properties of the substitution series Ce$_2$(Pd$_{1-x}$Ni$_x$)$_2$Sn
Authors:
J. G. Sereni,
G. Schmerber,
A. Braghta,
B. Chevalier,
J. P. Kappler
Abstract:
Structural, magnetization and heat capacity studies were performed on Ce$_2$(Pd$_{1-x}$Ni$_x$)$_2$Sn ($0 \leq x \leq 1$) alloys. The substitution of Pd atoms by isoelectronic Ni leads to a change in the crystallographic structure from tetragonal (for $x \leq 0.3$) to centered orthorhombic lattice (for $x \geq 0.4$). The volume contraction thorough the series is comparable to the expected from the…
▽ More
Structural, magnetization and heat capacity studies were performed on Ce$_2$(Pd$_{1-x}$Ni$_x$)$_2$Sn ($0 \leq x \leq 1$) alloys. The substitution of Pd atoms by isoelectronic Ni leads to a change in the crystallographic structure from tetragonal (for $x \leq 0.3$) to centered orthorhombic lattice (for $x \geq 0.4$). The volume contraction thorough the series is comparable to the expected from the atomic size ratio between transition metal components. The consequent weak increase of the Kondo temperature drives the two transitions observed in Ce$_2$Pd$_2$Sn to merge at $x = 0.25$. After about a 1% of volume collapse at the structural modification, the system behaves as a weakly magnetic heavy fermion with an enhanced degenerate ground state. Notably, an incipient magnetic transition arises on the Ni-rich size. This unexpected behavior is discussed in terms of an enhancement of the density of states driven by the increase of the $4f$-conduction band hybridization and the incipient contribution of the first excited crystal field doublet on the ground state properties.
△ Less
Submitted 1 March, 2011;
originally announced March 2011.
-
Suppression of the Shastry-Sutherland phase in Ce$_{2}$Pd$_{2}$Sn at a field induced critical point
Authors:
J. G. Sereni,
M. Gomez Berisso,
G. Schmerber,
A. Braghta,
J. P. Kappler
Abstract:
The magnetic phase diagram of Ce$_2$Pd$_2$Sn is investigated through the field dependence of thermal, transport and magnetic measurements performed at low temperature. The upper transition, $T_M=4.8$ K is practically not affected by magnetic field up to B=1 T, whereas the lower one $T_C(B)$ rapidly increases from 2.1 K joining $T_M$ in a critical point at $T_{cr}=(4.2\pm 0.3$)K for…
▽ More
The magnetic phase diagram of Ce$_2$Pd$_2$Sn is investigated through the field dependence of thermal, transport and magnetic measurements performed at low temperature. The upper transition, $T_M=4.8$ K is practically not affected by magnetic field up to B=1 T, whereas the lower one $T_C(B)$ rapidly increases from 2.1 K joining $T_M$ in a critical point at $T_{cr}=(4.2\pm 0.3$)K for $B_{cr}=(0.12\pm 0.03)$ T. At that point the intermediate phase, previously described as an unstable Shastry-Sutherland phase, is suppressed. A detailed analysis around the critical point reveals a structure in the maximum of the $\partial M/\partial B(B)$ derivative which could be related to the formation of a novel phase in that critical region.
△ Less
Submitted 27 January, 2010; v1 submitted 7 January, 2010;
originally announced January 2010.
-
Temperature dependence of the pre-edge structure in the Ti K-edge x-ray absorption spectrum of rutile
Authors:
Olivier Durmeyer,
Eric Beaurepaire,
Jean-Paul Kappler,
Christian Brouder,
François Baudelet
Abstract:
The temperature dependence of the pre-edge features in x-ray absorption spectroscopy is reviewed. Then, the temperature dependence of the pre-edge structure at the K-edge of titanium in rutile TiO2 is measured at low and room temperature. The first two peaks grow with temperature. The fact that these two peaks also correspond to electric quadrupole transitions is explained by a recently proposed…
▽ More
The temperature dependence of the pre-edge features in x-ray absorption spectroscopy is reviewed. Then, the temperature dependence of the pre-edge structure at the K-edge of titanium in rutile TiO2 is measured at low and room temperature. The first two peaks grow with temperature. The fact that these two peaks also correspond to electric quadrupole transitions is explained by a recently proposed theory.
△ Less
Submitted 12 December, 2009;
originally announced December 2009.
-
Unstable Shastry-Sutherland phase in Ce$_{2}$Pd$_{2}$Sn
Authors:
J. G. Sereni,
M. Gomez-Berisso,
A. Braghta,
G. Schmerber,
J. P. Kappler
Abstract:
Thermal ($C_P$), magnetic ($M$ and $χ_{ac}$) and transport ($ρ$) measurements on Ce$_{2}$Pd$_{2}$Sn are reported. High temperature properties are well described by the presence of two excited crystal field levels at $(65\pm 5)$K and $(230\pm 20)$K, with negligible hybridization (Kondo) effects. According to literature, two transitions were observed at $T_M=4.8$K and $T_C =2.1$K respectively. The…
▽ More
Thermal ($C_P$), magnetic ($M$ and $χ_{ac}$) and transport ($ρ$) measurements on Ce$_{2}$Pd$_{2}$Sn are reported. High temperature properties are well described by the presence of two excited crystal field levels at $(65\pm 5)$K and $(230\pm 20)$K, with negligible hybridization (Kondo) effects. According to literature, two transitions were observed at $T_M=4.8$K and $T_C =2.1$K respectively. The upper transition cannot be considered as a standard anti-ferromagnetic because of frustration effects in a triangular network of Ce-atoms and the positive sign of the paramagnetic temperature $θ_P^{LT} = 4.4$K. The nature of the this intermediated phase is described accounting for the formation of ferromagnetic (F) Ce-dimers disposed in a quasi-2D square lattice, resembling a Shastry-Sutherland pattern. According to hysteretic features in $ρ(T)$ and $χ_{ac}(T)$, the lower F-transition is of first order, with $C_P(T<T_C)$ revealing a gap of anisotropy $E_g\approx 7$K.
△ Less
Submitted 12 March, 2009;
originally announced March 2009.
-
Paramagnetism of the Co sublattice in ferromagnetic Zn$_{1-x}$Co$_{x}$O films
Authors:
A. Barla,
G. Schmerber,
E. Beaurepaire,
A. Dinia,
H. Bieber,
S. Colis,
F. Scheurer,
J. -P. Kappler,
P. Imperia,
F. Nolting,
F. Wilhelm,
A. Rogalev,
D. Muller,
J. J. Grob
Abstract:
Using the spectroscopies based upon x-ray absorption, we have studied the structural and magnetic properties of Zn$_{1-x}$Co$_{x}$O films ($x$ = 0.1 and 0.25) produced by reactive magnetron sputtering. These films show ferromagnetism with a Curie temperature $T_{\mathrm{C}}$ above room temperature in bulk magnetization measurements. Our results show that the Co atoms are in a divalent state and…
▽ More
Using the spectroscopies based upon x-ray absorption, we have studied the structural and magnetic properties of Zn$_{1-x}$Co$_{x}$O films ($x$ = 0.1 and 0.25) produced by reactive magnetron sputtering. These films show ferromagnetism with a Curie temperature $T_{\mathrm{C}}$ above room temperature in bulk magnetization measurements. Our results show that the Co atoms are in a divalent state and in tetrahedral coordination, thus substituting Zn in the wurtzite-type structure of ZnO. However, x-ray magnetic circular dichroism at the Co \textit{L}$_{2,3}$ edges reveals that the Co 3\textit{d} sublattice is paramagnetic at all temperatures down to 2 K, both at the surface and in the bulk of the films. The Co 3\textit{d} magnetic moment at room temperature is considerably smaller than that inferred from bulk magnetisation measurements, suggesting that the Co 3\textit{d} electrons are not directly at the origin of the observed ferromagnetism.
△ Less
Submitted 7 December, 2006;
originally announced December 2006.
-
Probing the gamma-alpha Transition in Bulk Ce under Pressure: A Direct Investigation by Resonant Inelastic X-ray Scattering
Authors:
Jean-Pascal Rueff,
Jean-Paul Itié,
M. Taguchi,
C. F. Hague,
J. -M. Mariot,
R. Delaunay,
J. -P. Kappler,
N. Jaouen
Abstract:
We report on the most complete investigation to date of the 4f-electron properties at the gamma-alpha transition in elemental Ce by resonant inelastic x-ray scattering (RIXS). The Ce 2p3d-RIXS spectra were measured directly in the bulk material as a function of pressure through the transition. The spectra were simulated within the Anderson impurity model. The occupation number nf was derived fro…
▽ More
We report on the most complete investigation to date of the 4f-electron properties at the gamma-alpha transition in elemental Ce by resonant inelastic x-ray scattering (RIXS). The Ce 2p3d-RIXS spectra were measured directly in the bulk material as a function of pressure through the transition. The spectra were simulated within the Anderson impurity model. The occupation number nf was derived from the calculations in both gamma- and alpha-phases in the ground state along with the f doubleoccupancy. We find that the electronic structure changes result mainly from band formation of 4f electrons which concurs with reduced electron correlation and increased Kondo screening at high pressure.
△ Less
Submitted 7 February, 2006;
originally announced February 2006.
