Molecular Dynamics Simulation of Methane Hydrate Dissociation Process in the Presence of Thermodynamic Inhibitor

The dissociation of methane hydrate in the presence of ethylene glycol (11.45 mol.L-1) at 277.0 K was studied using canonical ensemble (NVT) molecular dynamics simulations. Results show that hydrate dissociation starts from the surface layer of the solid hydrate and then gradually expands to the internal layer. Thus, the solid structure gradually shrinks until it disappears. A distortion of the hydrate lattice structure occurs first and then the hydrate evolves from a fractured frame to a fractional fragment. Finally, water molecules in the hydrate construction exist in the liquid state. The inner dissociating layer is, additionally, coated by a liquid film formed from outer dissociated water molecules outside. This film inhibits the mass transfer performance of the inner molecules during the hydrate dissociation process.

Self-assembled quantum dots, wires and quantum-dot lasers

Molecular beam epitaxy-grown self-assembled In(Ga)As/GaAs and InAs/InAlAs/InP quantum dots (QDs) and quantum wires (QWRs) have been studied. By adjusting growth conditions, surprising alignment. preferential elongation, and pronounced sequential coalescence of dots and wires under specific condition are realized. The lateral ordering of QDs and the vertical anti-correlation of QWRs are theoretically discussed. Room-temperature (RT) continuous-wave (CW) lasing at the wavelength of 960 nm with output power of 3.6 W from both uncoated facets is achieved fi-om vertical coupled InAs/GaAs QDs ensemble. The RT threshold current density is 218 A/cm(2). A RT CW output power of 0.6 W/facet ensures at least 3570 h lasing (only drops 0.83 dB). (C) 2001 Elsevier Science B.V, All rights reserved.

Influence of interdot electronic coupling on photoluminescence spectra of self-assembled InAs/GaAs quantum dots

The influence of interdot electronic coupling on photoluminescence (PL) spectra of self-assembled InAs/GaAs quantum dots (QDs) has been systematically investigated combining with the measurement of transmission electron microscopy. The experimentally observed fast red-shift of PL energy and an anomalous reduction of the linewidth with increasing temperature indicate that the QD ensemble can be regarded as a coupled system. The study of multilayer vertically coupled QD structures shows that a red-shift of PL peak energy and a reduction of PL linewidth are expected as the number of QD layers is increased. On the other hand, two layer QDs with different sizes have been grown according to the mechanism of a vertically correlated arrangement. However, only one PL peak related to the large QD ensemble has been observed due to the strong coupling in InAs pairs. A new possible mechanism to reduce the PL linewidth of QD ensemble is also discussed.

Molecular dynamics study on the phase diagrams of linear and branched chain molecules

The particle transfer molecular dynamics is used to study the phase equilibria of linear and branched chain molecules. The scaling of the critical temperature versus chain length is obtained and the critical densities are found to decrease with increasing chain length, which are in agreement with the results of experiment and theory. The phase diagrams of the linear and the branched chain molecules nearly overlap with each other. Moreover, the radial distribution functions of linear and branched chain molecules in gas phase are very similar, but in the liquid phase, they are different for different kinds of chains.

Effect of annealing on the viscoelastic and viscoplastic responses of low-density polyethylene

Two series of tensile tests with constant crosshead speeds (ranging from 5 to 200 mm/min) and tensile relaxation tests (at strains from 0.03 to 0.09) were performed on low-density polyethylene in the subyield region of deformations at room temperature. Mechanical tests were carried out on nonannealed specimens and on samples annealed for 24 h at the temperatures T = 50, 60, 70, 80, and 100 degreesC. Constitutive equations were derived for the time-dependent response of semicrystalline polymers at isothermal deformations with small strains. A polymer is treated as an equivalent heterogeneous network of chains bridged by temporary junctions (entanglements, physical crosslinks, and lamellar blocks). The network is thought of as an ensemble of mesoregions linked with each other. The viscoelastic behavior of a polymer is modeled as a thermally induced rearrangement of strands (separation of active strands from temporary junctions and merging of dangling strands with the network). The viscoplastic response reflects sliding of junctions in the network with respect to their reference positions driven by macrostrains. Stress-strain relations involve five material constants that were found by fitting the observations.

The viscoelastic and viscoplastic behavior of low-density polyethylene

Three series of tensile tests with constant cross-head speeds (ranging from 5 to 200 mm/min), tensile relaxation tests (at strains from 0.03 to 0.09) and tensile creep tests (at stresses from 2.0 to 6.0 MPa) are performed on low-density polyethylene at room temperature. Constitutive equations are derived for the time-dependent response of semicrystalline polymers at isothermal deformation with small strains. A polymer is treated as an equivalent heterogeneous network of chains bridged by temporary junctions (entanglements, physical cross-links and lamellar blocks). The network is thought of as an ensemble of meso-regions linked with each other. The viscoelastic behavior of a polymer is modelled as thermally-induced rearrangement of strands (separation of active strands from temporary junctions and merging of dangling strands with the network). The viscoplastic response reflects mutual displacement of meso-domains driven by macro-strains. Stress-strain relations for uniaxial deformation are developed by using the laws of thermodynamics. The governing equations involve five material constants that are found by fitting the observations. Fair agreement is demonstrated between the experimental data and the results of numerical simulation.

Single-molecule dynamics reveals cooperative binding-folding in protein recognition

The study of associations between two biomolecules is the key to understanding molecular function and recognition. Molecular function is often thought to be determined by underlying structures. Here, combining a single-molecule study of protein binding with an energy-landscape-inspired microscopic model, we found strong evidence that biomolecular recognition is determined by flexibilities in addition to structures. Our model is based on coarse-grained molecular dynamics on the residue level with the energy function biased toward the native binding structure ( the Go model). With our model, the underlying free-energy landscape of the binding can be explored. There are two distinct conformational states at the free-energy minimum, one with partial folding of CBD itself and significant interface binding of CBD to Cdc42, and the other with native folding of CBD itself and native interface binding of CBD to Cdc42. This shows that the binding process proceeds with a significant interface binding of CBD with Cdc42 first, without a complete folding of CBD itself, and that binding and folding are then coupled to reach the native binding state.

Temperature dependence of the distribution of the first passage time: Results from discontinuous molecular dynamics simulations of an all-atom model of the second beta-hairpin fragment of protein G

More than 22 000 folding kinetic simulations were performed to study the temperature dependence of the distribution of first passage time (FPT) for the folding of an all-atom Go-like model of the second beta-hairpin fragment of protein G. We find that the mean FPT (MFPT) for folding has a U (or V)-shaped dependence on the temperature with a minimum at a characteristic optimal folding temperature T-opt*. The optimal folding temperature T-opt* is located between the thermodynamic folding transition temperature and the solidification temperature based on the Lindemann criterion for the solid. Both the T-opt* and the MFPT decrease when the energy bias gap against nonnative contacts increases. The high-order moments are nearly constant when the temperature is higher than T-opt* and start to diverge when the temperature is lower than T-opt*. The distribution of FPT is close to a log-normal-like distribution at T* greater than or equal to T-opt*. At even lower temperatures, the distribution starts to develop long power-law-like tails, indicating the non-self-averaging intermittent behavior of the folding dynamics. It is demonstrated that the distribution of FPT can also be calculated reliably from the derivative of the fraction not folded (or fraction folded), a measurable quantity by routine ensemble-averaged experimental techniques at dilute protein concentrations.

Studies of electrochemical quantized capacitance charging of surface ensembles of silver nanoparticles

CTAB-stabilized silver nanoparticles were synthesized by NaBH4 reduction. The as-prepared nanoparticles can be self-assembled on 3-mercaptopropionic acid (MPA) modified gold electrode, which was supported strongly by XPS measurements. Exceptional long-term stability of the as-prepared colloidal silver aqueous solution and the desorption of silver nanoparticle ensemble on MPA after alcohol rinsing proved that these CTAB molecules adsorbed on silver core formed interdigitated bilayer structure. DPV and differential capacitance measurements were performed to characterize the as-prepared silver nanoparticle ensemble. and the interesting quantized capacitance charging behaviors were observed.

Effective response of a non-linear composite in external AC electric field

A general effective response is proposed for nonlinear composite media, which obey a current field relation of the form J = sigmaE + chi\E\(2) E when an external alternating current (AC) electrical field is applied. For a sinusoidal applied field with finite frequency omega, the effective constitutive relation between the current density and electric field can be defined as, = sigma(e) + chi(e) <\E(x, omega, t)\(2) E(x, omega, t)> + (. . .), where sigma(e) and chi(e) are the general effective linear and nonlinear conductive responses, respectively. The angled brackets <(. . .)> denotes the ensemble average. As two examples, we have investigated the cylindrical and spherical inclusions embedded in a host and also derived the formulae of the general effective linear and nonlinear conductive responses in dilute limit. For higher volume fraction of inclusions, we have proposed a nonlinear effective medium approximation (EMA) method to estimate the general effective response of nonlinear composites in external AC field. Furthermore, the effective nonlinear responses at harmonics are predicted by using the general effective response. (C) 2002 Elsevier Science B.V. All rights reserved.

A hybrid coupled ocean-atmosphere model and ENSO prediction study

A hybrid coupled ocean-atmosphere model is designed, which consists of a global AGCM and a simple anomaly ocean model in the tropical Pacific. Retroactive experimental predictions initiated in each season from 1979 to 1994 are performed. Analyses indicate that: (1) The overall predictive capability of this model for SSTA over the central-eastern tropical Pacific can reach one year, and the error is not larger than 0.8 degrees C. (2) The prediction skill depends greatly on the season when forecasts start. However, the phenomenon of SPB (spring prediction barrier) is not found in the model. (3) The ensemble forecast method can effectively improve prediction results. A new initialization scheme is discussed.

Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.

东中国海区域海浪同化系统设计与研究

针对目前海浪同化中没有适合东中国海区的业务应用系统，以及缺乏集合卡尔曼滤波方法（Ensemble Kalman Filter，EnKF）的应用的现状，设计了东中国海区域海浪同化系统。首先比较了目前常用的各种资料同化方法，指出各种同化方法的实质都是一种滤波过程，并选取最优插值方法（Optimal Interpolation，OI）和EnKF方法开展同化试验，动力模式选取WAVEWATCH III，观测资料为Topex/Poseidon卫星高度计观测波高。然后用观测法研究了模式预报误差协方差的统计性质，指出误差相关距离尺度在3°至6°之间。最后用22001号浮标观测资料验证了两个同化系统2000年8月的有效波高计算值，结果表明OI方案的同化系统对有效波高的均方根误差减少了9.0%，系统运行稳定，可应用于业务化部门；EnKF方案的同化系统在集合样本数为50的情况下，对有效波高的均方根误差减少了6.0%，用EnKF方法同化有效波高是可行的。

Making Connections.

The piece is written for an ensemble of alto/bass flute, alto saxophone, viola, vibraphone/bongos and piano. Written as part of a telematic research project connecting musicians and dancers in remote sites, the work involved a collaboration with musicians at Edinburgh Napier University and dancers at Liverpool John Moores University and Nova Southeastern University, Florida.