Quasi-classical model of non-destructive wavepacket manipulation by intense ultrashort nonresonant laser pulses

A quasi-classical model (QCM) of nuclear wavepacket generation, modification and imaging by three intense ultrafast near-infrared laser pulses has been developed. Intensities in excess of 10(13) W cm(-2) are studied, the laser radiation is non-resonant and pulse durations are in the few-cycle regime, hence significantly removed from the conditions typical of coherent control and femtochemistry. The 1s sigma ground state of the D-2 precursor is projected onto the available electronic states in D-2(+) (1s sigma(g) ground and 2p sigma(u) dissociative) and D+ + D+ (Coulomb explosion) by tunnel ionization by an ultrashort 'pump' pulse, and relative populations are found numerically. A generalized non-adiabatic treatment allows the dependence of the initial vibrational population distribution on laser intensity to be calculated. The wavepacket is approximated as a classical ensemble of particles moving on the 1s sigma(g) potential energy surface (PES), and hence follow trajectories of different amplitudes and frequencies depending on the initial vibrational state. The 'control' pulse introduces a time-dependent polarization of the molecular orbital, causing the PES to be modified according to the dynamic Stark effect and the transition dipole. The trajectories adjust in amplitude, frequency and phase-offset as work is done on or by the resulting force; comparing the perturbed and unperturbed trajectories allows the final vibrational state populations and phases to be determined. The action of the 'probe' pulse is represented by a discrete internuclear boundary, such that elements of the ensemble at a larger internuclear separation are assumed to be photodissociated. The vibrational populations predicted by the QCM are compared to recent quantum simulations (Niederhausen and Thumm 2008 Phys. Rev. A 77 013404), and a remarkable agreement has been found. The applicability of this model to femtosecond and attosecond time-scale experiments is discussed and the relation to established femtochemistry and coherent control techniques are explored.

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

Characterization and parametrization in terms of atomic number of x-ray emission from K-shell filling during ion-surface interactions

When highly charged ions are incident on a surface, part of their potential energy is emitted as characteristic radiation. The energies and yields of these characteristic x rays have been measured for a series of elements at the Tokyo electron-beam ion trap. These data have been used to develop a simple model of the relaxation of the hollow atoms which are formed as the ion approaches the surface, as well as a set of semiempirical scaling laws, which allow for the ready calculation of the K-shell x-ray spectrum which would be produced by an arbitrary slow bare or hydrogenlike ion on a surface. These semiempirical scaling laws can be used to assess the merit of highly charged ion fluorescence x-ray generation in a wide range of applications.

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.

The Role of Arg228 in the phosphorylation of galactokinase: The mechanism of GHMP kinases by QM/MM studies

GHMP kinases are a group of structurally-related small molecule kinases. They have been found in all kingdoms of life and are mostly responsible for catalysing the ATP-dependent phosphorylation of intermediary metabolites. Although the GHMP kinases are of clinical, pharmaceutical and biotechnological importance, the mechanism of GHMP-kinases is controversial. A catalytic base mechanism was suggested for mevalonate kinase that has a structural feature of the ?-phosphate of ATP close to an aspartate residue; however, for one GHMP member, homoserine kinase, where the residue acting as general base is absent, a direct phosphorylation mechanism was suggested. Furthermore, it has been proposed by some authors that all the GHMP kinases function via the direct phosphorylation mechanism. This controversy in mechanism has limited our ability to exploit these enzymes as drug targets and in biotechnology. Here the phosphorylation reaction mechanism of the human galactokinase, a member of GHMP kinase was investigated using molecular dynamics simulations and density functional theory-based QM/MM calculations (B3LYP-D/AMBER99). The reaction coordinates were localized by potential energy scan using adiabatic mapping method. Our results indicate that a highly conserved Glu174 captures Arg105 to the proximity of the a-phosphate of ATP forming a H-bond network, therefore the mobility of ATP in the large oxyanion hole is restricted. Arg228 functions to stabilize the negative charge developed at the ß,?-bridging oxygen of the ATP during bond cleavage. The reaction occurs via direct phosphorylation mechanism and the Asp186 in proximity of ATP does not directly participate in the reaction pathway. Since Arg228 is not conserved among GHMP kinases, reagents which form interactions with Arg228, and therefore can interrupt its function in phosphorylation may be developed into potential selective inhibitors for galactokinase.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.

Diamond machining of silicon: A review of advances in molecular dynamics simulation

Molecular dynamics (MD) simulation has enhanced our understanding about ductile-regime machining of brittle materials such as silicon and germanium. In particular, MD simulation has helped understand the occurrence of brittle–ductile transition due to the high-pressure phase transformation (HPPT), which induces Herzfeld–Mott transition. In this paper, relevant MD simulation studies in conjunction with experimental studies are reviewed with a focus on (i) the importance of machining variables: undeformed chip thickness, feed rate, depth of cut, geometry of the cutting tool in influencing the state of the deviatoric stresses to cause HPPT in silicon, (ii) the influence of material properties: role of fracture toughness and hardness, crystal structure and anisotropy of the material, and (iii) phenomenological understanding of the wear of diamond cutting tools, which are all non-trivial for cost-effective manufacturing of silicon. The ongoing developmental work on potential energy functions is reviewed to identify opportunities for overcoming the current limitations of MD simulations. Potential research areas relating to how MD simulation might help improve existing manufacturing technologies are identified which may be of particular interest to early stage researchers.

Electrical conductivity and glass formation in nitrile-functionalized pyrrolidinium bis(trifluoromethylsulfonyl)imide ionic liquids: chain length and odd–even effects of the alkyl spacer between the pyrrolidinium ring and the nitrile group

The electrical conductivity of a series of pyrrolidinium bis(trifluoromethylsulfonyl)imide ionic liquids, functionalized with a nitrile (cyano) group at the end of an alkyl chain attached to the cation, was studied in the temperature range between 173 K and 393 K. The glass formation of the ionic liquids is influenced by the length of the alkyl spacer separating the nitrile function from the pyrrolidinium ring. The electrical conductivity and the viscosity do not show a monotonic dependence on the alkyl spacer length, but rather an odd-even effect. An explanation for this behavior is given, including the potential energy landscape picture for the glass transition.

Cyclic nanoindentation and nano-impact fatigue characterisation of functionally graded TiN/TiNi shape memory film

This paper investigates the mechanism of nanoscale fatigue of functionally graded TiN/TiNi films using nano-impact and multiple-loading-cycle nanoindentation tests. The functionally graded films were deposited on silicon substrate, in which TiNi films maintain shape memory and pseudo elastic behavior, while a modified TiN surface layer provides tribological and anti-corrosion properties. Nanomechanical tests were performed to comprehend the localized film performance and failure modes of the functionally graded film using NanoTestTM equipped with Berkovich and conical indenter between 100 μN to 500 mN loads. The loading mechanism and load history are critical to define film failure modes (i.e. backward depth deviation) including the shape memory effect of the functionally graded layer. The results are sensitive to the applied load, loading type (e.g. semi-static, dynamic) and probe geometry. Based on indentation force-depth profiles, depth-time data and post-test surface observations of films, it is concluded that the shape of the nanoindenter is critical in inducing the localized indentation stress and film failure, including shape recovery at the lower load range. Elastic-plastic finite element (FE) simulation during nanoindentation loading indicated that the location of subsurface maximum stress near the interface influences the backward depth deviation type of film failure. A standalone, molecular dynamics simulation was performed with the help of a long range potential energy function to simulate the tensile test of TiN nanowire with two different aspect ratios to investigate the theory of its failure mechanism.

Ultracold, radiative charge transfer in hybrid Yb ion – Rb atom traps

Ultracold hybrid ion–atom traps offer the possibility of microscopic manipulation of quantum coherences in the gas using the ion as a probe. However, inelastic processes, particularly charge transfer can be a significant process of ion loss and has been measured experimentally for the ${\rm Y}{{{\rm b}}^{+}}$ ion immersed in a Rb vapour. We use first-principles quantum chemistry codes to obtain the potential energy curves and dipole moments for the lowest-lying energy states of this complex. Calculations for the radiative decay processes cross sections and rate coefficients are presented for the total decay processes; ${\rm Y}{{{\rm b}}^{+}}(6{\rm s}{{\;}^{2}}{\rm S})+{\rm Rb}(5{\rm s}{{\;}^{2}}{\rm S})\to {\rm Yb}(6{{{\rm s}}^{2}}{{\;}^{1}}{\rm S})+{\rm R}{{{\rm b}}^{+}}(4{{{\rm p}}^{6}}{{\;}^{1}}{\rm S})+h\nu $ and ${\rm Y}{{{\rm b}}^{+}}(6{\rm s}{{\;}^{2}}{\rm S})+{\rm Rb}(5{\rm s}{{\;}^{2}}{\rm S})\to {\rm YbR}{{{\rm b}}^{+}}({{X}^{1}}{{\Sigma }^{+}})+h\nu $. Comparing the semi-classical Langevin approximation with the quantum approach, we find it provides a very good estimate of the background at higher energies. The results demonstrate that radiative decay mechanisms are important over the energy and temperature region considered. In fact, the Langevin process of ion–atom collisions dominates cold ion–atom collisions. For spin-dependent processes [1] the anisotropic magnetic dipole–dipole interaction and the second-order spin–orbit coupling can play important roles, inducing coupling between the spin and the orbital motion. They measured the spin-relaxing collision rate to be approximately five orders of magnitude higher than the charge-exchange collision rate [1]. Regarding the measured radiative charge transfer collision rate, we find that our calculation is in very good agreement with experiment and with previous calculations. Nonetheless, we find no broad resonances features that might underly a strong isotope effect. In conclusion, we find, in agreement with previous theory that the isotope anomaly observed in experiment remains an open question.

Radiative charge transfer in cold and ultracold sulphur atoms colliding with protons

Radiative decay processes at cold and ultra cold temperatures for sulfur atoms colliding with protons are investigated. The MOLPRO quantum chemistry suite of codes was used to obtain accurate potential energies and transition dipole moments, as a function of internuclear distance, between low-lying states of the SH+ molecular cation. A multi-reference configuration-interaction approximation together with the Davidson correction is used to determine the potential energy curves and transition dipole moments, between the states of interest, where the molecular orbitals are obtained from state-averaged multi-configuration-self-consistent field calculations. The collision problem is solved approximately using an optical potential method to obtain radiative loss, and a fully two-channel quantum approach for radiative charge transfer. Cross sections and rate coefficients are determined for the first time for temperatures ranging from 10 μK up to 10 000 K. Results are obtained for all isotopes of sulfur, colliding with H+ and D+ ions and comparison is made to a number of other collision systems.