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.

Response of soil microbial biomass to 1,2-dichlorobenzene addition in the presence of plant residues

The impact of 1,2-dichlorobenzene on soil microbial biomass in the presence and absence of fresh plant residues (roots) was investigated by assaying total vital bacterial counts, vital fungel hyphal length, total culturable bacterial counts, and culturable fluorescent pseudomonads. Diversity of the fluorescent pseudomonads was investigated using fatty acid methyl ester (FAME) characterization in conjunction with metabolic profiling of the sampled culturable community (Biolog). Mineralization of [¹⁴C]1,2- dichlorobenzene was also assayed. Addition of fresh roots stimulated 1,2- dichlorobenzene mineralization by over 100%, with nearly 20% of the label mineralized in root-amended treatments by the termination of the experiment. Presence of roots also buffered any impacts of 1,2-dichlorobenzene on microbial numbers. In the absence of roots, 1,2-dichlorobenzene greatly stimulated total culturable bacteria and culturable pseudomonads in a concentration-dependent manner. 1,2-Dichlorobenzene, up to concentrations of 50 μg/g soil dry weight had little or no deleterious effects on microbial counts. The phenotypic diversity of the fluorescent pseudomonad population was unaffected by the treatments, even though fluorescent pseudomonad numbers were greatly stimulated by both roots and 1,2-dichlorobenzene. The presence of roots had no detectable impact on the bacterial community composition. No phenotypic shifts in the natural population were required to benefit from the presence of roots and 1,2-dichlorobenzene. The metabolic capacity of the culturable bacterial community was altered in the presence of roots but not in the presence of 1,2-dichlorobenzene. It is argued that the increased microbial biomass and shifts in metabolic capacity of the microbial biomass are responsible for enhanced degradation of 1,2-dichlorobenzene in the presence of decaying plant roots.

Toxic interactions of metal ions (Cd2+, Pb2+, Zn2+ and Sb3- on in vitro biomass production of ectomycorrhizal fungi

A number of ectomycorrhizal (ECM) fungi, from sites uncontaminated by toxic metals, were investigated to determine their sensitivity to Cd^2-, Pb²⁺, Zn²⁺ and Sb^3-, measured as an inhibition of fungal biomass production. Isolates were grown in liquid media amended with the metals, individually (over a range of concentrations) and in combination (at single concentrations) to determine any significant interactions between the metals. Significant interspecific variation in sensitivity to Cd²⁺ and Zn²⁺ was recorded, while Pb²⁺ and Sb^3- individually had little effect. The presence of Pb²⁺ and Sb^3- in the media did however, ameliorate Cd²⁺ and Zn²⁺ toxicity in some circumstances. Interactions between Cd²⁺ and Zn²⁺ were investigated further over a range of concentrations. Zn²⁺ was found to significantly ameliorate the toxicity of Cd²⁺ to three of the four isolates tested. The influence of Zn²⁺ varied between ECM species and with the concentrations of metals tested.

Biomass allocation, phosphorus nutrition and vesicular-arbuscular mycorrhizal infection in clones of Yorkshire Fog, Holcus lanatus L. (Poaceae) that differ in their phosphate uptake kinetics and tolerance to arsenate

Biomass and phosphorus allocation were determined in arsenate tolerant and non-tolerant clones of the grass Holcus lanatus L. in both solution culture and in soil. Arsenate is a phosphate analogue and is taken up by the phosphate uptake system. Tolerance to arsenate in this grass is achieved by suppression of arsenate (and phosphate) influx. When clones differing in their arsenate tolerance were grown in solution culture with a range of phosphate levels, a tolerant clone did not fare as well as a non-tolerant at low levels of phosphate nutrition in that it had reduced shoot biomass production, increased biomass allocation to the roots and lower shoot phosphorus concentration. At a higher level of phosphate nutrition there was little or no difference in these parameters, suggesting that differences at lower levels of phosphate nutrition were due solely to differences in the rates of phosphate accumulation. In experiments in sterile soil (potting compost) the situation was more complicated with tolerant plants having lower growth rates but higher phosphorus concentrations. The gene for arsenate tolerance is polymorphic in arsenate uncontaminated populations. When phosphorus concentration of tolerant phenotypes was determined in one such population, again tolerants had a higher phosphorus status than non-tolerants. Tolerants also had higher rates of vesicular-arbuscular mycorrhizal (VAM) infection. The ecological implications of these results are that it appears that suppression of the high affinity uptake system, is at least in part, compensated by increased mycorrhizal infection. © 1994 Kluwer Academic Publishers.

Techno-Economic and Life Cycle Assessment on Lignocellulosic Biomass Thermochemical Conversion Technologies: A Review

Bioenergy derived from biomass provides a promising energy alternative and can reduce the greenhouse gas (GHG) emissions generated from fossil fuels. Biomass-based thermochemical conversion technologies have been acknowledged as apt options to convert bioresources into bioenergy; this bioenergy includes electricity, heat, and fuels/chemicals in solid, liquid, and gaseous phases. In this review, the techno-economic and life cycle assessment of these technologies (combustion, gasification, pyrolysis, liquefaction, carbonization, and co-firing) are summarized. Specific indicators (production costs in a techno-economic analysis, functional units and environmental impacts in a life cycle analysis) for different technologies were compared. Finally, gaps in research and future trends in biomass thermochemical conversion were identified. This review could be used to guide future research related to economic and environmental benefits of bioenergy.

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.

Estimating biomass, fishing mortality, and "total allowable discards" for surveyed non-target fish

Demersal fisheries targeting a few high-value species often catch and discard other "non-target" species. It is difficult to quantify the impact of this incidental mortality when population biomass of a non-target species is unknown. We calculate biomass for 14 demersal fish species in ICES Area VIIg (Celtic Sea) by applying species-and length-based catchability corrections to catch records from the Irish Groundfish Survey (IGFS). We then combine these biomass estimates with records of commercial discards (and landings for marketable non-target species) to calculate annual harvesting rates (HR) for each study species. Uncertainty is incorporated into estimates of both biomass andHR. Our survey-based HR estimates for cod and whiting compared well with HR-converted fishing mortality (F) estimates from analytical assessments for these two stocks. Of the non-target species tested, red gurnard (Chelidonichthys cuculus) recorded some annual HRs greater than those for cod or whiting; challenging "Pope's postulate" that F on non-target stocks in an assemblage will not exceed that on target stocks. We relate HR for each species to two corresponding maximum sustainable yield (MSY) reference levels; six non-target species (including three ray species) show annual HRs >= HRMSY. This result suggests that it may not be possible to conserve vulnerable non-target species when F is coupled to that of target species. Based on biomass, HR, and HRMSY, we estimate "total allowable catch" for each non-target species.

Changes in epilithic biomasses and invertebrate community structure over a deposit metal concentration gradient in upland headwater streams

Stream bed metal deposits affect the taxon richness, density and taxonomic diversity of primary and secondary producers by a variety of direct or indirect abiotic and biotic processes but little is known about the relative importance of these processes over a deposit metal concentration gradient. Inorganic matter (IM), algal and non-photosynthetic detrital (NPD) dry biomasses were estimated for 10 monthly samples, between 2007 and 2008, from eight sites differing in deposit density. Invertebrate abundance, taxon richness and composition were also determined. Relations between these variables were investigated by canonical correspondence analysis (CCA), generalized estimating equation models and path analysis. The first CCA axis correlates with deposit density and invertebrate abundance, with lumbriculids and chironomids increasing in abundance with deposit density and all other taxa declining. Community structure changes significantly above a deposit density of approximately 8 mg cm, when algal biomass, invertebrate richness and diversity decline. Invertebrate richness and diversity were determined by direct effects of NPD biomass and indirect effects of IM. Algal biomass only had an effect on invertebrate abundance. Possible pH, oxygen, food and ecotoxicological effects of NPD biomass on the biota are discussed.

Biomass-derived oxymethylene ethers as diesel additives: A techno-economic assessment

Conversion of agricultural biomass such as wood chips, wheat straw and forest residue for the production of fuels can help in reducing GHG emissions since they are considered as nearly carbon neutral. Around the world there is a significant amount of forest and agricultural-biomass available which could be used for the production of liquid fuels that can be blended with the petroleum-based diesel. Oxymethylene ethers (OMEs) can be derived from biomass via gasification, water-gas shift reaction and methanol production. The addition of OMEs to conventional diesel fuel has great potential to reduce soot formation during the combustion in diesel engines. Unlike methanol and dimethyl ether (DMM) which can also reduce soot formation, the physical properties of OMEs allow the use in modern diesel engines without significant change of the engines infrastructure. In this study, a detailed and data intensive process simulation model was developed to simulate all the unit operations involved in the production of OMEs from biomass. The unit operation considered include biomass drying, gasification, gas cleaning, water gas shift reaction, methanol production and OMEs synthesis. The simulation results were then utilized to conduct a detailed techno-economic assessment study of the whole biomass conversion chain to determine the most attractive pathways for OMEs production. Our recent study shows that the key parameters affecting the OMEs production are equivalence ratio, H2/CO ratio and optimal air flow. Overall, the cost of production ($/liter) of OMEs from different biomass feedstock in Alberta will be determined

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.