Aging to Equilibrium Dynamics of SiO2

Molecular dynamics computer simulations are used to study the aging dynamics of SiO2 (modeled by the BKS model). Starting from fully equilibrated configurations at high temperatures T_i =5000K/3760K the system is quenched to lower temperatures T_f=2500K, 2750K, 3000K, 3250K and observed after a waiting time t_w. Since the simulation runs are long enough to reach equilibrium at T_f, we are able to study the transition from out-of-equilibrium to equilibrium dynamics. We present results for the partial structure factors, for the generalized incoherent intermediate scattering function C_q(t_w, t_w+t), and for the mean square displacement msd(t_w,t_w+t). We conclude that there are three different t_w regions: (I) At very short waiting times, C_q(t_w, t_w+t) decays very fast without forming a plateau. Similarly msd(t_w,t_w+t) increases without forming a plateau. (II) With increasing t_w a plateau develops in C_q(t_w, t_w+t) and msd(t_w,t_w+t). For intermediate waiting times the plateau height is independent of t_w and T_i. Time superposition applies, i.e. C_q=C_q(t/t_r) where t_r=t_r(t_w) is a waiting time dependent decay time. Furthermore C_q=C(q,t_w,t_w+t) scales as C_q=C(q,z(t_w,t) where z is a function of t_w and t only, i.e. independent of q. (III) At large t_w the system reaches equilibrium, i.e. C_q(t_w,t_w+t) and msd(t_w,t_w+t) are independent of t_w and T_i. For C_q(t_w,t_w+t) we find that the time superposition of intermediate waiting times (II) includes the equilibrium curve (III).Comment: 9 pages, 11 figures, submission to PR

Some new results concerning the vacuum in Dirac Hole Theory

In Dirac's hole theory the vacuum state is generally believed to be the state of minimum energy. It will be shown that this is not, in fact, the case and that there must exist states in hole theory with less energy than the vacuum state. It will be shown that energy can be extracted from the hole theory vacuum state through the application of an electric field.Comment: Accepted by Physica Scripta, 19 page

The global electroweak fit at NNLO and prospects for the LHC and ILC

For a long time, global fits of the electroweak sector of the Standard Model (SM) have been used to exploit measurements of electroweak precision observables at lepton colliders (LEP, SLC), together with measurements at hadron colliders (Tevatron, LHC), and accurate theoretical predictions at multi-loop level, to constrain free parameters of the SM, such as the Higgs and top masses. Today, all fundamental SM parameters entering these fits are experimentally determined, including information on the Higgs couplings, and the global fits are used as powerful tools to assess the validity of the theory and to constrain scenarios for new physics. Future measurements at the Large Hadron Collider (LHC) and the International Linear Collider (ILC) promise to improve the experimental precision of key observables used in the fits. This paper presents updated electroweak fit results using newest NNLO theoretical predictions, and prospects for the LHC and ILC. The impact of experimental and theoretical uncertainties is analysed in detail. We compare constraints from the electroweak fit on the Higgs couplings with direct LHC measurements, and examine present and future prospects of these constraints using a model with modified couplings of the Higgs boson to fermions and bosons.Comment: 26 pages, 9 figure

Symmetry adapted ro-vibrational basis functions for variational nuclear motion calculations: TROVE approach

We present a general, numerically motivated approach to the construction of symmetry adapted basis functions for solving ro-vibrational Schr\"{o}dinger equations. The approach is based on the property of the Hamiltonian operator to commute with the complete set of symmetry operators and hence to reflect the symmetry of the system. The symmetry adapted ro-vibrational basis set is constructed numerically by solving a set of reduced vibrational eigenvalue problems. In order to assign the irreducible representations associated with these eigenfunctions, their symmetry properties are probed on a grid of molecular geometries with the corresponding symmetry operations. The transformation matrices are re-constructed by solving over-determined systems of linear equations related to the transformation properties of the corresponding wavefunctions on the grid. Our method is implemented in the variational approach TROVE and has been successfully applied to a number of problems covering the most important molecular symmetry groups. Several examples are used to illustrate the procedure, which can be easily applied to different types of coordinates, basis sets, and molecular systems

Dust emission at 8-mic and 24-mic as Diagnostics of HII Region Radiative Transfer

We use the Spitzer SAGE survey of the Magellanic Clouds to evaluate the relationship between the 8-mic PAH emission, 24-mic hot dust emission, and HII region radiative transfer. We confirm that in the higher-metallicity Large Magellanic Cloud, PAH destruction is sensitive to optically thin conditions in the nebular Lyman continuum: objects identified as optically thin candidates based on nebular ionization structure show 6 times lower median 8-mic surface brightness (0.18 mJy arcsec^-2) than their optically thick counterparts (1.2 mJy arcsec^-2). The 24-mic surface brightness also shows a factor of 3 offset between the two classes of objects (0.13 vs 0.44 mJy arcsec^-2, respectively), which is driven by the association between the very small dust grains and higher density gas found at higher nebular optical depths. In contrast, PAH and dust formation in the low-metallicity Small Magellanic Cloud is strongly inhibited such that we find no variation in either 8-mic or 24-mic emission between our optically thick and thin samples. This is attributable to extremely low PAH and dust production together with high, corrosive UV photon fluxes in this low-metallicity environment. The dust mass surface densities and gas-to-dust ratios determined from dust maps using Herschel HERITAGE survey data support this interpretation.Comment: Accepted to ApJ, May 15, 2017. 10 pages, 9 figure

Dynamical mechanism of atrial fibrillation: a topological approach

While spiral wave breakup has been implicated in the emergence of atrial fibrillation, its role in maintaining this complex type of cardiac arrhythmia is less clear. We used the Karma model of cardiac excitation to investigate the dynamical mechanisms that sustain atrial fibrillation once it has been established. The results of our numerical study show that spatiotemporally chaotic dynamics in this regime can be described as a dynamical equilibrium between topologically distinct types of transitions that increase or decrease the number of wavelets, in general agreement with the multiple wavelets hypothesis. Surprisingly, we found that the process of continuous excitation waves breaking up into discontinuous pieces plays no role whatsoever in maintaining spatiotemporal complexity. Instead this complexity is maintained as a dynamical balance between wave coalescence -- a unique, previously unidentified, topological process that increases the number of wavelets -- and wave collapse -- a different topological process that decreases their number.Comment: 15 pages, 14 figure