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]

Embracing Integrability in Stellar and Planetary Dynamics

Hamiltonian systems in stellar and planetary dynamics are typically near integrable. For example, Solar System planets are almost in two-body orbits, and in simulations of the Galaxy, the orbits of stars seem regular. For such systems, sophisticated numerical methods can be developed through integrable approximations. Following this theme, we discuss three distinct problems. We start by considering numerical integration techniques for planetary systems. Perturbation methods (that utilize the integrability of the two-body motion) are preferred over conventional "blind" integration schemes. We introduce perturbation methods formulated with Cartesian variables. In our numerical comparisons, these are superior to their conventional counterparts, but, by definition, lack the energy-preserving properties of symplectic integrators. However, they are exceptionally well suited for relatively short-term integrations in which moderately high positional accuracy is required. The next exercise falls into the category of stability questions in solar systems. Traditionally, the interest has been on the orbital stability of planets, which have been quantified, e.g., by Liapunov exponents. We offer a complementary aspect by considering the protective effect that massive gas giants, like Jupiter, can offer to Earth-like planets inside the habitable zone of a planetary system. Our method produces a single quantity, called the escape rate, which characterizes the system of giant planets. We obtain some interesting results by computing escape rates for the Solar System. Galaxy modelling is our third and final topic. Because of the sheer number of stars (about 10^11 in Milky Way) galaxies are often modelled as smooth potentials hosting distributions of stars. Unfortunately, only a handful of suitable potentials are integrable (harmonic oscillator, isochrone and Stäckel potential). This severely limits the possibilities of finding an integrable approximation for an observed galaxy. A solution to this problem is torus construction; a method for numerically creating a foliation of invariant phase-space tori corresponding to a given target Hamiltonian. Canonically, the invariant tori are constructed by deforming the tori of some existing integrable toy Hamiltonian. Our contribution is to demonstrate how this can be accomplished by using a Stäckel toy Hamiltonian in ellipsoidal coordinates.

Origin of a trend in ECMWF wind analysis and its objective removal

The suitability of the European Centre for Medium Range Weather Forecasting (ECMWF) operational wind analysis for the period 1980-1991 for studying interannual variability is examined. The changes in the model and the analysis procedure are shown to give rise to a systematic and significant trend in the large scale circulation features. A new method of removing the systematic errors at all levels is presented using multivariate EOF analysis. Objectively detrended analysis of the three-dimensional wind field agrees well with independent Florida State University (FSU) wind analysis at the surface. It is shown that the interannual variations in the detrended surface analysis agree well in amplitude as well as spatial patterns with those of the FSU analysis. Therefore, the detrended analyses at other levels as well are expected to be useful for studies of variability and predictability at interannual time scales. It is demonstrated that this trend in the wind field is due to the shift in the climatologies from the period 1980-1985 to the period 1986-1991.

Time scales for transitions between free-energy minima of a dense hard-sphere system

Time scales associated with activated transitions between glassy metastable states of a free-energy functional appropriate for a dense hard-sphere system are calculated by using a new Monte Carlo method for the local density variables. In particular, we calculate the time the system, initially placed in a shallow glassy minimum of the free-energy, spends in the neighborhood of this minimum before making a transition to the basin of attraction of another free-energy minimum. This time scale is found to increase as the average density is increased. We find a crossover density near which this time scale increases very sharply and becomes longer than the longest times accessible in our simulation. This time scale does not show any evidence of increasing with sample size

Resent radar observations of asteroid 1566 Icarus

We report Doppler-only radar observations of Icarus at Goldstone at a transmitter frequency of 8510 MHz (3.5 cm wavelength) during 8-10 June 1996, the first radar detection of the object since 1968. Optimally filtered and folded spectra achieve a maximum opposite-circular (OC) polarization signal-to-noise ratio of about 10 and help to constrain Icarus' physical properties. We obtain an OC radar cross section of 0.05 km(2) (with a 35% uncertainty), which is less than values estimated by Goldstein (1969) and by Pettengill et al. (1969), and a circular polarization (SC/OC) ratio of 0.5+/-0.2. We analyze the echo power spectrum with a model incorporating the echo bandwidth B and a spectral shape parameter it, yielding a coupled constraint between B and n. We adopt 25 Hz as the lower bound on B, which gives a lower bound on the maximum pole-on breadth of about 0.6 km and upper bounds on the radar and optical albedos that are consistent with Icarus' tentative QS classification. The observed circular polarization ratio indicates a very rough near-surface at spatial scales of the order of the radar wavelength. (C) 1999 Elsevier Science Ltd. All rights reserved.

Multiple scales without center manifold reductions for delay differential equations near Hopf bifurcations

We study small perturbations of three linear Delay Differential Equations (DDEs) close to Hopf bifurcation points. In analytical treatments of such equations, many authors recommend a center manifold reduction as a first step. We demonstrate that the method of multiple scales, on simply discarding the infinitely many exponentially decaying components of the complementary solutions obtained at each stage of the approximation, can bypass the explicit center manifold calculation. Analytical approximations obtained for the DDEs studied closely match numerical solutions.

Role of length scales on microwave thawing dynamics in 2D cylinders

Microwave (MW) thawing of 2D frozen cylinders exposed to uniform plane waves from one face, is modeled using the effective heat capacity formulation with the MW power obtained from the electric field equations. Computations are illustrated for tylose (23% methyl cellulose gel) which melts over a range of temperatures giving rise to a mushy zone. Within the mushy region the dielectric properties are functions of the liquid volume fraction. The resulting coupled, time dependent non-linear equations are solved using the Galerkin finite element method with a fixed mesh. Our method efficiently captures the multiple connected thawed domains that arise due to the penetration of MWs in the sample. For a cylinder of diameter D, the two length scales that control the thawing dynamics are D/D-p and D/lambda(m), where D-p and lambda(m) are the penetration depth and wavelength of radiation in the sample respectively. For D/D-p, D/lambda(m) much less than 1 power absorption is uniform and thawing occurs almost simultaneously across the sample (Regime I). For D/D-p much greater than 1 thawing is seen to occur from the incident face, since the power decays exponentially into the sample (Regime III). At intermediate values, 0.2 < D/D-p, D/lambda(m) < 2.0 (Regime II) thawing occurs from the unexposed face at smaller diameters, from both faces at intermediate diameters and from the exposed and central regions at larger diameters. Average power absorption during thawing indicates a monotonic rise in Regime I and a monotonic decrease in Regime III. Local maxima in the average power observed for samples in Regime II are due to internal resonances within the sample. Thawing time increases monotonically with sample diameter and temperature gradients in the sample generally increase from Regime I to Regime III. (C) 2002 Elsevier Science Ltd. All rights reserved.

Electronic structure of ${\rm Ca}_{1-x}{\rm Sr}_x{\rm VO}_3$: A tale of two energy scales

We investigate the electronic structure of Ca1-xSrxVO3 using photoemission spectroscopy. Core level spectra establish an electronic phase separation at the surface, leading to a distinctly different surface electronic structure compared to the bulk. Analysis of the photoemission spectra of this system allowed us to separate the surface and bulk contributions. These results help us to understand properties related to two vastly differing energy scales, namely the low-energy scale of thermal excitations ( $\sim\!k_{\rm B}T$) and the high-energy scale related to Coulomb and other electronic interactions.

On breaks of the Indian monsoon

For over a century, the term break has been used for spells in which the rainfall over the Indian monsoon zone is interrupted. The phenomenon of 'break monsoon' is of great interest because long intense breaks are often associated with poor monsoon seasons. Such breaks have distinct circulation characteristics (heat trough type circulation) and have a large impact on rainfed agriculture. Although interruption of the monsoon rainfall is considered to be the most important feature of the break monsoon, traditionally breaks have been identified on the basis of the surface pressure and wind patterns over the Indian region. We have defined breaks (and active spells) on the basis of rainfall over the monsoon zone. The rainfall criteria are chosen so as to ensure a large overlap with the traditional breaks documented by Ramamurthy (1969) and De et al (1998). We have identified these rainbreaks for 1901-89. We have also identified active spells on the basis of rainfall over the Indian monsoon zone. We have shown that the all-India summer monsoon rainfall is significantly negatively correlated with the number of rainbreak days (correlation coefficient -0.56) and significantly positively correlated with the number of active days (correlation coefficient 0.47). Thus the interannual variation of the all-India summer monsoon rainfall is shown to be related to the number of days of rainbreaks and active spells identified here. There have been several studies of breaks (and also active spells in several cases) identified on the basis of different criteria over regions differing in spatial scales (e.g., Webster et al 1998; Krishnan et al 2000; Goswami and Mohan 2000; and Annamalai and Slingo 2001). We find that there is considerable overlap between the rainbreaks we have identified and breaks based on the traditional definition. There is some overlap with the breaks identified by Krishnan et al (2000) but little overlap with breaks identified by Webster et al (1998). Further, there are three or four active-break cycles in a season according to Webster et al (1998) which implies a time scale of about 40 days for which Goswami and Mohan (2000), and Annamalai and Slingo'(2001) have studied breaks and active minus break fluctuations. On the other hand, neither the traditional breaks (Ramamurthy 1969; and De et al 1998) nor the rainbreaks occur every year. This suggests that the 'breaks' in these studies axe weak spells of the intraseasonal variation of the monsoon, which occur every year. We have derived the OLR and circulation patterns associated with rainbreaks and active spells and compared them with the patterns associated with breaks/active minus break spells from these studies. Inspite of differences in the patterns over the Indian region, there is one feature which is seen in the OLR anomaly patterns of breaks identified on the basis of different criteria as well as the rainbreaks identified in this paper viz., a quadrapole over the Asia-west Pacific region arising from anomalies opposite (same) in sign to those over the Indian region occurring over the equatorial Indian Ocean and northern tropical (equatorial) parts of the west Pacific. Thus it appears that this quadrapole is a basic feature of weak spells of the intraseasonal,variation over the Asia-west Pacific region. Since the rainbreaks are intense weak spells, this basic feature is also seen in the composite patterns of these breaks. We find that rainbreaks (active spells) are also associated with negative (positive) anomalies over a part of the cast Pacific suggesting that the convection over the Indian region is linked to that over the east Pacific not only on the interannual scale (as evinced by the link between the Indian summer monsoon rainfall and ENSO) but on the intraseasonal scale as well.

Influence of Spatial Scales on the Performance of Saturated Subsurface Flow and Transport Models

Numerical modeling of saturated subsurface flow and transport has been widely used in the past using different numerical schemes such as finite difference and finite element methods. Such modeling often involves discretization of the problem in spatial and temporal scales. The choice of the spatial and temporal scales for a modeling scenario is often not straightforward. For example, a basin-scale saturated flow and transport analysis demands larger spatial and temporal scales than a meso-scale study, which in turn has larger scales compared to a pore-scale study. The choice of spatial-scale is often dictated by the computational capabilities of the modeler as well as the availability of fine-scale data. In this study, we analyze the impact of different spatial scales and scaling procedures on saturated subsurface flow and transport simulations.

On the linkage of mesospheric planetary waves with those of the lower atmosphere and ionosphere: A case study from Indian low latitudes

In this paper we study the planetary-scale wave features using concurrent observations of mesospheric wind and temperature, ionospheric h'F, and tropospheric wind from Tirunelveli, Gadanki, and Kolhapur, all located in the Indian low latitudes, made during February 2009. Our investigations reveal that 3 to 5 day periodicity, characterized as ultrafast Kelvin (UFK) waves, was persistent throughout the atmosphere during this period. These waves show clear signatures of upward wave propagation from troposphere to the upper mesosphere, linking the ionosphere through a clear correlation between mesospheric winds and h'F variations. We also note that the amplitude of this wave decreased as we moved away from the equator. These results are the first of their kind from Indian sector, portraying the vertical as well as latitudinal characteristics of the 3 to 5 day UFK waves simultaneously from the troposphere to the ionosphere.

Postcards from oases in the desert: phenomenology of SUSY with intermediate scales

The presence of new matter fields charged under the Standard Model gauge group at intermediate scales below the Grand Unification scale modifies the renormalization group evolution of the gauge couplings. This can in turn significantly change the running of the Minimal Supersymmetric Standard Model parameters, in particular the gaugino and the scalar masses. In the absence of new large Yukawa couplings we can parameterise all the intermediate scale models in terms of only two parameters controlling the size of the unified gauge coupling. As a consequence of the modified running, the low energy spectrum can be strongly affected with interesting phenomenological consequences. In particular, we show that scalar over gaugino mass ratios tend to increase and the regions of the parameter space with neutralino Dark Matter compatible with cosmological observations get drastically modified. Moreover, we discuss some observables that can be used to test the intermediate scale physics at the LHC in a wide class of models.

Space-Time Scales of Northward Propagation of Convection during Boreal Summer

This study uses precipitation estimates from the Tropical Rainfall Measuring Mission to quantify the spatial and temporal scales of northward propagation of convection over the Indian monsoon region during boreal summer. Propagating modes of convective systems in the intraseasonal time scales such as the Madden-Julian oscillation can interact with the intertropical convergence zone and bring active and break spells of the Indian summer monsoon. Wavelet analysis was used to quantify the spatial extent (scale) and center of these propagating convective bands, as well as the time period associated with different spatial scales. Results presented here suggest that during a good monsoon year the spatial scale of this oscillation is about 30 degrees centered around 10 degrees N. During weak monsoon years, the scale of propagation decreases and the center shifts farther south closer to the equator. A strong linear relationship is obtained between the center/scale of convective wave bands and intensity of monsoon precipitation over Indian land on the interannual time scale. Moreover, the spatial scale and its center during the break monsoon were found to be similar to an overall weak monsoon year. Based on this analysis, a new index is proposed to quantify the spatial scales associated with propagating convective bands. This automated wavelet-based technique developed here can be used to study meridional propagation of convection in a large volume of datasets from observations and model simulations. The information so obtained can be related to the interannual and intraseasonal variation of Indian monsoon precipitation.

On the mechanism and the length scales involved in the ductile fracture of a bulk metallic glass

Some bulk metallic glasses (BMGs) exhibit high crack initiation toughness due to shear band mediated plastic flow at the crack tip and yet do not display additional resistance to crack growth due to the lack of a microstructure. Thus, at crack initiation, the fracture behavior of BMGs transits from that of ductile alloys to that of brittle ceramics. In this paper, we attempt to understand the physics behind the characteristic length from the notch root at which this transition occurs, through testing of four-point bend specimens made of a nominally ductile Zr-based BMG in three different structural states. In the as-cast state, both symmetric (mode I) and asymmetric (mixed mode) bend specimens are tested. The process of shear band mediated plastic flow followed by crack initiation at the notch root was monitored through in situ imaging. Results show that stable crack growth occurs inside a dominant shear band through a distance of, similar to 60 mu m, irrespective of the structural state and mode mixity, before attaining criticality. Detailed finite element simulations show that this length corresponds to the distance from the notch root over which a positive hydrostatic stress gradient prevails. The mean ridge heights on fractured surfaces are found to correlate with the toughness of the BMG. The Argon and Salama model, which is based on the meniscus instability phenomenon at the notch root, is modified to explain the experimentally observed physics of fracture in ductile BMGs. (C) 2012 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.