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.

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.

Cosmological Constraints from Measurements of Type Ia Supernovae Discovered During the First 1.5 Yr of the Pan-STARRS1 Survey

We present griz_P1 light curves of 146 spectroscopically confirmed Type Ia supernovae (SNe Ia; 0.03 < z < 0.65) discovered during the first 1.5 yr of the Pan-STARRS1 Medium Deep Survey. The Pan-STARRS1 natural photometric system is determined by a combination of on-site measurements of the instrument response function and observations of spectrophotometric standard stars. We find that the systematic uncertainties in the photometric system are currently 1.2% without accounting for the uncertainty in the Hubble Space Telescope Calspec definition of the AB system. A Hubble diagram is constructed with a subset of 113 out of 146 SNe Ia that pass our light curve quality cuts. The cosmological fit to 310 SNe Ia (113 PS1 SNe Ia + 222 light curves from 197 low-z SNe Ia), using only supernovae (SNe) and assuming a constant dark energy equation of state and flatness, yields w = -1.120^+0.360_-0.206(Stat)^+0.269_0.291(Sys). When combined with BAO+CMB(Planck)+H₀, the analysis yields Ω_M = 0.280^+0.013_0.012 and w = -1.166^+0.072_-0.069 including all identified systematics. The value of w is inconsistent with the cosmological constant value of -1 at the 2.3σ level. Tension endures after removing either the baryon acoustic oscillation (BAO) or the H₀ constraint, though it is strongest when including the H₀ constraint. If we include WMAP9 cosmic microwave background (CMB) constraints instead of those from Planck, we find w = -1.124^+0.083_-0.065, which diminishes the discord to <2σ. We cannot conclude whether the tension with flat ΛCDM is a feature of dark energy, new physics, or a combination of chance and systematic errors. The full Pan-STARRS1 SN sample with ∼three times as many SNe should provide more conclusive results.

Solutions of fourth-order parabolic equation modeling thin film growth

In this paper we study the well-posedness for a fourth-order parabolic equation modeling epitaxial thin film growth. Using Kato's Method [1], [2] and [3] we establish existence, uniqueness and regularity of the solution to the model, in suitable spaces, namelyC⁰([0,T];L^p(Ω)) where  with 1<α<2, n∈N and n≥2. We also show the global existence solution to the nonlinear parabolic equations for small initial data. Our main tools are L_p–L_q-estimates, regularization property of the linear part of e^−tΔ2 and successive approximations. Furthermore, we illustrate the qualitative behavior of the approximate solution through some numerical simulations. The approximate solutions exhibit some favorable absorption properties of the model, which highlight the stabilizing effect of our specific formulation of the source term associated with the upward hopping of atoms. Consequently, the solutions describe well some experimentally observed phenomena, which characterize the growth of thin film such as grain coarsening, island formation and thickness growth.

On the numerical solution of the 2d wave equation with compact fdtd schemes

This paper discusses compact-stencil finite difference time domain (FDTD) schemes for approximating the 2D wave equation in the context of digital audio. Stability, accuracy, and efficiency are investigated and new ways of viewing and interpreting the results are discussed. It is shown that if a tight accuracy constraint is applied, implicit schemes outperform explicit schemes. The paper also discusses the relevance to digital waveguide mesh modelling, and highlights the optimally efficient explicit scheme.