Molecular Dynamics of Thermotropic Liquid Crystals: Anomalous relaxation dynamics of calamitic and discotic liquid crystals

Recent optical kerr effect (OKE) studies have demonstrated that orientational relaxation of rod-like nematogens exhibits temporal power law decay at intermediate times not only near the isotropic–nematic (I–N) phase boundary but also in the nematic phase. Such behaviour has drawn an intriguing analogy with supercooled liquids. We have investigated both collective and single-particle orientational dynamics of a family of model system of thermotropic liquid crystals using extensive computer simulations. Several remarkable features of glassy dynamics are on display including non-exponential relaxation, dynamical heterogeneity, and non-Arrhenius temperature dependence of the orientational relaxation time. Over a temperature range near the I–N phase boundary, the system behaves remarkably like a fragile glass-forming liquid. Using proper scaling, we construct the usual relaxation time versus inverse temperature plot and explicitly demonstrate that one can successfully define a density dependent fragility of liquid crystals. The fragility of liquid crystals shows a temperature and density dependence which is remarkably similar to the fragility of glass forming supercooled liquids. Energy landscape analysis of inherent structures shows that the breakdown of the Arrhenius temperature dependence of relaxation rate occurs at a temperature that marks the onset of the growth of the depth of the potential energy minima explored by the system. A model liquid crystal, consisting of disk-like molecules, has also been investigated in molecular dynamics simulations for orientational relaxation along two isobars starting from the high temperature isotropic phase. The isobars have been so chosen that the phase sequence isotropic (I)–nematic (N)–columnar (C) appears upon cooling along one of them and the sequence isotropic (I)–columnar(C) along the other. While the orientational relaxation in the isotropic phase near the I–N phase transition shows a power law decay at short to intermediate times, such power law relaxation is not observed in the isotropic phase near the I–C phase boundary. The origin of the power law decay in the single-particle second-rank orientational time correlation function (OTCF) is traced to the growth of the orientational pair distribution functions near the I–N phase boundary. As the system settles into the nematic phase, the decay of the single-particle second-rank orientational OTCF follows a pattern that is similar to what is observed with calamitic liquid crystals and supercooled molecular liquids.

Wheeler-Feynman absorber revisited: a useful technique to calculate decay rates and lifetimes in small scale optical systems

The Wheeler-Feynman (WF) absorber theory of radiation though no more of interest in explaining self interaction of an electron, can be very useful in today's research in small scale optical systems. The significance of the WF absorber is the use of time-symmetrical solution of Maxwell's equations as opposed to only the retarded solution. The radiative coupling of emitters to nano wires in the near field and change in their lifetimes due to small mode volume enclosures have been elucidated with the retarded solutions before. These solutions have also been shown to agree with quantum electrodynamics, thus allowing for classical electromagnetic approaches in such problems. It is here assumed that the radiative coupling of the emitter with a body is in proportion to its contribution to the classical force of radiative reaction as derived in the WF absorber theory. Representing such nano structures as a partial WF absorber acting on the emitter makes the computations considerably easier than conventional electromagnetic solutions for full boundary conditions.

Probing the indefinite CP nature of the Higgs boson through decay distributions in the process e(+)e(-) -> t(t)over-bar Phi

The recently discovered scalar resonance at the Large Hadron Collider is now almost confirmed to be a Higgs boson, whose CP properties are yet to be established. At the International Linear Collider with and without polarized beams, it may be possible to probe these properties at high precision. In this work, we study the possibility of probing departures from the pure CP-even case, by using the decay distributions in the process e(+)e(-) -> t (t) over bar Phi, with Phi mainly decaying into a b (b) over bar pair. We have compared the case of a minimal extension of the Standard Model case (model I) with an additional pseudoscalar degree of freedom, with a more realistic case namely the CP-violating two-Higgs doublet model (model II) that permits a more general description of the couplings. We have considered the International Linear Collider with root s = 800 GeV and integrated luminosity of 300 fb(-1). Our main findings are that even in the case of small departures from the CP-even case, the decay distributions are sensitive to the presence of a CP-odd component in model II, while it is difficult to probe these departures in model I unless the pseudoscalar component is very large. Noting that the proposed degrees of beam polarization increase the statistics, the process demonstrates the effective role of beam polarization in studies beyond the Standard Model. Further, our study shows that an indefinite CP Higgs would be a sensitive laboratory to physics beyond the Standard Model.

Implications of Higgs boson to diphoton decay rate in the bilinear R-parity violating supersymmetric model

The Large Hadron Collider has recently discovered a Higgs-like particle having a mass around 125 GeVand also indicated that there is an enhancement in the Higgs to diphoton decay rate as compared to that in the standard model. We have studied implications of these discoveries in the bilinear R-parity violating supersymmetric model, whose main motivation is to explain the nonzero masses for neutrinos. The R-parity violating parameters in this model are epsilon and b(epsilon), and these parameters determine the scale of neutrino masses. If the enhancement in the Higgs to diphoton decay rate is true, then we have found epsilon greater than or similar to 0.01 GeV and b epsilon similar to 1 GeV2 in order to be compatible with the neutrino oscillation data. Also, in the above mentioned analysis, we can determine the soft masses of sleptons (m(L)) and CP-odd Higgs boson mass (mA). We have estimated that m(L) greater than or similar to 300 GeV and m(A) greater than or similar to 700 GeV. We have also commented on the allowed values of epsilon and b(epsilon), in case there is no enhancement in the Higgs to diphoton decay rate. Finally, we present a model to explain the smallness of epsilon and b(epsilon).

Photoluminescence decay rate engineering of CdSe quantum dots in ensemble arrays embedded with gold nano-antennae

We discuss experimental results on the ability to significantly tune the photoluminescence decay rates of CdSe quantum dots embedded in an ordered template, using lightly doped small gold nanoparticles (nano-antennae), of relatively low optical efficiency. We observe both enhancement and quenching of photoluminescence intensity of the quantum dots varying monotonically with increasing volume fraction of added gold nanoparticles, with respect to undoped quantum dot arrays. However, the corresponding variation in lifetime of photoluminescence spectra decay shows a hitherto unobserved, non-monotonic variation with gold nanoparticle doping. We also demonstrate that Purcell effect is quite effective for the larger (5 nm) gold nano-antenna leading to more than four times enhanced radiative rate at spectral resonance, for largest doping and about 1.75 times enhancement for off-resonance. Significantly for spectral off-resonance samples, we could simultaneously engineer reduction of non-radiative decay rate along with increase of radiative decay rate. Non-radiative decay dominates the system for the smaller (2 nm) gold nano-antenna setting the limit on how small these plasmonic nano-antennae could be to be effective in engineering significant enhancement in radiative decay rate and, hence, the overall quantum efficiency of quantum dot based hybrid photonic assemblies.

Exponential compact higher order scheme and their stability analysis for unsteady convection-diffusion equations

Exponential compact higher-order schemes have been developed for unsteady convection-diffusion equation (CDE). One of the developed scheme is sixth-order accurate which is conditionally stable for the Peclet number 0 <= Pe <= 2.8 and the other is fourth-order accurate which is unconditionally stable. Schemes for two-dimensional (2D) problems are made to use alternate direction implicit (ADI) algorithm. Example problems are solved and the numerical solutions are compared with the analytical solutions for each case.

Contribution of forest fire ash and plant litter decay on stream dissolved composition in a sub-humid tropical watershed (Mule Hole, Southern India)

The current understanding of wildfire effects on water chemistry is limited by the quantification of the elemental dissolution rates from ash and element release rate from the plant litter, as well as quantification of the specific ash contribution to stream water chemistry. The main objective of the study was to provide such knowledge through combination of experimental modelling, field data and end-member mixing analysis (EMMA) of wildfire impact on a watershed scale. The study concerns watershed effects of fire in the Indian subcontinent, a region that is typically not well represented in the fire science literature. In plant litter ash, major elements are either hosted in readily-soluble phases (K, Mg) such as salts, carbonates and oxides or in less-soluble carrier-phases (Si, Ca) such as amorphous silica, quartz and calcite. Accordingly, elemental release rates, inferred from ash leaching experiments in batch reactor, indicated that the element release into solution followed the order K > Mg > Na > Si > Ca. Experiments on plant litter leaching in mixed-flow reactor indicated two dissolution regimes: rapid, over the week and slower over the month. The mean dissolution rates at steady-state (R-ss) indicated that the release of major elements from plant litter followed the order Ca > Si > Cl > Mg > K > Na. R-ss for Si and Ca for tree leaves and herbaceous species are similar to those reported for boreal and European tree species and are higher than that from the dissolution of soil clay minerals. This identifies tropical plant litters as important source of Si and Ca for tropical surface waters. In the wildfire-impacted year 2004, the EMMA indicated that the streamflow composition (Ca, K, Mg, Na, Si, Cl) was controlled by four main sources: rainwater, throughfall, ash leaching and soil solution. The influence of the ash end-member was maximal early in the rainy season (the two first storm events) and decreased later in the rainy season, when the stream was dominated by the throughfall end-member. The contribution of plant litter decay to the streamwater composition for a year not impacted by wildfire is significant with estimated solute fluxes originating from this decay greatly exceed, for most major elements, the annual elemental dissolved fluxes at the Mule Hole watershed outlet. This highlighted the importance of solute retention and vegetation back uptake processes within the soil profile. Overall, the fire increased the mobility and export of major elements from the soils to the stream. It also shifted the vegetation-related contribution to the elemental fluxes at the watershed outlet from long-term (seasonal) to short-term (daily to monthly). (C) 2014 Elsevier B.V. All rights reserved.

Anomalous power law decay in solvation dynamics of DNA: a mode coupling theory analysis of ion contribution

Several time dependent fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long-range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect, yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Finally, we propose several experiments to check the prediction of the present theoretical work.

Diffusing-wave spectroscopy in an inhomogeneous object: Local viscoelastic spectra from ultrasound-assisted measurement of correlation decay arising from the ultrasound focal volume

We demonstrate diffusing-wave spectroscopy (DWS) in a localized region of a viscoelastically inhomogeneous object by measurement of the intensity autocorrelation g(2)(tau)] that captures only the decay introduced by the temperature-induced Brownian motion in the region. The region is roughly specified by the focal volume of an ultrasound transducer which introduces region specific mechanical vibration owing to insonification. Essential characteristics of the localized non-Markovian dynamics are contained in the decay of the modulation depth M(tau)], introduced by the ultrasound forcing in the focal volume selected, on g(2)(tau). The modulation depth M(tau(i)) at any delay time tau(i) can be measured by short-time Fourier transform of g(2)(tau) and measurement of the magnitude of the spectrum at the ultrasound drive frequency. By following the established theoretical framework of DWS, we are able to connect the decay in M(tau) to the mean-squared displacement (MSD) of scattering centers and the MSD to G*(omega), the complex viscoelastic spectrum. A two-region composite polyvinyl alcohol phantom with different viscoelastic properties is selected for demonstrating local DWS-based recovery of G*(omega) corresponding to these regions from the measured region specific M(tau(i))vs tau(i). The ultrasound-assisted measurement of MSD is verified by simulating, using a generalized Langevin equation (GLE), the dynamics of the particles in the region selected as well as by the usual DWS experiment without the ultrasound. It is shown that whereas the MSD obtained by solving the GLE without the ultrasound forcing agreed with its experimental counterpart covering small and large values of tau, the match was good only in the initial transients in regard to experimental measurements with ultrasound.

Photoluminescence, photocatalysis and Judd-Ofelt analysis of Eu3+-activated layered BiOCl phosphors

Eu3+-activated layered BiOCl phosphors were synthesized by the conventional solid-state method at relatively low temperature and shorter duration (400 degrees C for 1 h). All the samples were crystallized in the tetragonal structure with the space group P4/nmm (no. 129). Field emission scanning electron microscopy (FE-SEM) studies confirmed the plate-like morphology. Photoluminescence spectra exhibit characteristic luminescent D-5(0) -> F-7(J) (J = 0-4) intra-4f shell Eu3+ ion transitions. The electric dipole transition located at 620 nm (D-5(0) -> F-7(2)) was stronger than the magnetic dipole transition located at 594 nm (D-5(0) -> F-7(1)). The evaluated Commission International de l'Eclairage (CIE) color coordinates of Eu3+-activated BiOCl phosphors were close to the commercial Y2O3:Eu3+ and Y2O2S:Eu3+ red phosphors. Intensity parameters (Omega(2), Omega(4)) and various radiative properties such as transition probability (A(tot)), radiative lifetime (tau(rad)), stimulated emission cross-section (sigma(e)), gain bandwidth (sigma(e) x Delta lambda(eff)) and optical gain (sigma(e) x tau(rad)) were calculated using the Judd-Ofelt theory. The experimental decay curves of the D-5(0) level in Eu3+-activated BiOCl have a single exponential profile. In comparison with other Eu3+ doped materials, Eu3+-activated BiOCl phosphors have a long lifetime (tau(exp)), low non-radiative relaxation rate (W-NR), high quantum efficiency (eta) and better optical gain (sigma(e) x tau(rad)). The determined radiative properties revealed the usefulness of Eu3+-activated BiOCl in developing red lasers as well as optical display devices. Further, these samples showed efficient photocatalytic activity for the degradation of rhodamine B (RhB) dye under visible light irradiation. These photocatalysts are useful for the removal of toxic and non-biodegradable organic pollutants in water.

Mode coupling theory analysis of electrolyte solutions: Time dependent diffusion, intermediate scattering function, and ion solvation dynamics

A self-consistent mode coupling theory (MCT) with microscopic inputs of equilibrium pair correlation functions is developed to analyze electrolyte dynamics. We apply the theory to calculate concentration dependence of (i) time dependent ion diffusion, (ii) intermediate scattering function of the constituent ions, and (iii) ion solvation dynamics in electrolyte solution. Brownian dynamics with implicit water molecules and molecular dynamics method with explicit water are used to check the theoretical predictions. The time dependence of ionic self-diffusion coefficient and the corresponding intermediate scattering function evaluated from our MCT approach show quantitative agreement with early experimental and present Brownian dynamic simulation results. With increasing concentration, the dispersion of electrolyte friction is found to occur at increasingly higher frequency, due to the faster relaxation of the ion atmosphere. The wave number dependence of intermediate scattering function, F(k, t), exhibits markedly different relaxation dynamics at different length scales. At small wave numbers, we find the emergence of a step-like relaxation, indicating the presence of both fast and slow time scales in the system. Such behavior allows an intriguing analogy with temperature dependent relaxation dynamics of supercooled liquids. We find that solvation dynamics of a tagged ion exhibits a power law decay at long times-the decay can also be fitted to a stretched exponential form. The emergence of the power law in solvation dynamics has been tested by carrying out long Brownian dynamics simulations with varying ionic concentrations. The solvation time correlation and ion-ion intermediate scattering function indeed exhibit highly interesting, non-trivial dynamical behavior at intermediate to longer times that require further experimental and theoretical studies. (c) 2015 AIP Publishing LLC.

Correlation Between Decay Rate and Amplitude of Solar Cycles as Revealed from Observations and Dynamo Theory

Using different proxies of solar activity, we have studied the following features of the solar cycle: i) The linear correlation between the amplitude of cycle and its decay rate, ii) the linear correlation between the amplitude of cycle and the decay rate of cycle , and iii) the anti-correlation between the amplitude of cycle and the period of cycle . Features ii) and iii) are very useful because they provide precursors for future cycles. We have reproduced these features using a flux-transport dynamo model with stochastic fluctuations in the Babcock-Leighton effect and in the meridional circulation. Only when we introduce fluctuations in meridional circulation, are we able to reproduce different observed features of the solar cycle. We discuss the possible reasons for these correlations.

Looking for BSM physics using top-quark polarization and decay-lepton kinematic asymmetries

We explore beyond-standard-model (BSM) physics signatures in the l + jets channel of the t (t) over bar pair production process at the Tevatron and the LHC. We study the effects of BSM physics scenarios on the top-quark polarization and on the kinematics of the decay leptons. To this end, we construct asymmetries using the lepton energy and angular distributions. Further, we find their correlations with the top polarization, net charge asymmetry and top forward-backward asymmetry. We show that when used together, these observables can help discriminate effectively between SM and different BSM scenarios, which can lead to varying degrees of top polarization at the Tevatron as well as the LHC. We use two types of colored mediator models to demonstrate the effectiveness of proposed observables, an s-channel axigluon and a u-channel diquark.

An Exponential Lower Bound for Homogeneous Depth Four Arithmetic Formulas

We show here a 2(Omega(root d.log N)) size lower bound for homogeneous depth four arithmetic formulas. That is, we give an explicit family of polynomials of degree d on N variables (with N = d(3) in our case) with 0, 1-coefficients such that for any representation of a polynomial f in this family of the form f = Sigma(i) Pi(j) Q(ij), where the Q(ij)'s are homogeneous polynomials (recall that a polynomial is said to be homogeneous if all its monomials have the same degree), it must hold that Sigma(i,j) (Number of monomials of Q(ij)) >= 2(Omega(root d.log N)). The above mentioned family, which we refer to as the Nisan-Wigderson design-based family of polynomials, is in the complexity class VNP. Our work builds on the recent lower bound results 1], 2], 3], 4], 5] and yields an improved quantitative bound as compared to the quasi-polynomial lower bound of 6] and the N-Omega(log log (N)) lower bound in the independent work of 7].

Structure and Dynamics of Octamethylcyclotetrasiloxane Confined between Mica Surfaces

Using a molecular model for octamethylcydotetrasiloxane (OMCTS), molecular dynamics simulations are carried out to probe the phase state of OMCTS confined between two mica surfaces in equilibrium With a reservoir. Molecular dynamics simulations are carried out for elevations ranging from 5 to 35 K above the melting point for the OMCTS model used in this study. The Helmholtz free energy is, computed for a specific confinement using the :two-phase thermodynamic (2PT) method. Analysis of the in-plane pair correlation functions did not reveal signatures of freezing even under an extreme confinement of two layers. OMCTS is found to orient with a wide distribution of orientations with respect to the mica surface, with a distinct preference for the surface parallel configuration in the contact layers. The self-intermediate scattering function is found to decay with increasing relaxation times as the surface separation is decreased, and the two-step relaxation in the scattering function, a signature of glassy dynamics, distinctly evolves as the temperature is lowered. However, even at 5 K above the melting point, we did not observe a freezing transition and the self-intermediate scattering functions relax within 200 ps for the seven-layered confined system. The self diffusivity and relaxation times obtained from the Kohlrausch-Williams-Watts stretched exponential fits to the late alpha-relaxation exhibit power law scalings with the packing fraction as predicted by mode coupling theory. A distinct discontinuity in the Helmholtz free energy, potential energy, and a sharp change in the local bond order parameter, Q(4), was observed at 230 K for a five-layered system upon cooling, indicative of a first-order transition. A freezing point depression of about 30 K was observed for this five-layered confined system, and at the lower temperatures, contact layers were found to be disordered with long-range order present only in the inner layers. These dynamical signatures indicate that confined OMCTS undergoes a slowdown akin to a fluid approaching a glass transition upon increasing confinement, and freezing under confinement would require substantial subcooling below the bulk melting point of OMCTS.