Rigidity of Orientationally Ordered Domains of Short Chain Molecules

By molecular dynamics simulation, discovered is a strange rigid-like nature for a hexagonally packed domain of short chain molecules. In spite of the non-bonded short-range interaction potential (Lennard-Jones potential) among chain molecules, the packed domain gives rise to a resultant global moment of inertia. Accordingly, as two domains encounter obliquely, they rotate so as to be parallel to each other keeping their overall structures as if they were rigid bodies.Comment: 7 pages, 5 figures, and 2 table

Large Area Crop Inventory Experiment (LACIE). Intensive test site assessment report

There are no author-identified significant results in this report

Landscape Genetics of Raccoons (\u3ci\u3eProcyon lotor\u3c/i\u3e) Associated with Ridges and Valleys of Pennsylvania: Implications for Oral Rabies Vaccination Programs

Raccoons are the reservoir for the raccoon rabies virus variant in the United States. To combat this threat, oral rabies vaccination (ORV) programs are conducted in many eastern states. To aid in these efforts, the genetic structure of raccoons (Procyon lotor) was assessed in southwestern Pennsylvania to determine if select geographic features (i.e., ridges and valleys) serve as corridors or hindrances to raccoon gene flow (e.g., movement) and, therefore, rabies virus trafficking in this physiographic region. Raccoon DNA samples (n = 185) were collected from one ridge site and two adjacent valleys in southwestern Pennsylvania (Westmoreland, Cambria, Fayette, and Somerset counties). Raccoon genetic structure within and among these study sites was characterized at nine microsatellite loci. Results indicated that there was little population subdivision among any sites sampled. Furthermore, analyses using a model-based clustering approach indicated one essentially panmictic population was present among all the raccoons sampled over a reasonably broad geographic area (e.g., sites up to 36 km apart). However, a signature of isolation by distance was detected, suggesting that widths of ORV zones are critical for success. Combined, these data indicate that geographic features within this landscape influence raccoon gene flow only to a limited extent, suggesting that ridges of this physiographic system will not provide substantial long-term natural barriers to rabies virus trafficking. These results may be of value for future ORV efforts in Pennsylvania and other eastern states with similar landscapes

Algorithm for numerical integration of the rigid-body equations of motion

A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in terms of either principal axes or quaternions. Due to specific features of the algorithm, orthonormality and unit norms of the orientational variables are integrals of motion, despite an approximate character of the produced trajectories. It is shown that the method presented appears to be the most efficient among all known algorithms of such a kind.Comment: 4 pages, 1 figur

On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems

A consequent approach is proposed to construct symplectic force-gradient algorithms of arbitrarily high orders in the time step for precise integration of motion in classical and quantum mechanics simulations. Within this approach the basic algorithms are first derived up to the eighth order by direct decompositions of exponential propagators and further collected using an advanced composition scheme to obtain the algorithms of higher orders. Contrary to the scheme by Chin and Kidwell [Phys. Rev. E 62, 8746 (2000)], where high-order algorithms are introduced by standard iterations of a force-gradient integrator of order four, the present method allows to reduce the total number of expensive force and its gradient evaluations to a minimum. At the same time, the precision of the integration increases significantly, especially with increasing the order of the generated schemes. The algorithms are tested in molecular dynamics and celestial mechanics simulations. It is shown, in particular, that the efficiency of the new fourth-order-based algorithms is better approximately in factors 5 to 1000 for orders 4 to 12, respectively. The results corresponding to sixth- and eighth-order-based composition schemes are also presented up to the sixteenth order. For orders 14 and 16, such highly precise schemes, at considerably smaller computational costs, allow to reduce unphysical deviations in the total energy up in 100 000 times with respect to those of the standard fourth-order-based iteration approach.Comment: 23 pages, 2 figures; submitted to Phys. Rev.

Dynamical Properties and Plasmon Dispersion of a Weakly Degenerate Correlated One-Component Plasma

Classical Molecular Dynamics (MD) simulations for a one-component plasma (OCP) are presented. Quantum effects are included in the form of the Kelbg potential. Results for the dynamical structure factor are compared with the Vlasov and RPA (random phase approximation) theories. The influence of the coupling parameter $\Gamma$, degeneracy parameter $\rho \Lambda^3$ and the form of the pair interaction on the optical plasmon dispersion is investigated. An improved analytical approximation for the dispersion of Langmuir waves is presented.Comment: 23 pages, includes 7 ps/eps-figures and 2 table

Discovery of Bright Galactic R Coronae Borealis and DY Persei Variables: Rare Gems Mined from ACVS

We present the results of a machine-learning (ML) based search for new R Coronae Borealis (RCB) stars and DY Persei-like stars (DYPers) in the Galaxy using cataloged light curves from the All-Sky Automated Survey (ASAS) Catalog of Variable Stars (ACVS). RCB stars - a rare class of hydrogen-deficient carbon-rich supergiants - are of great interest owing to the insights they can provide on the late stages of stellar evolution. DYPers are possibly the low-temperature, low-luminosity analogs to the RCB phenomenon, though additional examples are needed to fully establish this connection. While RCB stars and DYPers are traditionally identified by epochs of extreme dimming that occur without regularity, the ML search framework more fully captures the richness and diversity of their photometric behavior. We demonstrate that our ML method can use newly discovered RCB stars to identify additional candidates within the same data set. Our search yields 15 candidates that we consider likely RCB stars/DYPers: new spectroscopic observations confirm that four of these candidates are RCB stars and four are DYPers. Our discovery of four new DYPers increases the number of known Galactic DYPers from two to six; noteworthy is that one of the new DYPers has a measured parallax and is m ~ 7 mag, making it the brightest known DYPer to date. Future observations of these new DYPers should prove instrumental in establishing the RCB connection. We consider these results, derived from a machine-learned probabilistic classification catalog, as an important proof-of-concept for the efficient discovery of rare sources with time-domain surveys.Comment: 18 pages, 2 new figures, accepted for publication in Ap

Atomic-scale modeling of the deformation of nanocrystalline metals

Nanocrystalline metals, i.e. metals with grain sizes from 5 to 50 nm, display technologically interesting properties, such as dramatically increased hardness, increasing with decreasing grain size. Due to the small grain size, direct atomic-scale simulations of plastic deformation of these materials are possible, as such a polycrystalline system can be modeled with the computational resources available today. We present molecular dynamics simulations of nanocrystalline copper with grain sizes up to 13 nm. Two different deformation mechanisms are active, one is deformation through the motion of dislocations, the other is sliding in the grain boundaries. At the grain sizes studied here the latter dominates, leading to a softening as the grain size is reduced. This implies that there is an ``optimal'' grain size, where the hardness is maximal. Since the grain boundaries participate actively in the deformation, it is interesting to study the effects of introducing impurity atoms in the grain boundaries. We study how silver atoms in the grain boundaries influence the mechanical properties of nanocrystalline copper.Comment: 10 pages, LaTeX2e, PS figures and sty files included. To appear in Mater. Res. Soc. Symp. Proc. vol 538 (invited paper). For related papers, see http://www.fysik.dtu.dk/~schiotz/publist.htm

Structure and relaxations in liquid and amorphous Selenium

We report a molecular dynamics simulation of selenium, described by a three-body interaction. The temperatures T_g and T_c and the structural properties are in agreement with experiment. The mean nearest neighbor coordination number is 2.1. A small pre-peak at about 1 AA^-1 can be explained in terms of void correlations. In the intermediate self-scattering function, i.e. the density fluctuation correlation, classical behavior, alpha- and beta-regimes, is found. We also observe the plateau in the beta-regime below T_g. In a second step, we investigated the heterogeneous and/or homogeneous behavior of the relaxations. At both short and long times the relaxations are homogeneous (or weakly heterogeneous). In the intermediate time scale, lowering the temperature increases the heterogeneity. We connect these different domains to the vibrational (ballistic), beta- and alpha-regimes. We have also shown that the increase in heterogeneity can be understood in terms of relaxations