Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

External validation of the 2008 Framingham cardiovascular risk equation for CHD and stroke events in a European population of middle-aged men. The PRIME study

To test the applicability of the sex-specific 2008 Framingham general cardiovascular risk equation for coronary heart disease (CHD) and stroke in European middle-aged men from Ireland and France.

Material Logistics in Urban Construction Sites: A Structural Equation Model

Aim: The aim of this paper is to identify best practice relating to the effective management of materials in an urban, confined construction site, using structural equation modelling.

Methodology: A literature review, case study analysis and questionnaire survey are employed, with the results scrutinised using confirmatory factor analysis in the form of structural equation modelling.

Results: The following are the leading strategies in the management of materials in a confined urban site environment; (1) Consult and review the project programme, (2) Effective communication and delivery, (3) Implement site safety management plans, and (4) Proactive spatial monitoring and control.

Implication for Practice: With the relentless expansion of urban centres and the increasing high cost of materials, any potential savings made on-site would translate into significant monetary concessions on completion of a development.

Originality/Value: As on-site project management professionals successfully identify and implement the various strategies in the management of plant and materials on a confined urban site, successful resource management in this restrictive environment is attainable.

Innovative Aspect of Paper: An empirical study of three different construction sites in three different countries (Ireland, England and USA) together with a questionnaire survey from the industry, investigating the managerial strategies in the management of plant and material in confined urban site environments

SOLVING THE INHOMOGENEOUS SCHRODINGER-EQUATION ON A BIDIRECTIONAL PIPELINE OF TRANSPUTERS

In this paper we describe the design of a parallel solution of the inhomogeneous Schrodinger equation, which arises in the construction of continuum orbitals in the R-matrix theory of atomic continuum processes. A prototype system is described which has been programmed in occam2 and implemented on a bi-directional pipeline of transputers. Some timing results for the prototype system are presented, and the development of a full production system is discussed.

Difference equation approach to two-thermocouple sensor characterization in constant velocity flow environments

Thermocouples are one of the most popular devices for temperature measurement due to their robustness, ease of manufacture and installation, and low cost. However, when used in certain harsh environments, for example, in combustion systems and engine exhausts, large wire diameters are required, and consequently the measurement bandwidth is reduced. This article discusses a software compensation technique to address the loss of high frequency fluctuations based on measurements from two thermocouples. In particular, a difference equation (DE) approach is proposed and compared with existing methods both in simulation and on experimental test rig data with constant flow velocity. It is found that the DE algorithm, combined with the use of generalized total least squares for parameter identification, provides better performance in terms of time constant estimation without any a priori assumption on the time constant ratios of the thermocouples.

Analysis of deformation behavior and workability of advanced 9Cr-Nb-V ferritic heat resistant steels

Hot compression tests were carried out on 9Cr–Nb–V heat resistant steels in the temperature range of 600–1200 °C and the strain rate range of 10−2–100 s−1 to study their deformation characteristics. The full recrystallization temperature and the carbon-free bainite phase transformation temperature were determined by the slope-change points in the curve of mean flow stress versus the inverse of temperature. The parameters of the constitutive equation for the experimental steels were calculated, including the stress exponent and the activation energy. The lower carbon content in steel would increase the fraction of precipitates by increasing the volume of dynamic strain-induced (DSIT) ferrite during deformation. The ln(εc) versus ln(Z) and the ln(σc) versus ln(Z) plots for both steels have similar trends. The efficiency of power dissipation maps with instability maps merged together show excellent workability from the strain of 0.05 to 0.6. The microstructure of the experimental steels was fully recrystallized upon deformation at low Z value owing to the dynamic recrystallization (DRX), and exhibited a necklace structure under the condition of 1050 °C/0.1 s−1 due to the suppression of the secondary flow of DRX. However, there were barely any DRX grains but elongated pancake grains under the condition of 1000 °C/1 s−1 because of the suppression of the metadynamic recrystallization (MDRX).

Constitutive modeling, microstructure evolution, and processing map for a nitride-strengthened heat-resistant steel

A constitutive equation was established to describe the deformation behavior of a nitride-strengthened (NS) steel through isothermal compression simulation test. All the parameters in the constitutive equation including the constant and the activation energy were precisely calculated for the NS steel. The result also showed that from the stress-strain curves, there existed two different linear relationships between critical stress and critical strain in the NS steel due to the augmentation of auxiliary softening effect of the dynamic strain-induced transformation. In the calculation of processing maps, with the change of Zener-Hollomon value, three domains of different levels of workability were found, namely excellent workability region with equiaxed-grain microstructure, good workability region with “stripe” microstructure, and the poor workability region with martensitic-ferritic blend microstructure. With the increase of strain, the poor workability region first expanded, then shrank to barely existing, but appeared again at the strain of 0.6.

Kinetics study on thermal dissociation of levoglucosan during cellulose pyrolysis

(Chemical Equation Presented) The mechanisms and kinetics studies of the levoglucosan (LG) primary decomposition during cellulose pyrolysis have been carried out theoretically in this paper. Three decomposition mechanisms (C-O bond scission, C-C bond scission, and LG dehydration) including nine pathways and 16 elementary reactions were studied at the B3LYP/6-31 + G(D,P) level based on quantum mechanics. The variational transi-tion- state rate constants for every elementary reaction and every pathway were calculated within 298-1550 K. The first-order Arrhenius expressions for these 16 elementary reactions and nine pathways were suggested. It was concluded that computational method using transition state theory (TST) without tunneling correction gives good description for LG decomposition by comparing with the experimental result. With the temperature range of 667-1327 K, one dehydration pathway, with one water molecule composed of a hydrogen atom from C3 and a hydroxyl group from C2, is a preferred LG decomposition pathway by fitting well with the experimental results. The calculated Arrhenius plot of C-O bond scission mechanism is better agreed with the experimental Arrhenius plot than that of C-C bond scission. This C-O bond scission mechanism starts with breaking of C1-O5 and C6-O1 bonds with formation of CO molecule (C1-O1) simultaneously. C-C bond scission mechanism is the highest energetic barrier pathway for LG decomposition. © 2013 Elsevier Ltd. All rights reserved.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Modeling relativistic soliton interactions in overdense plasmas: A perturbed nonlinear Schrodinger equation framework

We investigate the dynamics of localized solutions of the relativistic cold-fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed nonlinear Schrodinger equation that describes the evolution of the envelope of circularly polarized electromagnetic field. Retaining terms up to fifth order in the small perturbation parameter, we derive a self-consistent framework for the description of the plasma response in the presence of localized electromagnetic field. The formalism is applied to standing electromagnetic soliton interactions and the results are validated by simulations of the full cold-fluid model. To lowest order, a cubic nonlinear Schrodinger equation with a focusing nonlinearity is recovered. Classical quasiparticle theory is used to obtain analytical estimates for the collision time and minimum distance of approach between solitons. For larger soliton amplitudes the inclusion of the fifth-order terms is essential for a qualitatively correct description of soliton interactions. The defocusing quintic nonlinearity leads to inelastic soliton collisions, while bound states of solitons do not persist under perturbations in the initial phase or amplitude

GPU Acceleration of Partial Differential Equation Solvers

Differential equations are often directly solvable by analytical means only in their one dimensional version. Partial differential equations are generally not solvable by analytical means in two and three dimensions, with the exception of few special cases. In all other cases, numerical approximation methods need to be utilized. One of the most popular methods is the finite element method. The main areas of focus, here, are the Poisson heat equation and the plate bending equation. The purpose of this paper is to provide a quick walkthrough of the various approaches that the authors followed in pursuit of creating optimal solvers, accelerated with the use of graphical processing units, and comparing them in terms of accuracy and time efficiency with existing or self-made non-accelerated solvers.