How should biologists engage with controversial mathematical theory?

Mathematics is beautiful and precise and often necessary to understand complex biological phenomena. And yet biologists cannot always hope to fully understand the mathematical foundations of the theory they are using or testing. How then should biologists behave when mathematicians themselves are in dispute? Using the on-going controversy over Hamilton's rule as an example, I argue that biologists should be free to treat mathematical theory with a healthy dose of agnosticism. In doing so biologists should equip themselves with a disclaimer that publicly admits that they cannot entirely attest to the veracity of the mathematics underlying the theory they are using or testing. The disclaimer will only help if it is accompanied by three responsibilities - stay bipartisan in a dispute among mathematicians, stay vigilant and help expose dissent among mathematicians, and make the biology larger than the mathematics. I must emphasize that my goal here is not to take sides in the on-going dispute over the mathematical validity of Hamilton's rule, indeed my goal is to argue that we should refrain from taking sides.

Conformal perturbation theory and higher spin entanglement entropy on the torus

We study the free fermion theory in 1+1 dimensions deformed by chemical potentials for holomorphic, conserved currents at finite temperature and on a spatial circle. For a spin-three chemical potential mu, the deformation is related at high temperatures to a higher spin black hole in hs0] theory on AdS(3) spacetime. We calculate the order mu(2) corrections to the single interval Renyi and entanglement entropies on the torus using the bosonized formulation. A consistent result, satisfying all checks, emerges upon carefully accounting for both perturbative and winding mode contributions in the bosonized language. The order mu(2) corrections involve integrals that are finite but potentially sensitive to contact term singularities. We propose and apply a prescription for defining such integrals which matches the Hamiltonian picture and passes several non-trivial checks for both thermal corrections and the Renyi entropies at this order. The thermal corrections are given by a weight six quasi-modular form, whilst the Renyi entropies are controlled by quasi-elliptic functions of the interval length with modular weight six. We also point out the well known connection between the perturbative expansion of the partition function in powers of the spin-three chemical potential and the Gross-Taylor genus expansion of large-N Yang-Mills theory on the torus. We note the absence of winding mode contributions in this connection, which suggests qualitatively different entanglement entropies for the two systems.

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.

Novel Biobjective Clustering (BiGC) Based on Cooperative Game Theory

We propose a new approach to clustering. Our idea is to map cluster formation to coalition formation in cooperative games, and to use the Shapley value of the patterns to identify clusters and cluster representatives. We show that the underlying game is convex and this leads to an efficient biobjective clustering algorithm that we call BiGC. The algorithm yields high-quality clustering with respect to average point-to-center distance (potential) as well as average intracluster point-to-point distance (scatter). We demonstrate the superiority of BiGC over state-of-the-art clustering algorithms (including the center based and the multiobjective techniques) through a detailed experimentation using standard cluster validity criteria on several benchmark data sets. We also show that BiGC satisfies key clustering properties such as order independence, scale invariance, and richness.

Contrasting Electronic Requirements for C-H Binding and C-H Activation in d(6) Half-Sandwich Complexes of Rhenium and Tungsten

A computational study of the interaction half-sandwich metal fragments (metal=Re/W, electron count=d(6)), containing linear nitrosyl (NO+), carbon monoxide (CO), trifluorophosphine (PF3), N-heterocyclic carbene (NHC) ligands with alkanes are conducted using density functional theory employing the hybrid meta-GGA functional (M06). Electron deficiency on the metal increases with the ligand in the order NHC < CO < PF3 < NO+. Electron-withdrawing ligands like NO+ lead to more stable alkane complexes than NHC, a strong electron donor. Energy decomposition analysis shows that stabilization is due to orbital interaction involving charge transfer from the alkane to the metal. Reactivity and dynamics of the alkane fragment are facilitated by electron donors on the metal. These results match most of the experimental results known for CO and PF3 complexes. The study suggests activation of alkane in metal complexes to be facile with strong donor ligands like NHC. (C) 2015 Wiley Periodicals, Inc.

Anisotropic magnetic couplings and structure-driven canted to collinear transitions in Sr2IrO4 by magnetically constrained noncollinear DFT

We study the canted magnetic state in Sr2IrO4 using fully relativistic density functional theory (DFT) including an on-site Hubbard U correction. A complete magnetic phase diagram with respect to the tetragonal distortion and the rotation of IrO6 octahedra is constructed, revealing the presence of two types of canted to collinear magnetic transitions: a spin-flop transition with increasing tetragonal distortion and a complete quenching of the basal weak ferromagnetic moment below a critical octahedral rotation. Moreover, we put forward a scheme to study the anisotropic magnetic couplings by mapping magnetically constrained noncollinear DFT onto a general spin Hamiltonian. This procedure allows for the simultaneous account and direct control of the lattice, spin, and orbital interactions within a fully ab initio scheme. We compute the isotropic, single site anisotropy and Dzyaloshinskii-Moriya (DM) coupling parameters, and clarify that the origin of the canted magnetic state in Sr2IrO4 arises from the structural distortions and the competition between isotropic exchange and DM interactions.

Solid-state optical absorption from optimally tuned time-dependent range-separated hybrid density functional theory

We present a framework for obtaining reliable solid-state charge and optical excitations and spectra from optimally tuned range-separated hybrid density functional theory. The approach, which is fully couched within the formal framework of generalized Kohn-Sham theory, allows for the accurate prediction of exciton binding energies. We demonstrate our approach through first principles calculations of one- and two-particle excitations in pentacene, a molecular semiconducting crystal, where our work is in excellent agreement with experiments and prior computations. We further show that with one adjustable parameter, set to produce the known band gap, this method accurately predicts band structures and optical spectra of silicon and lithium fluoride, prototypical covalent and ionic solids. Our findings indicate that for a broad range of extended bulk systems, this method may provide a computationally inexpensive alternative to many-body perturbation theory, opening the door to studies of materials of increasing size and complexity.

Ultra-sensitive pressure dependence of bandgap of rutile-GeO2 revealed by many body perturbation theory

The reported values of bandgap of rutile GeO2 calculated by the standard density functional theory within local-density approximation (LDA)/generalized gradient approximation (GGA) show a wide variation (similar to 2 eV), whose origin remains unresolved. Here, we investigate the reasons for this variation by studying the electronic structure of rutile-GeO2 using many-body perturbation theory within the GW framework. The bandgap as well as valence bandwidth at Gamma-point of rutile phase shows a strong dependence on volume change, which is independent of bandgap underestimation problem of LDA/GGA. This strong dependence originates from a change in hybridization among O-p and Ge-(s and p) orbitals. Furthermore, the parabolic nature of first conduction band along X-Gamma-M direction changes towards a linear dispersion with volume expansion. (C) 2015 AIP Publishing LLC.

Evaluation of Precipitation Retrievals From Orbital Data Products of TRMM Over a Subtropical Basin in India

The spatial error structure of daily precipitation derived from the latest version 7 (v7) tropical rainfall measuring mission (TRMM) level 2 data products are studied through comparison with the Asian precipitation highly resolved observational data integration toward evaluation of the water resources (APHRODITE) data over a subtropical region of the Indian subcontinent for the seasonal rainfall over 6 years from June 2002 to September 2007. The data products examined include v7 data from the TRMM radiometer Microwave Imager (TMI) and radar precipitation radar (PR), namely, 2A12, 2A25, and 2B31 (combined data from PR and TMI). The spatial distribution of uncertainty from these data products were quantified based on performance metrics derived from the contingency table. For the seasonal daily precipitation over a subtropical basin in India, the data product of 2A12 showed greater skill in detecting and quantifying the volume of rainfall when compared with the 2A25 and 2B31 data products. Error characterization using various error models revealed that random errors from multiplicative error models were homoscedastic and that they better represented rainfall estimates from 2A12 algorithm. Error decomposition techniques performed to disentangle systematic and random errors verify that the multiplicative error model representing rainfall from 2A12 algorithm successfully estimated a greater percentage of systematic error than 2A25 or 2B31 algorithms. Results verify that although the radiometer derived 2A12 rainfall data is known to suffer from many sources of uncertainties, spatial analysis over the case study region of India testifies that the 2A12 rainfall estimates are in a very good agreement with the reference estimates for the data period considered.

Fluctuation theory of Rashba Fermi gases: Gaussian and beyond

Fermi gases with generalized Rashba spin-orbit coupling induced by a synthetic gauge field have the potential of realizing many interesting states, such as rashbon condensates and topological phases. Here, we address the key open problem of the fluctuation theory of such systems and demonstrate that beyond-Gaussian effects are essential to capture the finite temperature physics of such systems. We obtain their phase diagram by constructing an approximate non-Gaussian theory. We conclusively establish that spin-orbit coupling can enhance the exponentially small transition temperature (T-c) of a weakly attracting superfluid to the order of the Fermi temperature, paving a pathway towards high T-c superfluids.

Stability analysis of thick piezoelectric metal based FGM plate using first order and higher order shear deformation theory

This paper presents the stability analysis of functionally graded plate integrated with piezoelectric actuator and sensor at the top and bottom face, subjected to electrical and mechanical loading. The finite element formulation is based on first order and higher order shear deformation theory, degenerated shell element, von-Karman hypothesis and piezoelectric effect. The equation for static analysis is derived by using the minimum energy principle and solutions for critical buckling load is obtained by solving eigenvalue problem. The material properties of the functionally graded plate are assumed to be graded along the thickness direction according to simple power law function. Two types of boundary conditions are used, such as SSSS (simply supported) and CSCS (simply supported along two opposite side perpendicular to the direction of compression and clamped along the other two sides). Sensor voltage is calculated using present analysis for various power law indices and FG (functionally graded) material gradations. The stability analysis of piezoelectric FG plate is carried out to present the effects of power law index, material variations, applied mechanical pressure and piezo effect on buckling and stability characteristics of FG plate.

Pricing and Power allocation in Femtocells using Stackelberg Game Theory

We consider a system with multiple Femtocells operating in a Macrocell. The transmissions in one Femtocell interfere with its neighboring Femtocells as well as with the Macrocell Base Station. We model Femtocells as selfish nodes and the Macrocell Base Station protects itself by pricing subchannels for each usage. We use Stackelberg game model to study this scenario and obtain equilibrium policies that satisfy certain quality of service.

epsilon-expansions near three dimensions from conformal field theory

We formally extend the CFT techniques introduced in arXiv: 1505.00963, to phi(2d0/d0-2) theory in d = d(0) dimensions and use it to compute anomalous dimensions near d(0) = 3, 4 in a unified manner. We also do a similar analysis of the O(N) model in three dimensions by developing a recursive combinatorial approach for OPE contractions. Our results match precisely with low loop perturbative computations. Finally, using 3-point correlators in the CFT, we comment on why the phi(3) theory in d(0) = 6 is qualitatively different.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.