Time-lapse monitoring of climate effects on earthworks using surface waves

The UK’s transportation network is supported by critical geotechnical assets (cuttings/embankments/dams) that require sustainable, cost-effective management, while maintaining an appropriate service level to meet social, economic, and environmental needs. Recent effects of extreme weather on these geotechnical assets have highlighted their vulnerability to climate variations. We have assessed the potential of surface wave data to portray the climate-related variations in mechanical properties of a clay-filled railway embankment. Seismic data were acquired bimonthly from July 2013 to November 2014 along the crest of a heritage railway embankment in southwest England. For each acquisition, the collected data were first processed to obtain a set of Rayleigh-wave dispersion and attenuation curves, referenced to the same spatial locations. These data were then analyzed to identify a coherent trend in their spatial and temporal variability. The relevance of the observed temporal variations was also verified with respect to the experimental data uncertainties. Finally, the surface wave dispersion data sets were inverted to reconstruct a time-lapse model of S-wave velocity for the embankment structure, using a least-squares laterally constrained inversion scheme. A key point of the inversion process was constituted by the estimation of a suitable initial model and the selection of adequate levels of spatial regularization. The initial model and the strength of spatial smoothing were then kept constant throughout the processing of all available data sets to ensure homogeneity of the procedure and comparability among the obtained VS sections. A continuous and coherent temporal pattern of surface wave data, and consequently of the reconstructed VS models, was identified. This pattern is related to the seasonal distribution of precipitation and soil water content measured on site.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Path-integral study of a two-dimensional Lennard-Jones glass

The glass transition in a quantum Lennard-Jones mixture is investigated by constant-volume path-integral simulations. Particles are assumed to be distinguishable, and the strength of quantum effects is varied by changing h from zero (the classical case) to one (corresponding to a highly quantum-mechanical regime). Quantum delocalization and zero point energy drastically reduce the sensitivity of structural and thermodynamic properties to the glass transition. Nevertheless, the glass transition temperature T-g can be determined by analyzing the phase space mobility of path-integral centroids. At constant volume, the T-g of the simulated model increases monotonically with increasing h. Low temperature tunneling centers are identified, and the quantum versus thermal character of each center is analyzed. The relation between these centers and soft quasilocalized harmonic vibrations is investigated. Periodic minimizations of the potential energy with respect to the positions of the particles are performed to determine the inherent structure of classical and quantum glassy samples. The geometries corresponding to these energy minima are found to be qualitatively similar in all cases. Systematic comparisons for ordered and disordered structures, harmonic and anharmonic dynamics, classical and quantum systems show that disorder, anharmonicity, and quantum effects are closely interlinked.

Equilibrium polymerization of cyclic carbonate oligomers. II. Role of multiple active sites

Ring opening polymerization of bisphenol A polycarbonate is studied by Monte Carlo simulations of a model comprising a fixed number of Lennard-Jones particles and harmonic bonds [J. Chem. Phys. 115, 3895 (2001)]. Bond interchanges produced by a low concentration (0.10%less than or equal toc(a)less than or equal to0.36%) of chemically active particles lead to equilibrium polymerization. There is a continuous transition in both 2D and 3D from unpolymerized cyclic oligomers at low density to a system of linear chains at high density, and the polymeric phase is much more stable in three dimensions than in two. The steepness of the polymerization transition increases rapidly as c(a) decreases, suggesting that it is discontinuous in the limit c(a)-->0. The transition is entropy driven, since the average potential energy increases systematically upon polymerization, and there is a steady decline in the degree of polymerization as the temperature is lowered. The mass distribution functions for open chains and for rings are unimodal, with exponentially decaying tails that can be fitted by Zimm-Schulz functions and simpler exponential forms. (C) 2002 American Institute of Physics.

Equilibrium polymerization of cyclic carbonate oligomers

A model of the polymerization of ring oligomers of bisphenol A polycarbonate (BPA-PC) is used to investigate the influence of dimensionality (2D or 3D), density and temperature on the size distribution of the polymer chains. The polymerization step is catalyzed by a single active particle, conserves the number and type of the chemical bonds, and occurs without a significant gain in either potential energy or configurational entropy. Monte Carlo and molecular dynamics simulations show that polymerization of cyclic oligomers occurs readily at high density and is driven by the entropy associated with the distribution of interparticle bonds. Polymerization competes at lower densities with long range diffusion, which favors small molecular species, and is prevented if the system is sufficiently dilute. Polymerization occurs in 2D via a weakly first order transition as a function of density and is characterized by low hysteresis and large fluctuations in the size of polymer chains. Polymerization occurs more readily in 3D than in 2D, and is favored by increasing temperature, as expected for an entropy-driven process. (C) 2001 American Institute of Physics.

Simple models of complex aggregation: Vesicle formation by soft repulsive spheres with dipolelike interactions

Structural and thermodynamic properties of spherical particles carrying classical spins are investigated by Monte Carlo simulations. The potential energy is the sum of short range, purely repulsive pair contributions, and spin-spin interactions. These last are of the dipole-dipole form, with however, a crucial change of sign. At low density and high temperature the system is a homogeneous fluid of weakly interacting particles and short range spin correlations. With decreasing temperature particles condense into an equilibrium population of free floating vesicles. The comparison with the electrostatic case, giving rise to predominantly one-dimensional aggregates under similar conditions, is discussed. In both cases condensation is a continuous transformation, provided the isotropic part of the interatomic potential is purely repulsive. At low temperature the model allows us to investigate thermal and mechanical properties of membranes. At intermediate temperatures it provides a simple model to investigate equilibrium polymerization in a system giving rise to predominantly two-dimensional aggregates.

Charge exchange and dissociative processes in collisions of slow He2+ ions with H2O molecules

Experimental and theoretical studies of one-electron capture in collisions of He2+ ions with H2O molecules have been carried out in the range 0.025-12 keV amu(-1) corresponding to typical solar wind velocities of 70-1523 km s(-1). Translational energy spectroscopy (TES), photon emission spectroscopy (PES), and fragment ion spectroscopy were employed to identify and quantify the collision mechanisms involved. Cross sections for selective single electron capture into n=1, 2, and 3 states of the He+ ion were obtained using TES while PES provided cross sections for capture into the He+(2p) and He+(3p) states. Our model calculations show that He+(n=2) and He+(n=3) formation proceeds via a single-electron process governed by the nucleus-electron interaction. In contrast, the He+(1s) formation mechanism involves an exothermic two-electron process driven by the electron-electron interaction, where the potential energy released by the electron capture is used to remove a second electron thereby resulting in fragmentation of the H2O molecule. This process is found to become increasingly important as the collision energy decreases. The experimental cross sections are found to be in reasonable agreement with cross sections calculated using the Demkov and Landau-Zener models.

Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.

Absolute single and multiple charge exchange cross sections for highly charged C, O, and Ne ions on H2O, CO, and CO2

Reported herein are measured absolute single, double, and triple charge exchange (CE) cross sections for the highly charged ions (HCIs) Cq+ (q=5,6), Oq+ (q=6,7,8), and Neq+ (q=7,8) colliding with the molecular species H2O, CO, and CO2. Present data can be applied to interpreting observations of x-ray emissions from comets as they interact with the solar wind. As such, the ion impact energies of 7.0q keV (1.62–3.06 keV/amu) are representative of the fast solar wind, and data at 1.5q keV for O6+ (0.56 keV/amu) on CO and CO2 and 3.5q keV for O5+ (1.09 keV/amu) on CO provide checks of the energy dependence of the cross sections at intermediate and typical slow solar wind velocities. The HCIs are generated within a 14 GHz electron cyclotron resonance ion source. Absolute CE measurements are made using a retarding potential energy analyzer, with measurement of the target gas cell pressure and incident and final ion currents. Trends in the cross sections are discussed in light of the classical overbarrier model (OBM), extended OBM, and with recent results of the classical trajectory Monte Carlo theory.

Dust lattice wave dispersion relations in two-dimensional hexagonal crystals including the effect of dust charge polarization

A dusty plasma crystalline configuration with equal charge dust grains and mass is considered. Both charge and mass of each dust species are taken to be constant. Two differential equations for a two-dimensional hexagonal crystal on the basis of a Yukawa-type potential energy and a

Nonlinearly modulated dust lattice wave packets in two-dimensional hexagonal dust crystals

The amplitude modulation of dust lattice waves (DLWs) propagating in a two-dimensional hexagonal dust crystal is investigated in a continuum approximation, accounting for the effect of dust charge polarization (dressed interactions). A dusty plasma crystalline configuration with constant dust grain charge and mass is considered. The dispersion relation and the group velocity for DLWs are determined for wave propagation in both longitudinal and transverse directions. The reductive perturbation method is used to derive a (2+1)-dimensional nonlinear Schrodinger equation (NLSE). New expressions for the coefficients of the NLSE are derived and compared, for a Yukawa-type potential energy and for a

Intelligent demand response electric water heating to integrate wind power in Ireland

Dwindling fossil fuel resources and pressures to reduce greenhouse gas (GHG) emissions will result in a more diverse range of generation portfolios for future electricity systems. Irrespective of the portfolio mix the overarching requirement for all electricity suppliers and system operators is that supply instantaneously meets demand and that robust operating standards are maintained to ensure a consistent supply of high quality electricity to end-users. Therefore all electricity market participants will ultimately need to use a variety of tools to balance the power system. Thus the role of demand side management (DSM) with energy storage will be paramount to integrate future diverse generation portfolios. Electric water heating (EWH) has been studied previously, particularly at the domestic level to provide load control, peak shave and to benefit end-users financially with lower bills, particularly in vertically integrated monopolies. In this paper, a continuous Direct Load Control (DLC) EWH algorithm is applied in a liberalized market environment using actual historical electricity system and market data to examine the potential energy savings, cost reductions and electricity system operational improvements.

Measurements of fast electron scaling generated by petawatt laser systems

Fast electron energy spectra have been measured for a range of intensities between 10(18) and 10(21) W cm(-2) and for different target materials using electron spectrometers. Several experimental campaigns were conducted on petawatt laser facilities at the Rutherford Appleton Laboratory and Osaka University, where the pulse duration was varied from 0.5 to 5 ps relevant to upcoming fast ignition integral experiments. The incident angle was also changed from normal incidence to 40 degrees in p-polarized. The results confirm a reduction from the ponderomotive potential energy on fast electrons at the higher intensities under the wide range of different irradiation conditions.

Structure and interactions of ultracold Yb ions and Rb atoms

In order to study ultracold charge-transfer processes in hybrid atom-ion traps, we have mapped out the potential-energy curves and molecular parameters for several low-lying states of the Rb, Yb+ system. We employ both a multireference configuration interaction and a full configuration interaction (FCI) approach. Turning points, crossing points, potential minima, and spectroscopic molecular constants are obtained for the lowest five molecular states. Long-range parameters, including the dispersion coefficients, are estimated from our ab initio data. The separated-atom ionization potentials and atomic polarizability of the ytterbium atom (ad=128.4 atomic units) are in good agreement with experiment and previous calculations. We present some dynamical calculations for (adiabatic) scattering lengths for the two lowest (Yb, Rb+) channels that were carried out in our work. However, we find that the pseudopotential approximation is rather limited in validity and only applies to nK temperatures. The adiabatic scattering lengths for both the triplet and singlet channels indicate that both are large and negative in the FCI approximation.

Buckling mode transition in hat-stiffened composite panels loaded in uniaxial compression

A combined experimental and analytical study of a hat-stiffened carbon-fibre composite panel loaded in uniaxial compression was investigated. A buckling mode transition was observed in the panel's skin bay which was not captured using non-linear finite-element analysis. Good correlation between experimental and numerical strain and displacement results was achieved in the prebuckling and initial postbuckling region of the loading history. A Marguerre-type Rayleigh-Ritz energy method was applied to the skin bay using representative displacement functions of permissible mode shapes to explain the mode transition phenomenon. The central criterion of this method was based on the assumption that a change in mode shape occurred such that the total potential energy of the structure was maintained at a minimum. The ultimate strength of the panel was limited by the column buckling strength of the hat-stiffeners.