Criteria for local equilibrium in a system with transport of heat and mass

Nonequilibrium molecular dynamics is used to compute the coupled heat and mass transport in a binary isotope mixture of particles interacting with a Lennard-Jones/spline potential. Two different stationary states are studied, one with a fixed internal energy flux and zero mass flux, and the other with a fixed diffusive mass flux and zero temperature gradient. Computations are made for one overall temperature,T=2, and three overall number densities,n=0.1, 0.2, and 0.4. (All numerical values are given in reduced, Lennard-Jones units unless otherwise stated.) Temperature gradients are up to T=0.09 and weight-fraction gradients up to w ₁=0.007. The flux-force relationships are found to be linear over the entire range. All four transport coefficients (theL-matrix) are determined and the Onsager reciprocal relationship for the off-diagonal coefficients is verified. Four different criteria are used to analyze the concept of local equilibrium in the nonequilibrium system. The local temperature fluctuation is found to be T0.03T and of the same order as the maximum temperature difference across the control volume, except near the cold boundary. A comparison of the local potential energy, enthalpy, and pressure with the corresponding equilibrium values at the same temperature, density, and composition also verifies that local equilibrium is established, except near the boundaries of the system. The velocity contribution to the BoltzmannH-function agrees with its Maxwellian (equilibrium) value within 1%, except near the boundaries, where the deviation is up to 4%. Our results do not support the Eyring-type transport theory involving jumps across energy barriers; we find that its estimates for the heat and mass fluxes are wrong by at least one order of magnitude.     

Two-proton correlations in the target fragmentation region of nuclear collisions at 200A GeV

Correlations between protons are studied in the target fragmentation region of reactions of protons and¹⁶O with C, Cu, Ag, Au and of³²S with Al and Au at 200A GeV. The emitted protons were measured with the Plastic Ball detector in the WA80 experiment at the CERN SPS. The comparison of the correlation function with calculations, assuming a spherical, gaussian shaped source with a lifetime τ=0 fm/c, allows the extraction of radius parameters. The values are very close to those expected from the geometry of the target nuclei and increase with the target mass as αA _Target ^1/3 . Even in proton induced reactions the whole target nucleus is involved. The dependence of the radii on centrality, polar angleθ _lab, and energy, and their relation to measured proton yields are presented.     

Development of a very low energy μ+ beam at PSI

At PSI we are investigating the technique of decelerating an existing very intense secondary beam of surface ⁺ (4 MeV) to an energy of 10 eV using appropriate moderators. These ⁺ can then be used as a source of a tertiary beam of low energy muons with tunable kinetic energy between 10 eV and 10 keV.With a 1000 A layer of solid Argon deposited on an Al substrate we obtain a moderation efficiency (with respect to the number of incoming surface ⁺) of the order of 10^–4.Results of our investigations and the present status of the project are presented together with future plans and possibilities.     

Transient flow towards a well in an aquifer including the effect of fluid inertia

Transient propagation of weak pressure perturbations in a homogeneous, isotropic, fluid saturated aquifer has been studied. A damped wave equation for the pressure in the aquifer is derived using the macroscopic, volume averaged, mass conservation and momentum equations. The equation is applied to the case of a well in a closed aquifer and analytical solutions are obtained to two different flow cases. It is shown that the radius of influence propagates with a finite velocity. The results show that the effect of fluid inertia could be of importance where transient flow in porous media is studied.List of symbols b Thickness of the aquifer, m - c ₀ Wave velocity, m/s - k Permeability of the porous medium, m² - n Porosity of the porous medium - p( ,t) Pressure, N/m² - Q Volume flux, m³/s - r Radial coordinate, m - r _w Radius of the well, m - s Transform variable - S Storativity of the aquifer - S _d(r, t) Drawdown, m - t Time, s - T Transmissivity of the aquifer, m²/s - ( ,t) Velocity of the fluid, m/s - Coordinate vector, m - z Vertical coordinate, m - Coefficient of compressibility, m²/N - Coefficient of fluid compressibility, m²/N - Relaxation time, s - (r, t) Hydraulic potential, m - Dynamic viscosity of the fluid, Ns/m² - Dimensionless radius - Density of the fluid, Ns²/m⁴ - (, ) Dimensionless drawdown - Dimensionless time - , x Dummy variables - ₀, ₁ Auxilary functions     

Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations

Molecular dynamics simulations have been performed to gain insights into the catalytic mechanism of the hydrolysis of epoxides to vicinal diols by soluble epoxide hydrolase (sEH). The binding of a substrate, 1S,2S-trans-methylstyrene oxide, was studied in two conformations in the active site of the enzyme. It was found that only one is likely to be found in the active enzyme. In the preferred conformation the phenyl group of the substrate is pi-sandwiched between two aromatic residues, Tyr381 and His523, whereas the other conformation is pi-stacked with only one aromatic residue, Trp334. Two simulations were carried out to 1 ns for each conformation to evaluate the protonation state of active site residue His523. It was found that a protonated histidine is essential for keeping the active site from being disrupted. Long time scale, 4 ns, molecular dynamics simulation was done for the structure with the most likely combination of binding conformation and protonation state of His523. Near Attack Conformers (NACs) are present 5.3% of the time and nucleophilic attack on either epoxide carbon atom, approximately 75% on C(1) and approximately 25% on C(2), is found. A maximum of one hydrogen bond between the epoxide oxygen and either of the active site tyrosines, Tyr465 and Tyr381, is present, in agreement with experimental mutagenesis results that reveal a slight loss in activity if one tyrosine is mutated and essential loss of all activity upon double mutation of the two tyrosines in question. It was found that a hydrogen bond from Tyr465 to the substrate oxygen is essential for controlling the regioselectivity of the reaction. Furthermore, a relationship between the presence of this hydrogen bond and the separation of reactants was found. Two groups of amino acid segments were identified each as moving collectively. Furthermore, an overall anti-correlation was found between the movements of these two individually collectively moving groups, made up by parts of the cap-region, including the two tyrosines, and the site of the catalytic triad, respectively. This overall anti-correlated collective domain motion is, perhaps, involved in the conversion of E.NAC to E.TS.     

Toward separation of nuclear spin isomers with coherent light.

We propose an approach for separating nuclear spin isomers with coherent light and illustrate it by numerical calculations using fulvene as a model system. The scheme employs the equivalence of torsion and interchange of equivalent H-atoms in a class of molecules of which fulvene is a simple example. The exchange symmetry couples with the rotational symmetry to produce a spatial distinction between the two photo-excited nuclear spin isomers, and wavepacket interferometry is applied to separate the species.