Percolation on networks with antagonistic and dependent interactions

Drawing inspiration from real world interacting systems, we study a system consisting of two networks that exhibit antagonistic and dependent interactions. By antagonistic and dependent interactions we mean that a proportion of functional nodes in a network cause failure of nodes in the other, while failure of nodes in the other results in failure of links in the first. In contrast to interdependent networks, which can exhibit first-order phase transitions, we find that the phase transitions in such networks are continuous. Our analysis shows that, compared to an isolated network, the system is more robust against random attacks. Surprisingly, we observe a region in the parameter space where the giant connected components of both networks start oscillating. Furthermore, we find that for Erdos-Renyi and scale-free networks the system oscillates only when the dependence and antagonism between the two networks are very high. We believe that this study can further our understanding of real world interacting systems.

Role of boron diffusion in CoFeB/MgO magnetic tunnel junctions

Several scientific issues concerning the latest generation read heads for magnetic storage devices, based on CoFeB/MgO/CoFeBmagnetic tunnel junctions (MTJs) are known to be controversial, including such fundamental questions as to the behavior and the role of B in optimizing the physical properties of these devices. Quantitatively establishing the internal structures of several such devices with different annealing conditions using hard x-ray photoelectron spectroscopy, we resolve these controversies and establish that the B diffusion is controlled by the capping Ta layer, though Ta is physically separated from the layer with B by several nanometers. While explaining this unusual phenomenon, we also provide insight into why the tunneling magnetoresistance (TMR) is optimized at an intermediate annealing temperature, relating it to B diffusion, coupled with our studies based on x-ray diffraction and magnetic studies.

Metastable monoclinic and orthorhombic phases and electric field induced irreversible phase transformation at room temperature in the lead-free classical ferroelectric BaTiO3

For decades it has been a well-known fact that among the few ferroelectric compounds in the perovskite family, namely, BaTiO3, KNbO3, PbTiO3, and Na1/2Bi1/2TiO3, the dielectric and piezoelectric properties of BaTiO3 are considerably higher than the others in polycrystalline form at room temperature. Further, similar to ferroelectric alloys exhibiting morphotropic phase boundary, single crystals of BaTiO3 exhibit anomalously large piezoelectric response when poled away from the direction of spontaneous polarization at room temperature. These anomalous features in BaTiO3 remained unexplained so far from the structural standpoint. In this work, we have used high-resolution synchrotron x-ray powder diffraction, atomic resolution aberration-corrected transmission electron microscopy, in conjunction with a powder poling technique, to reveal that at 300 K (i) the equilibrium state of BaTiO3 is characterized by coexistence of metastable monoclinic Pm and orthorhombic (Amm2) phases along with the tetragonal phase, and (ii) strong electric field switches the polarization direction from the 001] direction towards the 101] direction. These results suggest that BaTiO3 at room temperature is within an instability regime, and that this instability is the fundamental factor responsible for the anomalous dielectric and piezoelectric properties of BaTiO3 as compared to the other homologous ferroelectric perovskite compounds at room temperature. Pure BaTiO3 at room temperature is therefore more akin to lead-based ferroelectric alloys close to the morphotropic phase boundary where polarization rotation and field induced ferroelectric-ferroelectric phase transformations play a fundamental role in influencing the dielectric and piezoelectric behavior.

Anisotropic isometric fluctuation relations in experiment and theory on a self-propelled rod

The isometric fluctuation relation (IFR) P. I. Hurtado et al., Proc. Natl. Acad. Sci. USA 108, 7704 (2011)] relates the relative probability of current fluctuations of fixed magnitude in different spatial directions. We test its validity in an experiment on a tapered rod, rendered motile by vertical vibration and immersed in a sea of spherical beads. We analyze the statistics of the velocity vector of the rod and show that they depart significantly from the IFR of Hurtado et al. Aided by a Langevin-equation model we show that our measurements are largely described by an anisotropic generalization of the IFR R. Villavicencio et al., Europhys. Lett. 105, 30009 (2014)], with no fitting parameters, but with a discrepancy in the prefactor whose origin may lie in the detailed statistics of the microscopic noise. The experimentally determined large-deviation function of the velocity vector has a kink on a curve in the plane.

Tunable mechanical and thermal properties of ZnS/CdS core/shell nanowires

Using all-atom molecular dynamics (MD) simulations, we have studied the mechanical properties of ZnS/CdS core/shell nanowires. Our results show that the coating of a few-atomic-layer CdS shell on the ZnS nanowire leads to a significant change in the stiffness of the core/shell nanowires compared to the stiffness of pure ZnS nanowires. The binding energy between the core and shell region decreases due to the lattice mismatch at the core-shell interface. This reduction in binding energy plays an important role in determining the stiffness of a core/shell nanowire. We have also investigated the effects of the shell on the thermal conductivity and melting behavior of the nanowires.

Descending from infinity: Convergence of tailed distributions

We investigate the relaxation of long-tailed distributions under stochastic dynamics that do not support such tails. Linear relaxation is found to be a borderline case in which long tails are exponentially suppressed in time but not eliminated. Relaxation stronger than linear suppresses long tails immediately, but may lead to strong transient peaks in the probability distribution. We also find that a delta-function initial distribution under stronger than linear decay displays not one but two different regimes of diffusive spreading.

Stability of the flow in a soft tube deformed due to an applied pressure gradient

A linear stability analysis is carried out for the flow through a tube with a soft wall in order to resolve the discrepancy of a factor of 10 for the transition Reynolds number between theoretical predictions in a cylindrical tube and the experiments of Verma and Kumaran J. Fluid Mech. 705, 322 (2012)]. Here the effect of tube deformation (due to the applied pressure difference) on the mean velocity profile and pressure gradient is incorporated in the stability analysis. The tube geometry and dimensions are reconstructed from experimental images, where it is found that there is an expansion and then a contraction of the tube in the streamwise direction. The mean velocity profiles at different downstream locations and the pressure gradient, determined using computational fluid dynamics, are found to be substantially modified by the tube deformation. The velocity profiles are then used in a linear stability analysis, where the growth rates of perturbations are calculated for the flow through a tube with the wall modeled as a neo-Hookean elastic solid. The linear stability analysis is carried out for the mean velocity profiles at different downstream locations using the parallel flow approximation. The analysis indicates that the flow first becomes unstable in the downstream converging section of the tube where the flow profile is more pluglike when compared to the parabolic flow in a cylindrical tube. The flow is stable in the upstream diverging section where the deformation is maximum. The prediction for the transition Reynolds number is in good agreement with experiments, indicating that the downstream tube convergence and the consequent modification in the mean velocity profile and pressure gradient could reduce the transition Reynolds number by an order of magnitude.

Study of motion around a static black hole in Einstein and Lovelock gravity

We study motion around a static Einstein and pure Lovelock black hole in higher dimensions. It is known that in higher dimensions bound orbits exist only for a pure Lovelock black hole in all even dimensions, D = 2N + 2, where N is the degree of Lovelock polynomial action. In particular, we compute periastron shift and light bending, and the latter is given by one of the transverse spatial components of the Riemann curvature tensor. We also consider the pseudo-Newtonian potentials and Kruskal coordinates.

Interphase anisotropy effects on lamellar eutectics: A numerical study

In directional solidification of binary eutectics, it is often observed that two-phase lamellar growth patterns grow tilted with respect to the direction z of the imposed temperature gradient. This crystallographic effect depends on the orientation of the two crystal phases alpha and beta with respect to z. Recently, an approximate theory was formulated that predicts the lamellar tilt angle as a function of the anisotropy of the free energy of the solid(alpha)-solid(beta) interphase boundary. We use two different numerical methods-phase field (PF) and dynamic boundary integral (BI)-to simulate the growth of steady periodic patterns in two dimensions as a function of the angle theta(R) between z and a reference crystallographic axis for a fixed relative orientation of alpha and beta crystals, that is, for a given anisotropy function (Wulff plot) of the interphase boundary. For Wulff plots without unstable interphase-boundary orientations, the two simulation methods are in excellent agreement with each other and confirm the general validity of the previously proposed theory. In addition, a crystallographic ``locking'' of the lamellae onto a facet plane is well reproduced in the simulations. When unstable orientations are present in the Wulff plot, it is expected that two distinct values of the tilt angle can appear for the same crystal orientation over a finite theta(R) range. This bistable behavior, which has been observed experimentally, is well reproduced by BI simulations but not by the PF model. Possible reasons for this discrepancy are discussed.

Role of supersymmetric heavy Higgs boson production in the self-coupling measurement of 125 GeV Higgs boson at the LHC

Measurement of the self-coupling of the 125 GeV Higgs boson is one of the most crucial tasks for a high luminosity run of the LHC, and it can only be measured in the di-Higgs final state. In the minimal supersymmetric standard model, heavy CP even Higgs (H) can decay into a lighter 125 GeV Higgs boson (h) and, therefore, can influence the rate of di-Higgs production. We investigate the role of single H production in the context of measuring the self-coupling of h. We have found that the H -> hh decay can change the value of Higgs (h) self-coupling substantially, in a low tan beta regime where the mass of the heavy Higgs boson lies between 250 and 600 GeV and, depending on the parameter space, it may be seen as an enhancement of the self-coupling of the 125 GeV Higgs boson.

Third post-Newtonian gravitational waveforms for compact binary systems in general orbits: Instantaneous terms

We compute the instantaneous contributions to the spherical harmonic modes of gravitational waveforms from compact binary systems in general orbits up to the third post-Newtonian (PN) order. We further extend these results for compact binaries in quasielliptical orbits using the 3PN quasi-Keplerian representation of the conserved dynamics of compact binaries in eccentric orbits. Using the multipolar post-Minkowskian formalism, starting from the different mass and current-type multipole moments, we compute the spin-weighted spherical harmonic decomposition of the instantaneous part of the gravitational waveform. These are terms which are functions of the retarded time and do not depend on the history of the binary evolution. Together with the hereditary part, which depends on the binary's dynamical history, these waveforms form the basis for construction of accurate templates for the detection of gravitational wave signals from binaries moving in quasielliptical orbits.

Two-dimensional homogeneous isotropic fluid turbulence with polymer additives

We carry out an extensive and high-resolution direct numerical simulation of homogeneous, isotropic turbulence in two-dimensional fluid films with air-drag-induced friction and with polymer additives. Our study reveals that the polymers (a) reduce the total fluid energy, enstrophy, and palinstrophy; (b) modify the fluid energy spectrum in both inverse-and forward-cascade regimes; (c) reduce small-scale intermittency; (d) suppress regions of high vorticity and strain rate; and (e) stretch in strain-dominated regions. We compare our results with earlier experimental studies and propose new experiments.

NMR relaxation studies in doped poly-3-methylthiophene

NMR relaxation rates (1/T-1), magnetic susceptibility, and electrical conductivity studies in doped poly-3-methylthiophene are reported in this paper. The magnetic susceptibility data show the contributions from both Pauli and Curie spins, with the size of the Pauli term depending strongly on the doping level. Proton and fluorine NMR relaxation rates have been studied as a function of temperature (3-300 K) and field (for protons at 0.9, 9.0, 16.4, and 23.4 T, and for fluorine at 9.0 T). The temperature dependence of T-1 is classified into three regimes: (a) For T < (g mu(B) B/2k(B)), the relaxation mechanism follows a modified Korringa relation due to electron-electron interactions and disorder. H-1-T-1 is due to the electron-nuclear dipolar interaction in addition to the contact term. (b) For the intermediate temperature range (g mu(B) B/2k(B)) < T < T-BPP (the temperature where the contribution from the reorientation motion to the T-1 is insignificant) the relaxation mechanism is via spin diffusion to the paramagnetic centers. (c) In the high-temperature regime and at low Larmor frequency the relaxation follows the modified Bloembergen, Purcell, and Pound model. T-1 data analysis has been carried out in light of these models depending upon the temperature and frequency range of study. Fluorine relaxation data have been analyzed and attributed to the PF6 reorientation. The cross relaxation among the H-1 and F-19 nuclei has been observed in the entire temperature range suggesting the role of magnetic dipolar interaction modulated by the reorientation of the symmetric molecular subgroups. The data analysis shows that the enhancement in the Korringa ratio is greater in a less conducting sample. Intra-and interchain hopping of charge carriers is found to be a dominant relaxation mechanism at low temperature. Frequency dependence of T-1(-1) on temperature shows that at low temperature T < (g mu(B) B/2k(B))] the system shows three dimensions and changes to quasi one dimension at high temperature. Moreover, a good correlation between electrical conductivity, magnetic susceptibility, and NMR T-1 data has been observed.

Spatial modulation of the composition of a binary liquid near a repulsive wall

When a binary liquid is confined by a strongly repulsive wall, the local density is depleted near the wall and an interface similar to that between the liquid and its vapor is formed. This analogy suggests that the composition of the binary liquid near this interface should exhibit spatial modulation similar to that near a liquid-vapor interface even if the interactions of the wall with the two components of the liquid are the same. The Guggenheim adsorption relation quantifies the concentrations of two components of a binary mixture near a liquid-vapor interface and qualitatively states that the majority (minority) component enriches the interface for negative (positive) mixing energy if the surface tensions of the two components are not very different. From molecular dynamics simulations of binary mixtures with different compositions and interactions we find that the Guggenheim relation is qualitatively satisfied at wall-induced interfaces for systems with negative mixing energy at all state points considered. For systems with positive mixing energy, this relation is found to be qualitatively valid at low densities, while it is violated at state points with high density where correlations in the liquid are strong. This observation is validated by a calculation of the density profiles of the two components of the mixture using density functional theory with the Ramakrishnan-Yussouff free-energy functional. Possible reasons for the violation of the Guggenheim relation are discussed.

Maximum group velocity in a one-dimensional model with a sinusoidally varying staggered potential

We use Floquet theory to study the maximum value of the stroboscopic group velocity in a one-dimensional tight-binding model subjected to an on-site staggered potential varying sinusoidally in time. The results obtained by numerically diagonalizing the Floquet operator are analyzed using a variety of analytical schemes. In the low-frequency limit we use adiabatic theory, while in the high-frequency limit the Magnus expansion of the Floquet Hamiltonian turns out to be appropriate. When the magnitude of the staggered potential is much greater or much less than the hopping, we use degenerate Floquet perturbation theory; we find that dynamical localization occurs in the former case when the maximum group velocity vanishes. Finally, starting from an ``engineered'' initial state where the particles (taken to be hard-core bosons) are localized in one part of the chain, we demonstrate that the existence of a maximum stroboscopic group velocity manifests in a light-cone-like spreading of the particles in real space.