BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
S. Aubry contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We propose a new and general formalism for elementary chemical reactions where quantum electronic variables are used as reaction coordinates. This formalism is in principle applicable to all kinds of chemical reactions ionic or covalent. Our theory reveals the existence of an intermediate situation between ionic and covalent which may be almost barrierless and isoenegetic and which should be of high interest for understanding biochemistry.
Recent molecular dynamics simulation methods have enabled thermal conductivity of bulk materials to be estimated. In these simulations, periodic boundary conditions are used to extend the system dimensions to the thermodynamic limit. Such a strategy cannot be used for nanostructures with finite dimensions which are typically much larger than it is possible to simulate directly. To bridge the length scales between the simulated and the actual nanostructures, we perform large-scale molecular dynamics calculations of thermal conductivities at different system dimensions to examine a recently developed conductivity vs. dimension scaling theory for both film and wire configurations. We demonstrate that by an appropriate application of the scaling law, reliable interpolations can be used to accurately predict thermal conductivity of films and wires as a function of film thickness or wire radius at realistic length scales from molecular dynamics simulations. We apply this method to predict thermal conductivities for GaN wurtzite nanostructures.
Significant differences exist among literature for thermal conductivity of various systems computed using molecular dynamics simulation. In some cases, unphysical results, for example, negative thermal conductivity, have been found. Using GaN as an example case and the direct non-equilibrium method, extensive molecular dynamics simulations and Monte Carlo analysis of the results have been carried out to quantify the uncertainty level of the molecular dynamics methods and to identify the conditions that can yield sufficiently accurate calculations of thermal conductivity. We found that the errors of the calculations are mainly due to the statistical thermal fluctuations. Extrapolating results to the limit of an infinite-size system tend to magnify the errors and occasionally lead to unphysical results. The error in bulk estimates can be reduced by performing longer time averages using properly selected systems over a range of sample lengths. If the errors in the conductivity estimates associated with each of the sample lengths are kept below a certain threshold, the likelihood of obtaining unphysical bulk values becomes insignificant. Using a Monte-Carlo approach developed here, we have determined the probability distributions for the bulk thermal conductivities obtained using the direct method. We also have observed a nonlinear effect that can become a source of significant errors. For the extremely accurate results presented here, we predict a [0001] GaN thermal conductivity of 185 $\rm{W/K \cdot m}$ at 300 K, 102 $\rm{W/K \cdot m}$ at 500 K, and 74 $\rm{W/K \cdot m}$ at 800 K. Using the insights obtained in the work, we have achieved a corresponding error level (standard deviation) for the bulk (infinite sample length) GaN thermal conductivity of less than 10 $\rm{W/K \cdot m}$, 5 $\rm{W/K \cdot m}$, and 15 $\rm{W/K \cdot m}$ at 300 K, 500 K, and 800 K respectively.
We investigate the long time behavior of a wavepacket initially localized at a single site $n_0$ in translationally invariant harmonic and anharmonic chains with random interactions. In the harmonic case, the energy profile $ \bar{< e_n(t)>}$ averaged on time and disorder decays for large $|n-n_0|$ as a power law $\bar{< e_n(t)>}\approx C|n-n_0|^{-η}$ where $η=5/2$ and 3/2 for initial displacement and momentum excitations, respectively. The prefactor $C$ depends on the probability distribution of the harmonic coupling constants and diverges in the limit of weak disorder. As a consequence, the moments $< m_ν(t)>$ of the energy distribution averaged with respect to disorder diverge in time as $t^{β(ν)}$ for $ν\geq 2$, where $β=ν+1-η$ for $ν>η-1$. Molecular dynamics simulations yield good agreement with these theoretical predictions. Therefore, in this system, the second moment of the wavepacket diverges as a function of time despite the wavepacket is not spreading. Thus, this only criteria often considered earlier as proving the spreading of a wave packet, cannot be considered as sufficient in any model. The anharmonic case is investigated numerically. It is found for intermediate disorder, that the tail of the energy profile becomes very close to those of the harmonic case. For weak and strong disorder, our results suggest that the crossover to the harmonic behavior occurs at much larger $|n-n_0|$ and larger time.
We investigate the dynamics of a macroscopic system which consists of an anharmonic subsystem embedded in an arbitrary harmonic lattice, including quenched disorder. Elimination of the harmonic degrees of freedom leads to a nonlinear Langevin equation for the anharmonic coordinates. For zero temperature, we prove that the support of the Fourier transform of the memory kernel and of the time averaged velocity-velocity correlations functions of the anharmonic system can not overlap. As a consequence, the asymptotic solutions can be constant, periodic,quasiperiodic or almost periodic, and possibly weakly chaotic. For a sinusoidal trajectory with frequency $Ω$ we find that the energy $E_T$ transferred to the harmonic system up to time $T$ is proportional to $T^α$. If $Ω$ equals one of the phonon frequencies $ω_ν$, it is $α=2$. We prove that there is a full measure set such that for $Ω$ in this set it is $α=0$, i.e. there is no energy dissipation. Under certain conditions there exists a zero measure set such that for $Ω\in this set the dissipation rate is nonzero and may be subdissipative $(0 \leq α< 1)$ or superdissipative $(1 <α\leq 2)$. Consequently, the harmonic bath does act as an anomalous thermostat. Intraband discrete breathers are such solutions which do not relax. We prove for arbitrary anharmonicity and small but finite coupling that intraband discrete breathers with frequency $Ω$ exist for all $Ω$ in a Cantor set $\mathcal{C}(k)$ of finite Lebesgue measure. This is achieved by estimating the contribution of small denominators appearing in the memory kernel. For $Ω\in\mathcal{C}(k)$ the small denominators do not lead to divergencies such that this kernel is a smooth and bounded function in $t$.
We present a practical algorithm for partially relaxing multiwell energy densities such as pertain to materials undergoing martensitic phase transitions. The algorithm is based on sequential lamination, but the evolution of the microstructure during a deformation process is required to satisfy a continuity constraint, in the sense that the new microstructure should be reachable from the preceding one by a combination of branching and pruning operations. All microstructures generated by the algorithm are in static and configurational equilibrium. Owing to the continuity constrained imposed upon the microstructural evolution, the predicted material behavior may be path-dependent and exhibit hysteresis. In cases in which there is a strict separation of micro and macrostructural lengthscales, the proposed relaxation algorithm may effectively be integrated into macroscopic finite-element calculations at the subgrid level. We demonstrate this aspect of the algorithm by means of a numerical example concerned with the indentation of an Cu-Al-Ni shape memory alloy by a spherical indenter.
We study the dynamics of the discrete nonlinear Schr{ö}dinger lattice initialized such that a very long transitory period of time in which standard Boltzmann statistics is insufficient is reached. Our study of the nonlinear system locked in this {\em non-Gibbsian} state focuses on the dynamics of discrete breathers (also called intrinsic localized modes). It is found that part of the energy spontaneously condenses into several discrete breathers. Although these discrete breathers are extremely long lived, their total number is found to decrease as the evolution progresses. Even though the total number of discrete breathers decreases we report the surprising observation that the energy content in the discrete breather population increases. We interpret these observations in the perspective of discrete breather creation and annihilation and find that the death of a discrete breather cause effective energy transfer to a spatially nearby discrete breather. It is found that the concepts of a multi-frequency discrete breather and of internal modes is crucial for this process. Finally, we find that the existence of a discrete breather tends to soften the lattice in its immediate neighborhood, resulting in high amplitude thermal fluctuation close to an existing discrete breather. This in turn nucleates discrete breather creation close to a already existing discrete breather.