Orbital evolution studies of two HMXB pulsars

High mass X-ray binary (HMXB) pulsars are of two types, persistent and transient. 4U1538-52 is a persistent HMXB whose orbit was previously measured to be circular but the RXTE observations revealed an eccentric orbit. We observed this system with RXTE-PCA in August 2003 and our timing analysis supports the eccentric orbit of the system. However, we do not find any evidence for orbital evolution. Rotational and tidal interactions between the stars of a closed binary system result in apsidal motion which can be measured in systems with eccentric orbit. 4U0115+63 is a Be-transient HMXB whose eccentric orbit was well-determined during its 1978 outburst. We report preliminary results from analysis of data obtained during the 1999 outburst of this source with the RXTE-PCA.

Orbital evolution studies of two HMXB pulsars

High mass X-ray binary (H M X B) pulsars are of two types, persistent and transient. 4U 1538-52 is a persistent HMXB whose orbit was previously measured to be circular but the RXTE observations revealed an eccentric orbit. We observed this system with RXTE-PCA in August 2003 and our timing analysis supports the eccentric orbit of the system. However, we do not find any evidence for orbital evolution. Rotational and tidal interactions between the stars of a closed binary system result in apsidal motion which can be measured in systems with eccentric orbit. 4U0115+63 is a Be-transient HMXB whose eccentric orbit was well-determined during its 1978 outburst. We report preliminary results from analysis of data obtained during the 1999 outburst of this source with the RXTE-PCA.

Time Scale Separated Nonlinear Partial Integrated Guidance and Control of Interceptors in the Terminal Phase

The paper proposes a time scale separated partial integrated guidance and control of an interceptor for engaging high speed targets in the terminal phase. In this two loop design, the outer loop is an optimal control formulation based on nonlinear model predictive spread control philosophies. It gives the commanded pitch and yaw rates whereas necessary roll-rate command is generated from a roll-stabilization loop. The inner loop tracks the outer loop commands using the dynamicinversion philosophy. However, unlike conventional designs, in both the loops the Six degree of freedom (Six-DOF) interceptor model is used directly. This intelligent manipulation preserves the inherent time scale separation property between the translational and rotational dynamics, and hence overcomes the deficiency of current IGC designs, while preserving its benefits. Six-DOF simulation studies have been carried out accounting for three dimensional engagement geometry. Different comparison studies were also conducted to measure the performance of the algorithm.

Structure and dynamics of benzene in one-dimensional channels

Molecular dynamics investigation of benzene in one-dimensional channel systems A1PO(4)-5, VPI-5, and carbon nanotube is reported. The results suggest that, in all the three host systems, the plane of benzene is almost perpendicular to the channel axis when the molecule is near the center of the channel and the plane of benzene is parallel to the channel axis when the molecule is near the wall of the channel. The density distribution of benzene as a function of channel length, z and the radial distance, r, from the channel axis is also different in the three host structures. Anisotropy in translational diffusion coefficient, calculated in body-fixed frame of benzene, suggests that benzene prefers to move with its plane parallel to the direction of motion in A1PO(4)-5 and VPI-5 whereas in carbon nanotube the motion occurs predominantly with the plane of the benzene perpendicular to the direction of motion.;Anisotropy associated with the rotational motion is seen to alter significantly in confinement as compared to liquid benzene. In A1PO(4)-5, the rotational anisotropy is reversed as compared to liquid benzene thereby suggesting that anisotropy arising out of molecular geometry can be reduced. Reorientational correlation times for C-6 and C-2 axes Of benzene are reported. Apart from the inertial decay of reorientational correlation function due to free, rotation, two other distinct regimes of decay are observed in narrower channels (AIPO(4)-5 and carbon nanotube): (i) an initial fast decay (0.5-2 ps) and (ii) a slower decay (>2 ps) of reorientational correlation function where C-6 decays slower than C-2 Similar to what is observed in liquid benzene. In the initial fast decay, it is seen that the decay for C-6 is faster than C-2 which is in contrast to what is observed in liquid benzene or for benzene confined in VPI-5.

Determination of distances of sugar protons from Mn2+ in concanavalin A

1H NMR spin-lattice relaxation time (T1) measurements have been carried out with various sugars, viz. methyl alpha-D-glucopyranoside (alpha-MeGluP), methyl beta-D-lucopyranoside (beta-MeGluP), methyl alpha--annopyranoside (alpha-MeManP), maltose (4-O-alpha-D-glucopyranosyl--glucose), nigerose (3-O-alpha-D-glucopyranosyl-D-glucose), p-nitrophenyl alpha-maltoside (PNP-alpha-maltoside) and p-nitrophenyl beta-maltoside (PNP-beta-maltoside) to determine the distances of sugar protons from Mn2+ in concanavalin A (Con A). With a rotational correlation time of 1.58 x 10(-10) s determined, distances were calculated using Solomon-Bloembergen equation. The data obtained indicated differences in disposition of different groups in the binding site of Con A. An average value of about 10 A was obtained for the distances of sugar protons from Mn2+ in Con A. In the case of mono and disaccharides, the non-reducing end sugar unit was found to be closer to Mn2+ than the reducing end one.

Polarization relaxation, dielectric dispersion, and solvation dynamics in dense dipolar liquid

A unified treatment of polarization relaxation, dielectric dispersion and solvation dynamics in a dense, dipolar liquid is presented. It is shown that the information of solvent polarization relaxation that is obtained by macroscopic dielectric dispersion experiments is not sufficient to understand dynamics of solvation of a newly created ion or dipole. In solvation, a significant contribution comes from intermediate wave vector processes which depend critically on the short range (nearest‐neighbor) spatial and orientational order that are present in a dense, dipolar liquid. An analytic expression is obtained for the time dependent solvation energy that depends, in addition to the translational and rotational diffusion coefficients of the liquid, on the ratio of solute–solvent molecular sizes and on the microscopic structure of the polar liquid. Mean spherical approximation (MSA) theory is used to obtain numerical results for polarization relaxation, for wave vector and frequency dependent dielectric function and for time dependent solvation energy. We find that in the absence of translational contribution, the solvation of an ion is, in general, nonexponential. In this case, the short time decay is dominated by the longitudinal relaxation time but the long time decay is dominated by much slower large wave vector processes involving nearest‐neighbor molecules. The presence of a significant translational contribution drastically alters the decay behavior. Now, the long‐time behavior is given by the longitudinal relaxation time constant and the short time dynamics is controlled by the large wave vector processes. Thus, although the continuum model itself is conceptually wrong, a continuum model like result is recovered in the presence of a sizeable translational contribution. The continuum model result is also recovered in the limit of large solute to solvent size ratio. In the opposite limit of small solute size, the decay is markedly nonexponential (if the translational contribution is not very large) and a complete breakdown of the continuum model takes place. The significance of these results is discussed.

NMR study of ammonium halides under high pressure

Wide-line c.w. proton resonance investigations have been carried out on the ammonium halides, namely, ammonium chloride, ammonium bromide and ammonium iodide in the temperature range between 77 and 300 K and in the pressure range between 1 bar and 14 kbar. It has been found that the narrow iodide spectrum at 77 K broadens under the application of hydrostatic pressure. The barrier height for the ammonium ion motion in ammonium iodide under pressure has been estimated by carrying out a temperature variation study. The rotational potential for the motion of ammonium ion in ammonium iodide at 1 bar and 14 kbar has been calculated using earlier theoretical models and compared with values calculated for ammonium chloride and bromide. The barrier height in the case of ammonium iodide under pressure is found to be of the same order of magnitude as the value obtained in the case of ammonium bromide at atmospheric pressure indicating that the high pressure phase of ammonium iodide is likely to have the same structure as the low temperature ordered CsCl phase found in the case of the chloride and the bromide. The increase in the potential barrier height in the case of ammonium iodide under pressure indicates that the reorientational motion executed by the ammonium ions is inhibited by the application of pressure. This is also confirmed by the broadening of the spectral line at 77 K under the application of pressure.

Spin and mass content of linearized Poincaré gauge theories

A unified gauge theory of massless and massive spin-2 fields is of considerable current interest. The Poincaré gauge theories with quadratic Lagrangian are linearized, and the conditions on the parameters are found which will lead to viable linear theories with massive gauge particles. As well as the 2+ massless gravitons coming from the translational gauge potential, the rotational gauge potentials, in the linearized limit, give rise to 2+ and 2− particles of equal mass, as well as a massive pseudoscalar.

A kinetic Ising model study of dynamical correlations in confined fluids: Emergence of both fast and slow time scales

Experiments and computer simulation studies have revealed existence of rich dynamics in the orientational relaxation of molecules in confined systems such as water in reverse micelles, cyclodextrin cavities, and nanotubes. Here we introduce a novel finite length one dimensional Ising model to investigate the propagation and the annihilation of dynamical correlations in finite systems and to understand the intriguing shortening of the orientational relaxation time that has been reported for small sized reverse micelles. In our finite sized model, the two spins at the two end cells are oriented in the opposite directions to mimic the effects of surface that in real system fixes water orientation in the opposite directions. This produces opposite polarizations to propagate inside from the surface and to produce bulklike condition at the center. This model can be solved analytically for short chains. For long chains, we solve the model numerically with Glauber spin flip dynamics (and also with Metropolis single-spin flip Monte Carlo algorithm). We show that model nicely reproduces many of the features observed in experiments. Due to the destructive interference among correlations that propagate from the surface to the core, one of the rotational relaxation time components decays faster than the bulk. In general, the relaxation of spins is nonexponential due to the interplay between various interactions. In the limit of strong coupling between the spins or in the limit of low temperature, the nature of relaxation of the spins undergoes a qualitative change with the emergence of a homogeneous dynamics where decay is predominantly exponential, again in agreement with experiments. (C) 2010 American Institute of Physics. doi: 10.1063/1.3474948]

An accurate numerical integration scheme for finite rotations using rotation vector parametrization

A numerical integration procedure for rotational motion using a rotation vector parametrization is explored from an engineering perspective by using rudimentary vector analysis. The incremental rotation vector, angular velocity and acceleration correspond to different tangent spaces of the rotation manifold at different times and have a non-vectorial character. We rewrite the equation of motion in terms of vectors lying in the same tangent space, facilitating vector space operations consistent with the underlying geometric structure. While any integration algorithm (that works within a vector space setting) may be used, we presently employ a family of explicit Runge-Kutta algorithms to solve this equation. While this work is primarily motivated out of a need for highly accurate numerical solutions of dissipative rotational systems of engineering interest, we also compare the numerical performance of the present scheme with some of the invariant preserving schemes, namely ALGO-C1, STW, LIEMIDEA] and SUBCYC-M. Numerical results show better local accuracy via the present approach vis-a-vis the preserving algorithms. It is also noted that the preserving algorithms do not simultaneously preserve all constants of motion. We incorporate adaptive time-stepping within the present scheme and this in turn enables still higher accuracy and a `near preservation' of constants of motion over significantly longer intervals. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Theoretical gain‐optimization studies in CO2–N2 gasdynamic lasers. II. Results of parametric study

Based on a method presented in detail in a previous work by the authors, similar solutions have been obtained for the steady inviscid quasi‐one‐dimensional nonreacting flow in the supersonic nozzle of a CO2–N2 gasdynamic laser system, with either H2O or He as the catalyst. It has been demonstrated how these solutions could be used to optimize the small‐signal gain coefficient on a specified vibrational‐rotational transition. Results presented for a wide range of mixture compositions include optimum values for the small‐signal gain, area ratio, reservoir temperature, and a binary scaling parameter, which is the product of reservoir pressure and nozzle shape factor.

Dynamics of sheared inelastic dumbbells

We study the dynamical properties of the homogeneous shear flow of inelastic dumbbells in two dimensions as a first step towards examining the effect of shape on the properties of flowing granular materials. The dumbbells are modelled as smooth fused disks characterized by the ratio of the distance between centres (L) and the disk diameter (D), with an aspect ratio (L/D) varying between 0 and 1 in our simulations. Area fractions studied are in the range 0.1-0.7, while coefficients of normal restitution (e(n)) from 0.99 to 0.7 are considered. The simulations use a modified form of the event-driven methodology for circular disks. The average orientation is characterized by an order parameter S, which varies between 0 (for a perfectly disordered fluid) and 1 (for a fluid with the axes of all dumbbells in the same direction). We investigate power-law fits of S as a function of (L D) and (1 - e(n)(2)) There is a gradual increase in ordering as the area fraction is increased, as the aspect ratio is increased or as the coefficient of restitution is decreased. The order parameter has a maximum value of about 0.5 for the highest area fraction and lowest coefficient of restitution considered here. The mean energy of the velocity fluctuations in the flow direction is higher than that in the gradient direction and the rotational energy, though the difference decreases as the area fraction increases, due to the efficient collisional transfer of energy between the three directions. The distributions of the translational and rotational velocities are Gaussian to a very good approximation. The pressure is found to be remarkably independent of the coefficient of restitution. The pressure and dissipation rate show relatively little variation when scaled by the collision frequency for all the area fractions studied here, indicating that the collision frequency determines the momentum transport and energy dissipation, even at the lowest area fractions studied here. The mean angular velocity of the particles is equal to half the vorticity at low area fractions, but the magnitude systematically decreases to less than half the vorticity as the area fraction is increased, even though the stress tensor is symmetric.

Molecular theory of dielectric relaxation in a dense binary dipolar liquid

A molecular theory of dielectric relaxation in a dense binary dipolar liquid is presented. The theory takes into account the effects of intra- and interspecies intermolecular interactions. It is shown that the relaxation is, in general, nonexponential. In certain limits, we recover the biexponential form traditionally used to analyze the experimental data of dielectric relaxation in a binary mixture. However, the relaxation times are widely different from the prediction of the noninteracting rotational diffusion model of Debye for a binary system. Detailed numerical evaluation of the frequency-dependent dielectric function epsilon-(omega) is carried out by using the known analytic solution of the mean spherical approximation (MSA) model for the two-particle direct correlation function for a polar mixture. A microscopic expression for both wave vector (k) and frequency (omega) dependent dielectric function, epsilon-(k,omega), of a binary mixture is also presented. The theoretical predictions on epsilon-(omega) (= epsilon-(k = 0, omega)) have been compared with the available experimental results. In particular, the present theory offers a molecular explanation of the phenomenon of fusing of the two relaxation channels of the neat liquids, observed by Schallamach many years ago.

Rocking response of rectangular rigid blocks under random noise base excitations

The response of a rigid rectangular block resting on a rigid foundation and acted upon simultaneously by a horizontal and a vertical random white-noise excitation is considered. In the equation of motion, the energy dissipation is modeled through a viscous damping term. Under the assumption that the body does not topple, the steady-state joint probability density function of the rotation and the rotational velocity is obtained using the Fokker-Planck equation approach. Closed form solution is obtained for a specific combination of system parameters. A more general but approximate solution to the joint probability density function based on the method of equivalent non-linearization is also presented. Further, the problem of overturning of the block is approached in the framework of the diffusion methods for first passage failure studies. The overturning of the block is deemed incipient when the response trajectories in the phase plane cross the separatrix of the conservative unforced system. Expressions for the moments of first passage time are obtained via a series solution to the governing generalized Pontriagin-Vitt equations. Numerical results illustra- tive of the theoretical solutions are presented and their validity is examined through limited amount of digital simulations.

Temperature and concentration dependence of adsorption properties of methane in NaY: a molecular dynamics study

Molecular dynamics calculations on methane sorbed in NaY (Si/Al = 3.0) employing realistic methane-methane and methane-zeolite intermolecular potential functions at different temperatures (50, 150, 220, and 300 K) and concentrations (2, 4, 6, and 8 molecules/cage) are reported. The thermodynamic results are in agreement with the available experimental data. Guest-guest and guest-host radial distribution functions (rdfs), energy distribution functions, distribution of cage occupancy, center-of-cage-center-of-mass (coc-com) rdfs, velocity autocorrelation functions for com and angular motion and the Fourier transformed power spectra, and diffusion coefficients are presented as a function of temperature and concentration. At 50 K, methane is localized near the adsorption site. Site-site migration and essentially free rotational motion are observed at 150 K. Molecules preferentially occupy the region near the inner surface of the alpha-cage. The vibrational frequencies for the com of methane shift toward higher values with decreasing temperature and increasing adsorbate concentration. The observed frequencies for com motion are 36, 53, and 85 cm-1 and for rotational motion at 50 K, 95 and 150 cm-1 in agreement with neutron scattering data. The diffusion coefficients show a type I behavior as a function of loading in agreement with NMR measurements. Cage-to-cage diffusion is found to be always mediated by the surface.