title: What is the mass of alpha Cen B b?

abstract: We investigate the possibility of constraining the sin i degeneracy of alpha Cen B b -- with orbital period P=3.24 d; a = 0.042 AU; msini = 1.1 M_earth -- to estimate the true mass of the newly reported terrestrial exoplanet in the nearest stellar system to our Sun. We present detailed numerical simulations of the dynamical stability of the exoplanet in the alpha Cen AB binary system for a range of initial inclinations, eccentricities, and semi-major axes. The system represents a benchmark case for the interplay of the Kozai mechanism, general relativistic and tidal forces. From our simulations, there is only a small boundary in initial inclinations and initial semi-major axes that result in the migration via the Kozai mechanism of alpha Cen B b to its present location. Inside this boundary, the planet orbit is stable for up to 1 Gyr against the Kozai mechanism, and outside this boundary the planet collides with alpha Cen B or is ejected. In our three simulations where the planet migrates in towards the star via the Kozai mechanism, the final inclination is 46-53 degrees relative to the AB orbital plane, lower than the initial inclination of 75 degrees in each case. We discuss inclination constraints from the formation of alpha Cen B b in situ at its present location, migration in a proto-planetary disk, or migration in resonance with additional planets. We conclude that alpha Cen B b probably has a mass of less than 2.7 M_earth, implying a likely terrestrial composition warranting future confirmation.


title: Multi-Step Cascade Annihilations of Dark Matter and the Galactic Center Excess

abstract: If dark matter is embedded in a non-trivial dark sector, it may annihilate and decay to lighter dark-sector states which subsequently decay to the Standard Model. Such scenarios - with annihilation followed by cascading dark-sector decays - can explain the apparent excess GeV gamma-rays identified in the central Milky Way, while evading bounds from dark matter direct detection experiments. Each 'step' in the cascade will modify the observable signatures of dark matter annihilation and decay, shifting the resulting photons and other final state particles to lower energies and broadening their spectra. We explore, in a model-independent way, the effect of multi-step dark-sector cascades on the preferred regions of parameter space to explain the GeV excess. We find that the broadening effects of multi-step cascades can admit final states dominated by particles that would usually produce too sharply peaked photon spectra; in general, if the cascades are hierarchical (each particle decays to substantially lighter particles), the preferred mass range for the dark matter is in all cases 20-150 GeV. Decay chains that have nearly-degenerate steps, where the products are close to half the mass of the progenitor, can admit much higher DM masses. We map out the region of mass/cross-section parameter space where cascades (degenerate, hierarchical or a combination) can fit the signal, for a range of final states. In the current work, we study multi-step cascades in the context of explaining the GeV excess, but many aspects of our results are general and can be extended to other applications.


title: A quantum peek inside the black hole event horizon

abstract: We solve the Klein-Gordon equation for a scalar field, in the background geometry of a dust cloud collapsing to form a black hole, everywhere in the (1+1) spacetime: that is, both inside and outside the event horizon and arbitrarily close to the curvature singularity. This allows us to determine the regularized stress tensor expectation value, everywhere in the appropriate quantum state (viz., the Unruh vacuum) of the field. We use this to study the behaviour of energy density and the flux measured in local inertial frames for the radially freely falling observer at any given event. Outside the black hole, energy density and flux lead to the standard results expected from the Hawking radiation emanating from the black hole, as the collapse proceeds. Inside the collapsing dust ball, the energy densities of both matter and scalar field diverge near the singularity in both (1+1) and (1+3) spacetime dimensions; but the energy density of the field dominates over that of classical matter. In the (1+3) dimensions, the total energy (of both scalar field and classical matter) inside a small spatial volume around the singularity is finite (and goes to zero as the size of the region goes to zero) but the total energy of the quantum field still dominates over that of the classical matter. Inside the event horizon, but \textit{outside} the collapsing matter, freely falling observers find that the energy density and the flux diverge close to the singularity. In this region, even the integrated energy inside a small spatial volume enclosing the singularity diverges. This result holds in both (1+1) and (1+3) spacetime dimensions with a \emph{milder} divergence for the total energy inside a small region in (1+3) dimensions. These results suggest that the back-reaction effects are significant even in the region \emph{outside the matter but inside the event horizon}, close to the singularity.


title: The EFT of Large Scale Structures at All Redshifts: Analytical Predictions for Lensing

abstract: We study the prediction of the Effective Field Theory of Large Scale Structures (EFTofLSS) for the matter power spectrum at different redshifts. In previous work, we found that the two-loop prediction can match the nonlinear power spectrum measured from $N$-body simulations at redshift zero within approximately 2% up to $k\sim 0.6\,h\, {\rm Mpc}^{-1}$ after fixing a single free parameter, the so-called "speed of sound". We determine the time evolution of this parameter by matching the EFTofLSS prediction to simulation output at different redshifts, and find that it is well-described by a fitting function that only includes one additional parameter. After the two free parameters are fixed, the prediction agrees with nonlinear data within approximately 2% up to at least $k\sim 1\,h\, {\rm Mpc}^{-1}$ at $z\geq 1$, and also within approximately 5% up to $k\sim 1.2\,h\, {\rm Mpc}^{-1}$ at $z=1$ and $k\sim 2.3\,h\, {\rm Mpc}^{-1}$ at $z=3$, a major improvement with respect to other perturbative techniques. We also develop an accurate way to estimate where the EFTofLSS predictions at different loop orders should fail, based on the sizes of the next-order terms that are neglected, and find agreement with the actual comparisons to data. Finally, we use our matter power spectrum results to perform analytical calculations of lensing potential power spectra corresponding to both CMB and galaxy lensing. This opens the door to future direct applications of the EFTofLSS to observations of gravitational clustering on cosmic scales.


title: M dwarf luminosity, radius, and $\alpha$-enrichment from $I$-band spectral features

abstract: Despite the ubiquity of M dwarfs and their growing importance to studies of exoplanets, Galactic evolution, and stellar structure, methods for precisely measuring their fundamental stellar properties remain elusive. Existing techniques for measuring M dwarf luminosity, mass, radius, or composition are calibrated over a limited range of stellar parameters or require expensive observations. We find a strong correlation between the $K_S$-band luminosity ($M_K$), the observed strength of the $I$-band sodium doublet absorption feature, and [Fe/H] in M dwarfs without strong H$\alpha$ emission. We show that the strength of this feature, coupled with [Fe/H] and spectral type, can be used to derive M dwarf $M_K$ and radius without requiring parallax. Additionally, we find promising evidence that the strengths of the $I$-band sodium doublet and the nearby $I$-band calcium triplet may jointly indicate $\alpha$-element enrichment. The use of these $I$-band features requires only moderate-resolution near-infrared spectroscopy to provide valuable information about the potential habitability of exoplanets around M dwarfs, and surface gravity and distance for M dwarfs throughout the Galaxy. This technique has immediate applicability for both target selection and candidate planet host system characterization for exoplanet missions such as \textit{TESS} and \textit{K2}.


title: Natural Inflation and Low Energy Supersymmetry

abstract: Natural (axionic) inflation provides a well-motivated and predictive scheme for the description of the early universe. It leads to sizeable primordial tensor modes and thus a high mass scale of the inflationary potential. Naively this seems to be at odds with low (TeV) scale supersymmetry, especially when embedded in superstring theory. We show that low scale supersymmetry is compatible with natural (high scale) inflation. The mechanism requires the presence of two axions that are provided through the moduli of string theory.


title: A PRIMAL view of the Milky Way, made possible by Gaia and M2M modelling

abstract: We have developed our original made-to-measure (M2M) algorithm, PRIMAL, with the aim of modelling the Galactic disc from upcoming Gaia data. From a Milky Way like N-body disc galaxy simulation, we have created mock Gaia data using M0III stars as tracers, taking into account extinction and the expected Gaia errors. In PRIMAL, observables calculated from the N-body model are compared with the target stars, at the position of the target stars. Using PRIMAL, the masses of the N-body model particles are changed to reproduce the target mock data, and the gravitational potential is automatically adjusted by the changing mass of the model particles. We have also adopted a new resampling scheme for the model particles to keep the mass resolution of the N-body model relatively constant. We have applied PRIMAL to this mock Gaia data and we show that PRIMAL can recover the structure and kinematics of a Milky Way like barred spiral disc, along with the apparent bar structure and pattern speed of the bar despite the galactic extinction and the observational errors.


title: Measuring Magnetic Fields Near and Far with the SKA via the Zeeman Effect

abstract: The measurement of Zeeman splitting in spectral lines---both in emission and absorption---can provide direct estimates of the magnetic field strength and direction in atomic and molecular clouds, both in our own Milky Way and in external galaxies. This method will probe the magnetic field in the warm and cold neutral components of the interstellar medium, providing a complement to the extensive SKA Faraday studies planning to probe the field in the ionized components.


title: Isospin violating dark matter in St\"uckelberg portal scenarios

abstract: Hidden sector scenarios in which dark matter (DM) interacts with the Standard Model matter fields through the exchange of massive Z' bosons are well motivated by certain string theory constructions. In this work, we thoroughly study the phenomenological aspects of such scenarios and find that they present a clear and testable consequence for direct DM searches. We show that such string motivated St\"uckelberg portals naturally lead to isospin violating interactions of DM particles with nuclei. We find that the relations between the DM coupling to neutrons and protons for both, spin-independent (fn/fp) and spin-dependent (an/ap) interactions, are very flexible depending on the charges of the quarks under the extra U(1) gauge groups. We show that within this construction these ratios are generically different from plus and minus 1 (i.e. different couplings to protons and neutrons) leading to a potentially measurable distinction from other popular portals. Finally, we incorporate bounds from searches for dijet and dilepton resonances at the LHC as well as LUX bounds on the elastic scattering of DM off nucleons to determine the experimentally allowed values of fn/fp and an/ap.


title: Solar Space Density of the Red Clump Stars and the Scale-length of the Thin Disc

abstract: We estimated the scale-length of the thin disc with the J and W1 magnitudes of the most probable Red Clump (RC) stars in the Galactic plane, $-0\overset{^\circ}.5 \leq b \leq +0\overset{^\circ}.5$, in 19 equal sized fields with consecutive Galactic longitudes which cover the interval $90^\circ \leq l \leq 270^\circ$. Our results are constrained with respect to the solar space density ($D^*=5.95$), which indicates that the radial variation of the density is lower for higher Galactocentric distances. The scale-length of the thin disc is 2 kpc for the fields in the Galactic anticentre direction or close to this direction, while it decreases continuously in the second and third quadrants reaching to a lower limit of $h$ = 1.6 kpc at the Galactic longitudes $l$ = 90$^\circ$ and $l$ = 270$^\circ$. The distribution of the scale-length in 19 fields is consistent with the predictions from the Galaxia model and its variation with longitude is probably due to the inhomogeneity structure of the disc caused by the accreted material or other features such as warp and flare.


title: Impact factor for high-energy two and three jets diffractive production

abstract: We present the calculation of the impact factor for the photon to quark, antiquark and gluon transition within Balitsky's shock-wave formalism. We also rederive the impact factor for photon to quark and antiquark transition. These results provide the necessary building blocks for further phenomenological studies of inclusive diffractive deep inelastic scattering as well as for two and three jets diffractive production which go beyond approximations discussed in the literature.


title: Evidence for feedback in action from the molecular gas content in the z~1.6 outflowing QSO XID2028

abstract: Gas outflows are believed to play a pivotal role in shaping galaxies, as they regulate both star formation and black hole growth. Despite their ubiquitous presence, the origin and the acceleration mechanism of such powerful and extended winds is not yet understood. Direct observations of the cold gas component in objects with detected outflows at other wavelengths are needed to assess the impact of the outflow on the host galaxy interstellar medium (ISM). We observed with the Plateau de Bure Interferometer an obscured quasar at z~1.5, XID2028, for which the presence of an ionised outflow has been unambiguously signalled by NIR spectroscopy. The detection of CO(3-2) emission in this source allows us to infer the molecular gas content and compare it to the ISM mass derived from the dust emission. We then analyze the results in the context of recent insights on scaling relations, which describe the gas content of the overall population of star-forming galaxies at a similar redshifts. The Star formation efficiency (~100) and gas mass (M_gas=2.1-9.5x10^{10} M_sun) inferred from the CO(3-2) line depend on the underlying assumptions on the excitation of the transition and the CO-to-H2 conversion factor. However, the combination of this information and the ISM mass estimated from the dust mass suggests that the ISM/gas content of XID2028 is significantly lower than expected for its observed M$_\star$, sSFR and redshift, based on the most up-to-date calibrations (with gas fraction <20% and depletion time scale <340 Myr). Overall, the constraints we obtain from the far infrared and millimeter data suggest that we are observing QSO feedback able to remove the gas from the host


title: Lower Bounds of Potential Blow-Up Solutions of the Three-dimensional Navier-Stokes Equations in $\dot{H}^\frac{3}{2}$

abstract: We improve previous known lower bounds for Sobolev norms of potential blow-up solutions to the three-dimensional Navier-Stokes equations in $\dot{H}^\frac{3}{2}$. We also present an alternate proof for the lower bound for the $\dot{H}^\frac{5}{2}$ blow-up.


title: A Hubble Astrometry Initiative: Laying the Foundation for the Next-Generation Proper-Motion Survey of the Local Group

abstract: High-precision astrometry throughout the Local Group is a unique capability of the Hubble Space Telescope (HST), with potential for transformative science, including constraining the nature of dark matter, probing the epoch of reionization, and understanding key physics of galaxy evolution. While Gaia will provide unparalleled astrometric precision for bright stars in the inner halo of the Milky Way, HST is the only current mission capable of measuring accurate proper motions for systems at greater distances (> 80 kpc), which represents the vast majority of galaxies in the Local Group. The next generation of proper-motion measurements will require long time baselines, spanning many years to decades and possibly multiple telescopes, combining HST with the James Webb Space Telescope (JWST) or the Wide-Field Infrared Survey Telescope (WFIRST). However, the current HST allocation process is not conducive to such multi-cycle/multi-mission science, which will bear fruit primarily over many years. We propose an HST astrometry initiative to enable long-time-baseline, multi-mission science, which we suggest could be used to provide comprehensive kinematic measurements of all dwarf galaxies and high surface-density stellar streams in the Local Group with HST's Advanced Camera for Surveys (ACS) or Wide Field Camera 3 (WFC3). Such an initiative not only would produce forefront scientific results within the next 5 years of HST's life, but also would serve as a critical anchor point for future missions to obtain unprecedented astrometric accuracy, ensuring that HST leaves a unique and lasting legacy for decades to come.


title: Protected quasi-locality in quantum systems with long-range interactions

abstract: We study the out-of-equilibrium dynamics of quantum systems with long-range interactions. Two different models describing, respectively, interacting lattice bosons and spins are considered. Our study relies on a combined approach based on accurate many-body numerical calculations as well as on a quasiparticle microscopic theory. For sufficiently fast decaying long-range potentials, we find that the quantum speed limit set by the long-range Lieb-Robinson bounds is never attained and a purely ballistic behavior is found. For slowly decaying potentials, a radically different scenario is observed. In the bosonic case, a remarkable local spreading of correlations is still observed, despite the existence of infinitely fast traveling excitations in the system. This is in marked contrast to the spin case, where locality is broken. We finally provide a microscopic justification of the different regimes observed and of the origin of the protected locality in the bosonic model.


title: The GAPS Programme with HARPS-N at TNG VIII: Observations of the Rossiter-McLaughlin effect and characterisation of the transiting planetary systems HAT-P-36 and WASP-11/HAT-P-10

abstract: We determine the true and the projected obliquity of HAT-P-36 and WASP-11/HAT-P-10 systems, respectively, which are both composed of a relatively cool star and a hot-Jupiter planet. Thanks to the high-resolution spectrograph HARPS-N, we observed the Rossiter-McLaughlin effect for both the systems by acquiring precise radial-velocity measurements during planetary transit events. We also present photometric observations comprising six light curves covering five transit events, obtained using three medium-class telescopes and the telescope-defocussing technique. One transit of WASP-11/HAT-P-10 was followed simultaneously from two observatories. The three transit light curves of HAT-P-36b show anomalies that are attributable to starspot complexes on the surface of the parent star, in agreement with the analysis of its spectra that indicate a moderate activity. By analysing the complete HATNet data set of HAT-P-36, we estimated the stellar rotation period by detecting a periodic photometric modulation in the light curve caused by star spots, obtaining Prot=15.3 days, which implies that the inclination of the stellar rotational axis with respect to the line of sight is 65 degree. We used the new spectroscopic and photometric data to revise the main physical parameters and measure the sky-projected misalignment angle of the two systems. We found \lambda=-14 degree for HAT-P-36 and \lambda=7 degree for WASP-11/HAT-P-10, indicating in both cases a good spin-orbit alignment. In the case of HAT-P-36, we also measured its real obliquity, which turned out to be 25 degrees.


title: Attending to lifelong learning skills through guided reflection in a physics class

abstract: This paper describes a tool, the Guided Reflection Form (GRF), which was used to promote reflection in a modeling-based physics course. Each week, students completed a guided reflection and received feedback from their instructors. These activities were intended to help students become better at the process of reflection, developing skills that they could apply in their future learning. We analyzed student reflections: (1) to provide insight into the reflection process itself and (2) to describe common themes in student reflections. Most students were able to use the GRF to reflect on their learning in meaningful ways. Moreover, the themes present in student reflections provide insights into struggles commonly faced by physics students. We discuss the design of the GRF in detail, so that others may use it as a tool to support student reflections.


title: Updated NNLO QCD predictions for the weak radiative B-meson decays

abstract: Weak radiative decays of the B mesons belong to the most important flavor changing processes that provide constraints on physics at the TeV scale. In the derivation of such constraints, accurate standard model predictions for the inclusive branching ratios play a crucial role. In the current Letter we present an update of these predictions, incorporating all our results for the O(alpha_s^2) and lower-order perturbative corrections that have been calculated after 2006. New estimates of nonperturbative effects are taken into account, too. For the CP- and isospin-averaged branching ratios, we find B_{s gamma} = (3.36 +_ 0.23) * 10^-4 and B_{d gamma} = 1.73^{+0.12}_{-0.22} * 10^-5, for E_gamma > 1.6GeV. Both results remain in agreement with the current experimental averages. Normalizing their sum to the inclusive semileptonic branching ratio, we obtain R_gamma = ( B_{s gamma} + B_{d gamma})/B_{c l nu} = (3.31 +_ 0.22) * 10^-3. A new bound from B_{s gamma} on the charged Higgs boson mass in the two-Higgs-doublet-model II reads M_{H^+} > 480 GeV at 95%C.L.


title: Probing Properties of Hot and Dense QCD Matter with Heavy Flavor in the PHENIX Experiment at RHIC

abstract: Hadrons carrying heavy quarks, i.e. charm or bottom, are important probes of the hot and dense medium created in relativistic heavy ion collisions. Heavy quark-antiquark pairs are mainly produced in initial hard scattering processes of partons. While some of the produced pairs form bound quarkonia, the vast majority hadronize into particles carrying open heavy flavor. Heavy quark production has been studied by the PHENIX experiment at RHIC via measurements of single leptons from semi-leptonic decays in both the electron channel at mid-rapidity and in the muon channel at forward rapidity. A large suppression and azimuthal anisotropy of single electrons have been observed in Au+Au collisions at 200 GeV. These results suggest a large energy loss and flow of heavy quarks in the hot, dense matter. The PHENIX experiment has also measured ${J/\psi}$ production at $\sqrt{s_{nn}}$ = 200 GeV in p+p, d+Au, Cu+Cu and Au+Au collisions, both at mid- and forward-rapidities, and additionally Cu+Au and U+U at forward-rapidities. In the most energetic collisions, more suppression is observed at forward rapidity than at central rapidity. This can be interpreted either as a sign of quark recombination, or as a hint of additional cold nuclear matter effects. The centrality dependence of nuclear modification factor, RAA, for ${J/\psi}$ in U+U collisions at $\sqrt{s_{nn}}$ = 193 GeV shows a similar trend to the lighter systems, Au+Au and Cu+Cu, at similar energy 200 GeV.


title: The (Q_7,Q_{1,2}) contribution to B -> X_s gamma at O(alpha_s^2)

abstract: Interference between the photonic dipole operator Q_7 and the current-current operators Q_{1,2} gives one of the most important QCD corrections to the B -> X_s gamma decay rate. So far, the O(alpha_s^2) part of this correction has been known in the heavy charm quark limit only (m_c >> m_b/2). Here, we evaluate this part at m_c=0, and use both limits in an updated phenomenological study. Our prediction for the CP- and isospin-averaged branching ratio in the Standard Model reads B^SM_{s gamma} = (3.36 +_ 0.23) * 10^-4 for E_gamma > 1.6 GeV.


title: The Atomic Hypothesis

abstract: A review is presented of the origin and development of the atomic hypothesis from antiquity till about the first millennium of the common era.


title: Correct-by-synthesis reinforcement learning with temporal logic constraints

abstract: We consider a problem on the synthesis of reactive controllers that optimize some a priori unknown performance criterion while interacting with an uncontrolled environment such that the system satisfies a given temporal logic specification. We decouple the problem into two subproblems. First, we extract a (maximally) permissive strategy for the system, which encodes multiple (possibly all) ways in which the system can react to the adversarial environment and satisfy the specifications. Then, we quantify the a priori unknown performance criterion as a (still unknown) reward function and compute an optimal strategy for the system within the operating envelope allowed by the permissive strategy by using the so-called maximin-Q learning algorithm. We establish both correctness (with respect to the temporal logic specifications) and optimality (with respect to the a priori unknown performance criterion) of this two-step technique for a fragment of temporal logic specifications. For specifications beyond this fragment, correctness can still be preserved, but the learned strategy may be sub-optimal. We present an algorithm to the overall problem, and demonstrate its use and computational requirements on a set of robot motion planning examples.


title: Inverse problem for Lagrangian systems on Lie algebroids and applications to reduction by symmetries

abstract: The language of Lagrangian submanifolds is used to extend a geometric characterization of the inverse problem of the calculus of variations on tangent bundles to regular Lie algebroids. Since not all closed sections are locally exact on Lie algebroids, the Helmholtz conditions on Lie algebroids are necessary but not sufficient, so they give a weaker definition of the inverse problem. As an application the Helmholtz conditions on Atiyah algebroids are obtained so that the relationship between the inverse problem and the reduced inverse problem by symmetries can be described. Some examples and comparison with previous approaches in the literature are provided.


title: The computer BESK and an early attempt to simulate galactic dynamics

abstract: The first N-body simulation of interacting galaxies, even producing spiral arms, was performed by Erik Holmberg in Lund (1941), not with a numerical computer, but by his arrangement of movable light-bulbs and photocells to measure the luminosity at each bulb and thereby estimate the gravitational force. A decade later, and with the arrival of the first programable computers, computations of galactic dynamics were performed, which were later transferred into a N-body simulation movie. I present here the background details for this work with a description of the important elements to note in the movie which may be retrieved at http://ttt.astro.su.se/~po .

categories: astro-ph.IM astro-ph.CO astro-ph.GA physics.hist-ph


title: A Meta-Algorithm for Creating Fast Algorithms for Counting ON Cells in Odd-Rule Cellular Automata

abstract: We develop a meta-algorithm that, given a polynomial (in one or more variables), and a prime p, produces a fast (logarithmic time) algorithm that takes a positive integer n and outputs the number of times each residue class modulo p appears as a coefficient when the polynomial is raised to the power n and the coefficients are read modulo p.


title: Hierarchy of Floquet gaps and edge states for driven honeycomb lattices

abstract: Electromagnetic driving in a honeycomb lattice can induce gaps and topological edge states with a structure of increasing complexity as the frequency of the driving lowers. While the high frequency case is the most simple to analyze we focus on the multiple photon processes allowed in the low frequency regime to unveil the hierarchy of Floquet edge-states. In the case of low intensities an analytical approach allows us to derive effective Hamiltonians and address the topological character of each gap in a constructive manner. At high intensities we obtain the net number of edge states, given by the winding number, with a numerical calculation of the Chern numbers of each Floquet band. Using these methods, we find a hierarchy that resembles that of a Russian nesting doll. This hierarchy classifies the gaps and the associated edge states in different orders according to the electron-photon coupling strength. For large driving intensities, we rely on the numerical calculation of the winding number, illustrated in a map of topological phase transitions. The hierarchy unveiled with the low energy effective Hamiltonians, alongside with the map of topological phase transitions discloses the complexity of the Floquet band structure in the low frequency regime. The proposed method for obtaining the effective Hamiltonian can be easily adapted to other Dirac Hamiltonians of two dimensional materials and even the surface of a 3D topological insulator.


title: Angular momentum evolution for galaxies

abstract: Using cosmological hydrodynamics simulations we study the angular momentum content of the simulated galaxies in relation with their morphological type. We found that not only the angular momentum of the disk component follow the expected theoretical relation, Mo, Mao & Whiye (1998), but also the spheroidal one, with a gap due to its lost of angular momentum, in agreement with Fall & Romanowsky (2013),. We also found that the galaxy size can plot in one general relation, despite the morphological type, as found by Kravtsov (2013).


title: Sums of four squares of primes

abstract: Let $E(N)$ denote the number of positive integers $n \le N$, with $n \equiv 4 \pmod{24}$, which cannot be represented as the sum of four squares of primes. We establish that $E(N)\ll N^{11/32}$, thus improving on an earlier result of Harman and the first author, where the exponent $7/20$ appears in place of $11/32$.


title: EmoNets: Multimodal deep learning approaches for emotion recognition in video

abstract: The task of the emotion recognition in the wild (EmotiW) Challenge is to assign one of seven emotions to short video clips extracted from Hollywood style movies. The videos depict acted-out emotions under realistic conditions with a large degree of variation in attributes such as pose and illumination, making it worthwhile to explore approaches which consider combinations of features from multiple modalities for label assignment. In this paper we present our approach to learning several specialist models using deep learning techniques, each focusing on one modality. Among these are a convolutional neural network, focusing on capturing visual information in detected faces, a deep belief net focusing on the representation of the audio stream, a K-Means based "bag-of-mouths" model, which extracts visual features around the mouth region and a relational autoencoder, which addresses spatio-temporal aspects of videos. We explore multiple methods for the combination of cues from these modalities into one common classifier. This achieves a considerably greater accuracy than predictions from our strongest single-modality classifier. Our method was the winning submission in the 2013 EmotiW challenge and achieved a test set accuracy of 47.67% on the 2014 dataset.


title: Weighted ${L^p}$-Liouville Theorems for Hypoelliptic Partial Differential Operators on Lie Groups

abstract: We prove weighted $L^p$-Liouville theorems for a class of second order hypoelliptic partial differential operators $\mathcal{L}$ on Lie groups $\mathbb{G}$ whose underlying manifold is $n$-dimensional space. We show that a natural weight is the right-invariant measure $\check{H}$ of $\mathbb{G}$. We also prove Liouville-type theorems for $C^2$ subsolutions in $L^p(\mathbb{G},\check{H})$. We provide examples of operators to which our results apply, jointly with an application to the uniqueness for the Cauchy problem for the evolution operator $\mathcal{L}-\partial_t$.


title: Game-theoretic approach to risk-sensitive benchmarked asset management

abstract: In this article we consider a game theoretic approach to the Risk-Sensitive Benchmarked Asset Management problem (RSBAM) of Davis and Lleo \cite{DL}. In particular, we consider a stochastic differential game between two players, namely, the investor who has a power utility while the second player represents the market which tries to minimize the expected payoff of the investor. The market does this by modulating a stochastic benchmark that the investor needs to outperform. We obtain an explicit expression for the optimal pair of strategies as for both the players.


title: On the Bumpy Metrics Theorem for Minimal Submanifolds

abstract: This paper proves several natural generalizations of the theorem that for a generic, $C^k$ Riemannian metric on a smooth manifold, there are no closed, embedded, minimal submanifolds with nontrivial jacobi fields.


title: Frequency Domain TOF: Encoding Object Depth in Modulation Frequency

abstract: Time of flight cameras may emerge as the 3-D sensor of choice. Today, time of flight sensors use phase-based sampling, where the phase delay between emitted and received, high-frequency signals encodes distance. In this paper, we present a new time of flight architecture that relies only on frequency---we refer to this technique as frequency-domain time of flight (FD-TOF). Inspired by optical coherence tomography (OCT), FD-TOF excels when frequency bandwidth is high. With the increasing frequency of TOF sensors, new challenges to time of flight sensing continue to emerge. At high frequencies, FD-TOF offers several potential benefits over phase-based time of flight methods.


title: Periodic and chaotic orbits of plane-confined micro-rotors in creeping flows

abstract: We explore theoretically the complex dynamics and emergent behaviors of spinning spheres immersed in viscous fluid. The particles are coupled to each-other via the fluid in which they are suspended: each particle disturbs the surrounding fluid with a rotlet field and that fluid flow affects the motion of the other particles. We notice the emergence of intricate periodic or chaotic trajectories that depend on the rotors initial position and separation. The point-rotor motions confined to a plane bear similarities the classic 2D point-vortex dynamics. Our analyses highlight the complexity of the interaction between just a few rotors and suggest richer behavior in denser populations. We discuss how the model gives insight into more complex systems and suggest possible extensions for future theoretical studies.

categories: physics.flu-dyn cond-mat.soft math.DS nlin.CD


title: Restricted linear congruences

abstract: In this paper, using properties of Ramanujan sums and of the discrete Fourier transform of arithmetic functions, we give an explicit formula for the number of solutions of the linear congruence $a_1x_1+\cdots +a_kx_k\equiv b \pmod{n}$, with $\gcd(x_i,n)=t_i$ ($1\leq i\leq k$), where $a_1,t_1,\ldots,a_k,t_k, b,n$ ($n\geq 1$) are arbitrary integers. As a consequence, we derive necessary and sufficient conditions under which the above restricted linear congruence has no solutions. The number of solutions of this kind of congruence was first considered by Rademacher in 1925 and Brauer in 1926, in the special case of $a_i=t_i=1$ $(1\leq i \leq k)$. Since then, this problem has been studied, in several other special cases, in many papers; in particular, Jacobson and Williams [{\it Duke Math. J.} {\bf 39} (1972), 521--527] gave a nice explicit formula for the number of such solutions when $(a_1,\ldots,a_k)=t_i=1$ $(1\leq i \leq k)$. The problem is very well-motivated and has found intriguing applications in several areas of mathematics, computer science, and physics, and there is promise for more applications/implications in these or other directions.


title: Non-spurious solutions to discrete boundary value problems through variational methods

abstract: Using direct variational method we consider the existence of non-spurious solutions to the following Dirichlet problem $\ddot{x}\left( t\right) =f\left( t,x\left( t\right) \right) $, $x\left( 0\right) =x\left( 1\right) =0 $ where $f:\left[ 0,1\right] \times \mathbb{R} \rightarrow \mathbb{R}$ is a jointly continuous function convex in $x$ which does not need to satisfy any further growth conditions.


title: Detecting temperature fluctuations at equilibrium

abstract: Gibbs and Boltzmann definitions of temperature agree only in the macroscopic limit. The ambiguity in identifying the equilibrium temperature of a finite sized `small' system exchanging energy with a bath is usually understood as a limitation of conventional statistical mechanics. We interpret this ambiguity as resulting from a stochastically fluctuating temperature coupled with the phase space variables giving rise to a broad temperature distribution. With this ansatz, we develop the equilibrium statistics and dynamics of small systems. Numerical evidence using an analytically tractable model shows that the effects of temperature fluctuations can be detected in equilibrium and dynamical properties of the phase space of the small system. Our theory generalizes statistical mechanics to small systems relevant to biophysics and nanotechnology.


title: A note on a global invertibility of mappings on $R^{n}$

abstract: We provide sufficient conditions for a mapping $f:R^{n}\rightarrow R^{n}$ to be a global diffeomorphism in case it is strictly (Hadamard) differentiable. We use classical local invertibility conditions together with the non-smooth critical point theory.


title: 1D topological chains with Majorana fermions in 2D non-topological optical lattices

abstract: The recent experimental realization of 1D equal Rashba-Dresselhaus spin-orbit coupling (ERD-SOC) for cold atoms provide a disorder-free platform for the implementation and observation of Majorana fermions (MFs), similar to the well studied solid state nanowire/superconductor heterostructures. However, the corresponding 1D chains of cold atoms possess strong quantum fluctuation, which may destroy the superfluids and MFs. In this Letter, we show that such 1D topological chains with MFs may be on demand generated in a 2D non-topological optical lattice with 1D ERD-SOC by modifying local potentials on target locations using experimentally already implemented atomic gas microscopes or patterned (e.g., double or triple well) optical lattices. All ingredients in our scheme have been experimentally realized and the combination of them may pave the way for the experimental observation of MFs in a clean system.


title: Optimally Combining Classifiers Using Unlabeled Data

abstract: We develop a worst-case analysis of aggregation of classifier ensembles for binary classification. The task of predicting to minimize error is formulated as a game played over a given set of unlabeled data (a transductive setting), where prior label information is encoded as constraints on the game. The minimax solution of this game identifies cases where a weighted combination of the classifiers can perform significantly better than any single classifier.


title: Ontology-Based Quality Evaluation of Value Generalization Hierarchies for Data Anonymization

abstract: In privacy-preserving data publishing, approaches using Value Generalization Hierarchies (VGHs) form an important class of anonymization algorithms. VGHs play a key role in the utility of published datasets as they dictate how the anonymization of the data occurs. For categorical attributes, it is imperative to preserve the semantics of the original data in order to achieve a higher utility. Despite this, semantics have not being formally considered in the specification of VGHs. Moreover, there are no methods that allow the users to assess the quality of their VGH. In this paper, we propose a measurement scheme, based on ontologies, to quantitatively evaluate the quality of VGHs, in terms of semantic consistency and taxonomic organization, with the aim of producing higher-quality anonymizations. We demonstrate, through a case study, how our evaluation scheme can be used to compare the quality of multiple VGHs and can help to identify faulty VGHs.


title: On some metabelian 2-group whose abelianization is of type (2, 2, 2) and applications

abstract: Let $G$ be some metabelian $2$-group satisfying the condition $G/G'\simeq \mathbb{Z}/2\mathbb{Z}\times\mathbb{Z}/2\mathbb{Z}\times\mathbb{Z}/2\mathbb{Z}$. In this paper, we construct all the subgroups of $G$ of index $2$ or $4$, we give the abelianization types of these subgroups and we compute the kernel of the transfer map. Then we apply these results to study the capitulation problem of the $2$-ideal classes of some fields $\mathbf{k}$ satisfying the condition $\mathrm{G}al(\mathbf{k}_2^{(2)}/\mathbf{k})\simeq G$, where $\mathbf{k}_2^{(2)}$ is the second Hilbert $2$-class field of $\mathbf{k}$.


title: Hierarchical Bond Graph Modelling of Biochemical Networks

abstract: The bond graph approach to modelling biochemical networks is extended to allow hierarchical construction of complex models from simpler components. This is made possible by representing the simpler components as thermodynamically open systems exchanging mass and energy via ports. A key feature of this approach is that the resultant models are robustly thermodynamically compliant: the thermodynamic compliance is not dependent on precise numerical values of parameters. Moreover, the models are reusable due to the well-defined interface provided by the energy ports. To extract bond graph model parameters from parameters found in the literature, general and compact formulae are developed to relate free-energy constants and equilibrium constants. The existence and uniqueness of solutions is considered in terms of fundamental properties of stoichiometric matrices. The approach is illustrated by building a hierarchical bond graph model of glycogenolysis in skeletal muscle.


title: Finding a Hamiltonian cycle by finding the global minimizer of a linearly constrained problem

abstract: It has been shown that a global minimizer of a smooth determinant of a matrix function corresponds to the largest cycle of a graph. When it exists, this is a Hamiltonian cycle. Finding global minimizers even of a smooth function is a challenge. The difficulty is often exacerbated by the existence of many global minimizers. One may think this would help, but in the case of Hamiltonian cycles the ratio of the number of global minimizers to the number of local minimizers is typically astronomically small. There are various equivalent forms of the problem and here we report on two. Although the focus is on finding Hamiltonian cycles, and this has an interest in and of itself, this is just a proxy for a class of problems that have discrete variables. The solution of relaxations of these problems is typically at a degenerate vertex, and in the neighborhood of the solution the Hessian is indefinite. The form of the Hamiltonian cycle problem we address has the virtue of being an ideal test problem for algorithms designed for discrete nonlinear problems in general. It is easy to generate problems of varying size and varying character, and they have the advantage of being able to determine if a global solution has been found. A feature of many discrete problems is that there are many solutions. For example, in the frequency assignment problem any permutation of a solution is also a solution. A consequence is that a common characteristic of the relaxed problems is that they have large numbers of global minimizers and even larger numbers of both local minimizers, and saddle points whose reduced Hessian has only a single negative eigenvalue. Efficient algorithms that seek to find global minimizers for this type of problem are described. Results using BONMIN, a solver for nonlinear problems with continuous and discrete variables, are also included.


title: On a surface formed by randomly gluing together polygonal discs

abstract: Starting with a collection of $n$ oriented polygonal discs, with an even number $N$ of sides in total, we generate a random oriented surface by randomly matching the sides of discs and properly gluing them together. Encoding the surface in a random permutation $\gamma$ of $[N]$, we use the Fourier transform on $S_N$ to show that $\gamma$ is asymptotic to the permutation distributed uniformly on the alternating group $A_N$ ($A_N^c$ resp.) if $N-n$ and $N/2$ are of the same (opposite resp.) parity. We use this to prove a local central limit theorem for the number of vertices on the surface, whence for its Euler characteristic $\chi$. We also show that with high probability the random surface consists of a single component, and thus has a well-defined genus $g=1-\chi/2$, which is asymptotic to a Gaussian random variable, with mean $(N/2-n-\log N)/2$ and variance $(\log N)/2$.


title: YFCC100M: The New Data in Multimedia Research

abstract: We present the Yahoo Flickr Creative Commons 100 Million Dataset (YFCC100M), the largest public multimedia collection that has ever been released. The dataset contains a total of 100 million media objects, of which approximately 99.2 million are photos and 0.8 million are videos, all of which carry a Creative Commons license. Each media object in the dataset is represented by several pieces of metadata, e.g. Flickr identifier, owner name, camera, title, tags, geo, media source. The collection provides a comprehensive snapshot of how photos and videos were taken, described, and shared over the years, from the inception of Flickr in 2004 until early 2014. In this article we explain the rationale behind its creation, as well as the implications the dataset has for science, research, engineering, and development. We further present several new challenges in multimedia research that can now be expanded upon with our dataset.


title: Stability of Discrete time Recurrent Neural Networks and Nonlinear optimization problems

abstract: We consider the method of Reduction of Dissipativity Domain to prove global Lyapunov stability of Discrete Time Recurrent Neural Networks. The standard and advanced criteria for Absolute Stability of these essentially nonlinear systems produce rather weak results. The method mentioned above is proved to be more powerful. It involves a multi-step procedure with maximization of special nonconvex functions over polytopes on every step. We derive conditions which guarantee an existence of at most one point of local maximum for such functions over every hyperplane. This nontrivial result is valid for wide range of neuron transfer functions.


title: Recovering Nonlocal Differential Pencils

abstract: Inverse problems for differential pencils with nonlocal conditions are investigated. Several uniqueness theorems of inverse problems from the Weyl-type function and spectra are proved, which are generalizations of the well-known Weyl function and Borg's inverse problem for the classical Sturm-Liouville operator.


title: Latent Hierarchical Model for Activity Recognition

abstract: We present a novel hierarchical model for human activity recognition. In contrast to approaches that successively recognize actions and activities, our approach jointly models actions and activities in a unified framework, and their labels are simultaneously predicted. The model is embedded with a latent layer that is able to capture a richer class of contextual information in both state-state and observation-state pairs. Although loops are present in the model, the model has an overall linear-chain structure, where the exact inference is tractable. Therefore, the model is very efficient in both inference and learning. The parameters of the graphical model are learned with a Structured Support Vector Machine (Structured-SVM). A data-driven approach is used to initialize the latent variables; therefore, no manual labeling for the latent states is required. The experimental results from using two benchmark datasets show that our model outperforms the state-of-the-art approach, and our model is computationally more efficient.


title: Attractors for Damped Semilinear Wave Equations with a Robin--Acoustic Boundary Perturbation

abstract: Under consideration is the damped semilinear wave equation \[ u_{tt}+u_t-\Delta u + u + f(u)=0 \] on a bounded domain $\Omega$ in $\mathbb{R}^3$ with a perturbation parameter $\varepsilon>0$ occurring in an acoustic boundary condition, limiting ($\varepsilon=0$) to a Robin boundary condition. With minimal assumptions on the nonlinear term $f$, the existence and uniqueness of global weak solutions is shown for each $\varepsilon\in[0,1]$. Also, the existence of a family of global attractors is shown to exist. After proving a general result concerning the upper-semicontinuity of a one-parameter family of sets, the result is applied to the family of global attractors. Because of the complicated boundary conditions for the perturbed problem, fractional powers of the Laplacian are not well-defined; moreover, because of the restrictive growth assumptions on $f$, the family of global attractors is obtained from the asymptotic compactness method developed by J. Ball for generalized semiflows. With more relaxed assumptions on the nonlinear term $f$, we are able to show the global attractors possess optimal regularity and prove the existence of an exponential attractor, for each $\varepsilon\in[0,1].$ This result insures that the corresponding global attractor inherits finite (fractal) dimension, however, the dimension is {\em{not}} necessarily uniform in $\varepsilon$.


title: A Noncommutative Borsuk-Ulam Theorem for Natsume-Olsen Spheres

abstract: Natsume-Olsen noncommutative spheres are C*-algebras which generalize C(S^k) when k is odd. These algebras admit natural actions by finite cyclic groups, and if one of these actions is fixed, any equivariant homomorphism between two Natsume-Olsen spheres of the same dimension induces a nontrivial map on odd K-theory. This result is an extended, noncommutative Borsuk-Ulam theorem in odd dimension, and just as in the topological case, this theorem has many (almost) equivalent formulations in terms of theta-deformed spheres of arbitrary dimension. In addition, we present theorems on graded Banach algebras, motivated by algebraic Borsuk-Ulam results of A. Taghavi.


title: Low-Frequency 1/f Noise in MoS2 Thin-Film Transistors: Comparison of Single and Multilayer Structures

abstract: We report on the transport and low-frequency noise measurements of MoS2 thin-film transistors with "thin" (2-3 atomic layers) and "thick" (15-18 atomic layers) channels. The back-gated transistors made with the relatively thick MoS2 channels have advantages of the higher electron mobility and lower noise level. The normalized noise spectral density of the low-frequency 1/f noise in "thick" MoS2 transistors is of the same level as that in graphene. The MoS2 transistors with the atomically thin channels have substantially higher noise levels. It was established that, unlike in graphene devices, the noise characteristics of MoS2 transistors with "thick" channels (15-18 atomic planes) could be described by the McWhorter model. Our results indicate that the channel thickness optimization is crucial for practical applications of MoS2 thin-film transistors.


title: Deep Clustered Convolutional Kernels

abstract: Deep neural networks have recently achieved state of the art performance thanks to new training algorithms for rapid parameter estimation and new regularization methods to reduce overfitting. However, in practice the network architecture has to be manually set by domain experts, generally by a costly trial and error procedure, which often accounts for a large portion of the final system performance. We view this as a limitation and propose a novel training algorithm that automatically optimizes network architecture, by progressively increasing model complexity and then eliminating model redundancy by selectively removing parameters at training time. For convolutional neural networks, our method relies on iterative split/merge clustering of convolutional kernels interleaved by stochastic gradient descent. We present a training algorithm and experimental results on three different vision tasks, showing improved performance compared to similarly sized hand-crafted architectures.


title: Review of Thermal Properties of Graphene and Few-Layer Graphene: Applications in Electronics

abstract: We review thermal properties of graphene and few-layer graphene, and discuss applications of these materials in thermal management of advanced electronics. The intrinsic thermal conductivity of graphene - among the highest of known materials - is dominated by phonons near the room temperature. The examples of thermal management applications include the few-layer graphene heat spreaders integrated near the heat generating areas of the high-power density transistors. It has been demonstrated that few-layer graphene heat spreaders can lower the hot-spot temperature during device operation resulting in improved performance and reliability of the devices.


title: The Semiclassical Einstein Equation on Cosmological Spacetimes

abstract: The subject of this thesis is the coupling of quantum fields to a classical gravitational background in a semiclassical fashion. It contains a thorough introduction into quantum field theory on curved spacetime with a focus on the stress-energy tensor and the semiclassical Einstein equation. Basic notions of differential geometry, topology, functional and microlocal analysis, causality and general relativity will be summarised, and the algebraic approach to QFT on curved spacetime will be reviewed. Apart from these foundations, the original research of the author and his collaborators will be presented: Together with Fewster, the author studied the up and down structure of permutations using their decomposition into so-called atomic permutations. The relevance of these results to this thesis is their application in the calculation of the moments of quadratic quantum fields. In a work with Pinamonti, the author showed the local and global existence of solutions to the semiclassical Einstein equation in flat cosmological spacetimes coupled to a scalar field by solving simultaneously for the quantum state and the Hubble function in an integral-functional equation. The theorem is proved with a fixed-point theorem using the continuous functional differentiability and boundedness of the integral kernel of the integral-functional equation. In another work with Pinamonti the author proposed an extension of the semiclassical Einstein equations which couples the moments of a stochastic Einstein tensor to the moments of the quantum stress-energy tensor. In a toy model of a Newtonianly perturbed exponentially expanding spacetime it is shown that the quantum fluctuations of the stress-energy tensor induce an almost scale-invariant power spectrum for the perturbation potential and that non-Gaussianties arise naturally.


title: ALMA Observations of the IRDC Clump G34.43+00.24 MM3: DNC/HNC Ratio

abstract: We have observed the clump G34.43+00.24 MM3 associated with an infrared dark cloud in DNC $J$=3--2, HN$^{13}$C $J$=3--2, and N$_2$H$^+$ $J$=3--2 with the Atacama Large Millimeter/submillimeter Array (ALMA). The N$_2$H$^+$ emission is found to be relatively weak near the hot core and the outflows, and its distribution is clearly anti-correlated with the CS emission. This result indicates that a young outflow is interacting with cold ambient gas. The HN$^{13}$C emission is compact and mostly emanates from the hot core, whereas the DNC emission is extended around the hot core. Thus, the DNC and HN$^{13}$C emission traces warm regions near the protostar differently. The DNC emission is stronger than the HN$^{13}$C emission toward most parts of this clump. The DNC/HNC abundance ratio averaged within a $15^{\prime\prime} \times 15^{\prime\prime}$ area around the phase center is higher than 0.06. This ratio is much higher than the value obtained by the previous single-dish observations of DNC and HN$^{13}$C $J$=1--0 ($\sim$0.003). It seems likely that the DNC and HNC emission observed with the single-dish telescope traces lower density envelopes, while that observed with ALMA traces higher density and highly deuterated regions. We have compared the observational results with chemical-model results in order to investigate the behavior of DNC and HNC in the dense cores. Taking these results into account, we suggest that the low DNC/HNC ratio in the high-mass sources obtained by the single-dish observations are at least partly due to the low filling factor of the high density regions.


title: Localizing spin dynamics in a spin-1 Bose-Einstein condensate via magnetic pulses

abstract: Spin exchange interaction between atoms in a spin-1 Bose-Einstein condensate causes atomic spin evolving periodically under the single spatial mode approximation in the mean field theory. By applying fast magnetic pulses according to a two-step or a four-step control protocol, we find analytically that the spin dynamics is significantly suppressed for an arbitrary initial state. Numerical calculations under single mode approximation are carried out to confirm the validity and robustness of these protocols. This localization method can be readily utilized to improve the sensitivity of a magnetometer based on spin-1 Bose-Einstein condensates.


title: Kepler Eclipsing Binary Stars. VI. Identification of Eclipsing Binaries in the K2 Campaign 0 Data-set

abstract: The original {\it Kepler} mission observed and characterized over 2400 eclipsing binaries in addition to its prolific exoplanet detections. Despite the mechanical malfunction and subsequent non-recovery of two reaction wheels used to stabilize the instrument, the {\it Kepler} satellite continues collecting data in its repurposed {\it K2} mission surveying a series of fields along the ecliptic plane. Here we present an analysis of the first full baseline {\it K2} data release: the Campaign 0 data-set. In the 7761 light curves, we have identified a total of 207 eclipsing binaries. Of these, 97 are new discoveries that were not previously identified. Our pixel-level analysis of these objects has also resulted in identification of several false positives (observed targets contaminated by neighboring eclipsing binaries), as well as the serendipitous discovery of two short period exoplanet candidates. We provide catalog cross-matched source identifications, orbital periods, morphologies and ephemerides for these eclipsing systems. We also describe the incorporation of the K2 sample into the Kepler Eclipsing Binary Catalog\footnote{\url{keplerebs.villanova.edu/k2}}, present spectroscopic follow-up observations for a limited selection of nine systems, and discuss prospects for upcoming {\it K2} campaigns.


title: Proton ordering dynamics of H2O ice

abstract: From high precision measurements of the complex dielectric constant of H2O ice, we identify the critical temperatures of the phase transition into and out of ice XI from ice Ih to occur at T_Ih-IX=58.9 K and T_IX-Ih=73.4 K. For D2O, T_Ih-IX=63.7 K and T_IX-Ih=78.2 K. A triple point is identified to exist at 0.07 GPa and 73.4 K for H2O and 0.08 GPa and 78.2 K for D2O where ices Ih, II and XI coexist. A first order phase transition with kinetic broadening associated to proton ordering dynamics is identified at 100 K.


title: Generalized Random Simplicial Complexes

abstract: We consider a multi-parameter model for randomly constructing simplicial complexes. This model interpolates between random clique complexes and Linial-Meshulam random $k$-dimensional complexes, two models that have been extensively studied. While these models asymptotically exhibit nontrivial cohomology in only one or two dimensions, we show that in this generalized setting nontrivial cohomology can occur in several dimensions simultaneously. We establish upper and lower thresholds for the appearance of nontrivial cohomology in a particular dimension, and in some instances characterize the behavior at criticality.


title: Linear Global Translation Estimation with Feature Tracks

abstract: This paper derives a novel linear position constraint for cameras seeing a common scene point, which leads to a direct linear method for global camera translation estimation. Unlike previous solutions, this method deals with collinear camera motion and weak image association at the same time. The final linear formulation does not involve the coordinates of scene points, which makes it efficient even for large scale data. We solve the linear equation based on $L_1$ norm, which makes our system more robust to outliers in essential matrices and feature correspondences. We experiment this method on both sequentially captured images and unordered Internet images. The experiments demonstrate its strength in robustness, accuracy, and efficiency.


title: Brauer Algebras of multiply laced Weyl type

abstract: In this paper, we will define the Brauer algebras of Weyl types, and describe some propositions of these algebras. Especially, we prove the result of type $G_2$ to accomplish our project of Brauer algebras of non-simply laced types.


title: Design of the Helium Purifier for IHEP-ADS Helium Purification System

abstract: Helium Purification System is an important sub-system in the Accelerator Driven Subcritical System of the Institute of High Energy Physics(IHEP ADS). The purifier is designed to work at the temperature of 77K. The purifier will work in a flow rate of 5g/s at 20MPa in continuous operation of 12 hours. The oil and moisture are removed by coalescing filters and a dryer, while nitrogen and oxygen are condensed by a phase separator and then adsorbed in several activated carbon adsorption cylinders. After purification, the purified helium has an impurity content of less than 5ppm.


title: Construction by bosonization of a fermion-phonon model

abstract: We discuss an extension of the (massless) Thirring model describing interacting fermions in one dimension which are coupled to phonons and where all interactions are local. This fermion-phonon model can be solved exactly by bosonization. We present a construction and solution of this model which is mathematically rigorous by treating it as a limit of a Luttinger-phonon model. A self-contained account of the mathematical results underlying bosonization is included, together with complete proofs.


title: The Nuclear Near-Infrared Spectral Properties of Nearby Galaxies

abstract: We present spectra of the nuclear regions of 50 nearby (D = 1 - 92 Mpc, median = 20 Mpc) galaxies of morphological types E to Sm. The spectra, obtained with the Gemini Near-IR Spectrograph on the Gemini North telescope, cover a wavelength range of approximately 0.85-2.5 microns at R~1300--1800. There is evidence that most of the galaxies host an active galactic nucleus (AGN), but the range of AGN luminosities (log (L2-10 keV [erg/s]) = 37.0-43.2) in the sample means that the spectra display a wide variety of features. Some nuclei, especially the Seyferts, exhibit a rich emission-line spectrum. Other objects, in particular the type 2 Low Ionisation Nuclear Emission Region galaxies, show just a few, weak emission lines, allowing a detailed view of the underlying stellar population. These spectra display numerous absorption features sensitive to the stellar initial mass function, as well as molecular bands arising in cool stars, and many other atomic absorption lines. We compare the spectra of subsets of galaxies known to be characterised by intermediate-age and old stellar populations, and find clear differences in their absorption lines and continuum shapes. We also examine the effect of atmospheric water vapor on the signal-to-noise ratio achieved in regions between the conventional NIR atmospheric windows, of potential interest to those planning observations of redshifted emission lines or other features affected by telluric H2O. Further exploitation of this data set is in progress, and the reduced spectra and data reduction tools are made available to the community.


title: Combinatorial rigidity and independence of generalized pinned subspace-incidence constraint systems

abstract: Given a hypergraph $H$ with $m$ hyperedges and a set $X$ of $m$ \emph{pins}, i.e.\ globally fixed subspaces in Euclidean space $\mathbb{R}^d$, a \emph{pinned subspace-incidence system} is the pair $(H, X)$, with the constraint that each pin in $X$ lies on the subspace spanned by the point realizations in $\mathbb{R}^d$ of vertices of the corresponding hyperedge of $H$. We are interested in combinatorial characterization of pinned subspace-incidence systems that are \emph{minimally rigid}, i.e.\ those systems that are guaranteed to generically yield a locally unique realization. As is customary, this is accompanied by a characterization of generic independence as well as rigidity. In a previous paper \cite{sitharam2014incidence}, we used pinned subspace-incidence systems towards solving the \emph{fitted dictionary learning} problem, i.e.\ dictionary learning with specified underlying hypergraph, and gave a combinatorial characterization of minimal rigidity for a more restricted version of pinned subspace-incidence system, with $H$ being a uniform hypergraph and pins in $X$ being 1-dimension subspaces. Moreover in a recent paper \cite{Baker2015}, the special case of pinned line incidence systems was used to model biomaterials such as cellulose and collagen fibrils in cell walls. In this paper, we extend the combinatorial characterization to general pinned subspace-incidence systems, with $H$ being a non-uniform hypergraph and pins in $X$ being subspaces with arbitrary dimension. As there are generally many data points per subspace in a dictionary learning problem, which can only be modeled with pins of dimension larger than $1$, such an extension enables application to a much larger class of fitted dictionary learning problems.


title: Encoding Source Language with Convolutional Neural Network for Machine Translation

abstract: The recently proposed neural network joint model (NNJM) (Devlin et al., 2014) augments the n-gram target language model with a heuristically chosen source context window, achieving state-of-the-art performance in SMT. In this paper, we give a more systematic treatment by summarizing the relevant source information through a convolutional architecture guided by the target information. With different guiding signals during decoding, our specifically designed convolution+gating architectures can pinpoint the parts of a source sentence that are relevant to predicting a target word, and fuse them with the context of entire source sentence to form a unified representation. This representation, together with target language words, are fed to a deep neural network (DNN) to form a stronger NNJM. Experiments on two NIST Chinese-English translation tasks show that the proposed model can achieve significant improvements over the previous NNJM by up to +1.08 BLEU points on average


title: The combinatorial algorithm for computing $\pi(x)$

abstract: This paper describes recent advances in the combinatorial method for computing $\pi(x)$, the number of primes $\leq x$. In particular, the memory usage has been reduced by a factor of $\log x$, and modifications for shared- and distributed-memory parallelism have been incorporated. The resulting method computes $\pi(x)$ with complexity $O(x^{2/3}\mathrm{log}^{-2}x)$ in time and $O(x^{1/3}\mathrm{log}^{2}x)$ in space. The algorithm has been implemented and used to compute $\pi(10^n)$ for $1 \leq n \leq 26$ and $\pi(2^m)$ for $1\leq m \leq 86$. The mathematics presented here is consistent with and builds on that of previous authors.


title: Spherical polytropic balls cannot mimic black holes

abstract: The so-called black hole shadow is a dark region which is expected to appear in a fine image of optical observation of black holes. It is essentially an absorption cross section of black hole, and the boundary of shadow is determined by unstable circular orbits of photons (UCOP). If there exists a compact object possessing UCOP but no black hole horizon, it can provide us with the same shadow image with black holes, and a detection of shadow image cannot be a direct evidence of black hole existence. Then, this paper examine whether or not such compact objects can exist under some suitable conditions. We investigate thoroughly the static spherical polytropic ball of perfect fluid with single polytrope index, and then investigate a representative example of the piecewise polytropic ball. Our result is that the spherical polytropic ball which we have investigated cannot possess UCOP, if the sound speed at center is subluminal (slower-than-light). This means that, if the polytrope treated in this paper is a good model of stellar matter in compact objects, the detection of shadow image is regarded as a good evidence of black hole existence. As a by-product, we have found the upper bound of the mass-to-radius radio (M/R) of polytropic ball with single index, M/R < 0.281, under the subluminal-sound-speed condition.


title: A Study of Multi-frequency Polarization Pulse Profiles of Millisecond Pulsars

abstract: We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $\lambda^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.


title: CEoptim: Cross-Entropy R Package for Optimization

abstract: The cross-entropy (CE) method is simple and versatile technique for optimization, based on Kullback-Leibler (or cross-entropy) minimization. The method can be applied to a wide range of optimization tasks, including continuous, discrete, mixed and constrained optimization problems. The new package CEoptim provides the R implementation of the CE method for optimization. We describe the general CE methodology for optimization and well as some useful modifications. The usage and efficacy of CEoptim is demonstrated through a variety of optimization examples, including model fitting, combinatorial optimization, and maximum likelihood estimation.


title: Host-feeding enhances stability of discrete-time host-parasitoid population dynamic models

abstract: Discrete-time models are the traditional approach for capturing population dynamics of a host-parasitoid system. Recent work has introduced a semi-discrete framework for obtaining model update functions that connect host-parasitoid population levels from year-to-year. In particular, this framework uses differential equations to describe the hosts-parasitoid interaction during the time of year where they come in contact, allowing specific behaviors to be mechanistically incorporated into the model. We use the semi-discrete approach to study the effects of host-feeding, which occurs when a parasitoid consumes a potential host larva without ovipositing. Our results show that host-feeding by itself cannot stabilize the system, and both the host and parasitoid populations exhibit diverging oscillations similar to the Nicholson-Bailey model. However, when combined with other stabilizing mechanisms such as density-dependent host mortality or density-dependent parasitoid attack rate, host-feeding expands the region of parameter space that allows for a stable host-parasitoid equilibrium. Finally, our results show that host-feeding causes inefficiency in the parasitoid population, which yields a higher population of hosts per generation. This suggests that host-feeding may have limited long-term impact in terms of suppressing host levels for biological control applications.


title: NMR evidence for inhomogeneous glassy behavior driven by nematic fluctuations in iron arsenide superconductors

abstract: We present $^{75}$As nuclear magnetic resonance spin-lattice and spin-spin relaxation rate data in Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$ and Ba(Fe$_{1-x}$Cu$_x$)$_2$As$_2$ as a function of temperature, doping and magnetic field. The relaxation curves exhibit a broad distribution of relaxation rates, consistent with inhomogeneous glassy behavior up to 100 K. The doping and temperature response of the width of the dynamical heterogeneity is similar to that of the nematic susceptibility measured by elastoresistance measurements. We argue that quenched random fields which couple to the nematic order give rise to a nematic glass that is reflected in the spin dynamics.


abstract: Hadronic decays of the highly excited $2D$ $D_s$ resonances have been studied in the $^3P_0$ model. Widths of all possible hadronic decay channels of the $2D$ $D_s$ have been computed. $D^*_{s1}(2700)$, $D^*_{s1}(2860)$, $D^*_{s3}(2860)$, $D(2600)$ and $D(2750)$ can be produced from hadronic decays of the $2D$ $D_s$, and relevant hadronic decay widths have been particularly paid attention to. The hadronic decay widths of $2D$ $D_s$ to $D(2600)$ or $D(2750)$ may be large, and the numerical results are different in different assignments of $D(2600)$ and $D(2750)$. The hadronic decay widths of $2D$ $D_s$ to $D^*_{s1}(2860)$, $D^*_{s3}(2860)$ or $D^*_{s1}(2700)$ are very small, and different in different assignments of $D^*_{s1}(2700)$.


title: HI emission and absorption in nearby, gas-rich galaxies

abstract: We present the results of a targeted search for intervening HI absorption in six nearby, gas-rich galaxies using the Australia Telescope Compact Array (ATCA). The sightlines searched have impact parameters of 10-20 kpc. By targeting nearby galaxies we are also able to map their HI emission, allowing us to directly relate the absorption-line detection rate to the extended HI distribution. The continuum sightlines intersect the HI disk in four of the six galaxies, but no intervening absorption was detected. Of these four galaxies, we find that three of the non-detections are the result of the background source being too faint. In the fourth case we find that the ratio of the spin temperature to the covering factor ($T_{\mathrm{S}}/f$) must be much higher than expected ($\gtrsim$5700 K) in order to explain the non-detection. We discuss how the structure of the background continuum sources may have affected the detection rate of HI absorption in our sample, and the possible implications for future surveys. Future work including an expanded sample, and VLBI observations, would allow us to better investigate the expected detection rate, and influence of background source structure, on the results of future surveys.


title: Estimation of the parameters of an infectious disease model using neural networks

abstract: In this paper, we propose a realistic mathematical model taking into account the mutual interference among the interacting populations. This model attempts to describe the control (vaccination) function as a function of the number of infective individuals, which is an improvement over the existing susceptible infective epidemic models. Regarding the growth of the epidemic as a nonlinear phenomenon we have developed a neural network architecture to estimate the vital parameters associated with this model. This architecture is based on a recently developed new class of neural networks known as co-operative and supportive neural networks. The application of this architecture to the present study involves preprocessing of the input data, and this renders an efficient estimation of the rate of spread of the epidemic. It is observed that the proposed new neural network outperforms a simple feed-forward neural network and polynomial regression.


title: SDP-based Joint Sensor and Controller Design for Information-regularized Optimal LQG Control

abstract: We consider a joint sensor and controller design problem for linear Gaussian stochastic systems in which a weighted sum of quadratic control cost and the amount of information acquired by the sensor is minimized. This problem formulation is motivated by situations where a control law must be designed in the presence of sensing, communication, and privacy constraints. We show that the optimal joint sensor-controller design is relatively easy when the sensing policy is restricted to be linear. Namely, an explicit form of the optimal linear sensor equation, the Kalman filter, and the certainty equivalence controller that jointly solves the problem can be efficiently found by semidefinite programming (SDP). Whether the linearity assumption in our design is restrictive or not is currently an open problem.


title: Anti-telephones in transformation optics: metamaterials with closed null geodesics

abstract: We apply the methods of transformation optics to theoretical descriptions of spacetimes that support closed null geodesic curves. The metric used is based on frame dragging spacetimes, such as the van Stockum dust or the Kerr black hole. Through transformation optics, this metric is analogous to a material that in theory should allow for communication between past and future. Presented herein is a derivation and description of the spacetime and the resulting permeability, permittivity, and magneto-electric couplings that a material would need in order for light in the material to follow closed null geodesics. We also address the paradoxical implications of such a material, and demonstrate why such a material would not actually result in a violation of causality. A full derivation of the Plebanski equations is also included.


title: What is User Experience Really: towards a UX Conceptual Framework

abstract: For more then a decade the term User Experience (UX) has been highly debated and defined in many ways. However, often UX remains as a vague concept and it may be hard to understand the very nature of it. In this paper we aimed at providing a better understanding of this concept. We explored the multi-faceted UX literature, reviewing the current state-of- the-art knowledge and emphasizing the multi-dimensional nature of the concept. Based on the literature review we built a conceptual framework of UX using the elements that are linked to it and reported in different studies. To show the potential use of the framework, we examined the UX delivered by different phone applications on different mobile devices using the elements in the framework. Several interesting insights have been obtained in terms of how the phone applications deliver different UX. Our study opens up a promising line of investigating the contemporary meaning of UX.


title: Cosmological constraints from Subaru weak lensing cluster counts

abstract: We present results of weak lensing cluster counts obtained from 11 sq.deg SuprimeCam data. Although the area is much smaller than previous work dealing with weak lensing peak statistics, the number density of galaxies usable for weak lensing analysis is about twice as large as those. The higher galaxy number density reduces the noise in the weak lensing mass maps, and thus increases the signal-to-noise ratio of peaks of the lensing signal due to massive clusters. This enables us to construct a weak lensing selected cluster sample by adopting a high threshold S/N, such that the contamination rate due to false signals is small. We find 6 peaks with S/N>5. For all the peaks, previously identified clusters of galaxies are matched within a separation of 1 arcmin, demonstrating good correspondence between the peaks and clusters of galaxies. We evaluate the statistical error using mock weak lensing data, and find Npeak=6+/-3.1 in an effective area of 9.0 sq.deg. We compare the measured weak lensing cluster counts with the theoretical model prediction based on halo models and place the constraint on Omega_m-sigma_8 plane which is found to be consistent with currently standard LCDM models. It is demonstrated that the weak lensing cluster counts can place a unique constraint on sigma_8-c_0 plane, where c_0 is the normalization of the dark matter halo mass-concentration relationship. Finally we discuss prospects for ongoing/future wide field optical galaxy surveys.


title: Log-stable law of energy dissipation as a framework of turbulence intermittency

abstract: To describe the small-scale intermittency of turbulence, a self-similarity is assumed for the probability density function of a logarithm of the rate of energy dissipation smoothed over a length scale among those in the inertial range. The result is an extension of Kolmogorov's classical theory in 1941, i.e., a one-parameter framework where the logarithm obeys some stable distribution. Scaling laws are obtained for the dissipation rate and for the two-point velocity difference. They are consistent with theoretical constraints and with the observed scaling laws. Also discussed is the physics that determines the value of the parameter.


title: Measurement of (alpha,n) reaction cross sections of erbium isotopes for testing astrophysical rate predictions

abstract: The $\gamma$-process in core-collapse and/or type Ia supernova explosions is thought to explain the origin of the majority of the so-called $p$ nuclei (the 35 proton-rich isotopes between Se and Hg). Reaction rates for $\gamma$-process reaction network studies have to be predicted using Hauser-Feshbach statistical model calculations. Recent investigations have shown problems in the prediction of $\alpha$-widths at astrophysical energies which are an essential input for the statistical model. It has an impact on the reliability of abundance predictions in the upper mass range of the $p$ nuclei. With the measurement of the $^{164,166}$Er($\alpha$,n)$^{167,169}$Yb reaction cross sections at energies close to the astrophysically relevant energy range we tested the recently suggested low energy modification of the $\alpha$+nucleus optical potential in a mass region where $\gamma$-process calculations exhibit an underproduction of the $p$ nuclei. Using the same optical potential for the $\alpha$-width which was derived from combined $^{162}$Er($\alpha$,n) and $^{162}$Er($\alpha$,$\gamma$) measurement makes it plausible that a low-energy modification of the optical $\alpha$+nucleus potential is needed.


title: Analysis of Discrete Choice Models: A Welfare-Based Framework

abstract: Based on the observation that many existing discrete choice models admit a welfare function of utilities whose gradient gives the choice probability vector, we propose a new representation of discrete choice model which we call the welfare-based choice model. The welfare-based choice model is meaningful on its own by providing a new way of constructing choice models. More importantly, it provides great analysis convenience for establishing connections among existing choice models. We prove by using convex analysis theory, that the welfare-based choice model is equivalent to the representative agent choice model and the semi-parametric choice model, establishing the equivalence of the latter two. We show that these three models are all strictly more general than the random utility model, while when there are only two alternatives, those four models are equivalent. In particular, we show that the distinction between the welfare-based choice model and the random utility model lies in the requirement of the higher-order derivatives of the welfare function. We then define a new concept in choice models: substitutability/complementarity between alternatives. We show that the random utility model only allows substitutability between different alternatives; while the welfare-based choice model allows flexible substitutability/complementarity patterns. We argue that such flexibility could be desirable in capturing certain practical choice patterns and expanding the scope of discrete choice models. Examples are given of new choice models proposed under our framework.


title: Vacuum Rabi spectra of a single quantum emitter

abstract: We report the observation of the vacuum Rabi splitting of a single quantum emitter by measuring its direct spontaneous emission into free space. We used a semiconductor quantum dot inside a photonic crystal nanocavity, in conjunction with an appropriate cavity design and filtering with a polarizer and an aperture, enabling the extraction of the inherently-weak emitter's signal. The emitter's vacuum Rabi spectra exhibit clear differences to those measured by detecting the cavity photon leakage. Moreover, we observed an asymmetric vacuum Rabi spectrum induced by interference between the emitter and cavity detection channels. Our observations lay the groundwork for accessing various cavity quantum electrodynamics phenomena that manifest themselves only in the emitter's direct spontaneous emission.


title: Discovery of a Disk Gap Candidate at 20 AU in TW Hydrae

abstract: We present a new Subaru/HiCIAO high-contrast H-band polarized intensity (PI) image of a nearby transitional disk associated with TW Hydrae. The scattered light from the disk was detected from 0.2" to 1.5" (11 - 81 AU) and the PI image shows a clear axisymmetric depression in polarized intensity at ~ 0.4" (~ 20 AU) from the central star, similar to the ~ 80 AU gap previously reported from HST images. Azimuthal polarized intensity profile also shows the disk beyond 0.2" is almost axisymmetric. We discuss two possible scenarios explaining the origin of the polarized intensity depression: 1) a gap structure may exist at ~ 20 AU from the central star because of shallow slope seen in the polarized intensity profile, and 2) grain growth may be occurring in the inner region of the disk. Multi-band observations at NIR and millimeter/sub-millimeter wavelengths play a complementary role in investigating dust opacity and may help reveal the origin of the gap more precisely.


title: Neutrino Masses and Flavor Mixing

abstract: We discuss the neutrino oscillations, using texture zero mass matrices for the leptons. The reactor mixing angle $\theta^{}_{l}$ is calculated. The ratio of the masses of two neutrinos is determined by the solar mixing angle. We can calculate the masses of the three neutrinos: $m_1$ $\approx$ 0.003 eV - $m_2$ $\approx$ 0.012 eV - $m_3$ $\approx$ 0.048 eV.


title: Classification of Charge Density Waves Based on Their Nature

abstract: The concept of a Charge Density Wave (CDW) permeates much of condensed matter physics and chemistry. Conceptually, CDWs have their origin rooted in the instability of a one-dimensional system described by Peierls. The extension of this concept to reduced dimensional systems has led to the concept of Fermi surface nesting (FSN), which dictates the wave vector q_cdw of the CDW and the corresponding lattice distortion. The idea is that segments of the Fermi contours are connected by q_cdw, resulting in the effective screening of phonons inducing Kohn Anomalies in their dispersion at q_cdw, driving a lattice restructuring at low temperatures. There is growing theoretical and experimental evidence that this picture fails in many real systems and in fact it is the momentum dependence of the electron-phonon coupling (EPC) matrix element that determines the characteristic of CDW phase (q_cdw). Here, based on the published results for the prototypical CDW system 2H-NbSe2, we show how well the q-dependent EPC matrix element, but not the FSN, can describe the origin of CDW. We further demonstrate a procedure of combing electronic band and phonon (dispersion and linewidth) measurements to extract the EPC matrix element, allowing the electronic states involved in the EPC to be identified. Thus we show that a large EPC does not necessarily induce the CDW phase, with Bi2Sr2CaCu2O8 (Bi2212) as the example, and the charge ordered phenomena observed in various cuprates are not driven by FSN or EPC. To experimentally resolve the microscopic picture of EPC will lead to a fundamental change in the way we think about, write about, and classify Charge Density Waves.


title: Constructing ab initio models of ultra-thin Al-AlOx-Al barriers

abstract: The microscopic structure of ultra-thin oxide barriers often plays a major role in modern nano-electronic devices. In the case of superconducting electronic circuits, their operation depends on the electrical non-linearity provided by one or more such oxide layers in the form of ultra-thin tunnel barriers (also known as Josephson junctions). Currently available fabrication techniques manufacture an amorphous oxide barrier, which is attributed as a major noise source within the device. The nature of this noise is currently an open question and requires both experimental and theoretical investigation. Here, we present a methodology for constructing atomic scale computational models of Josephson junctions using a combination of molecular mechanics, empirical and ab initio methods. These junctions consist of ultra-thin amorphous aluminium-oxide layers sandwiched between crystalline aluminium. The stability and structure of these barriers as a function of density and stoichiometry are investigated, which we compare to experimentally observed parameters


title: Using member galaxy luminosities as halo mass proxies of galaxy groups

abstract: Reliable halo mass estimation for a given galaxy system plays an important role both in cosmology and galaxy formation studies. Here we set out to find the way that can improve the halo mass estimation for those galaxy systems with limited brightest member galaxies been observed. Using four mock galaxy samples constructed from semi-analytical formation models, the subhalo abundance matching method and the conditional luminosity functions, respectively, we find that the luminosity gap between the brightest and the subsequent brightest member galaxies in a halo (group) can be used to significantly reduce the scatter in the halo mass estimation based on the luminosity of the brightest galaxy alone. Tests show that these corrections can significantly reduce the scatter in the halo mass estimations by $\sim 50\%$ to $\sim 70\%$ in massive halos depending on which member galaxies are considered. Comparing to the traditional ranking method, we find that this method works better for groups with less than five members, or in observations with very bright magnitude cut.


title: Standing spin wave mode in RFIM at T=0: Patterns and nonequilibrium phases

abstract: The dynamical responses of random field Ising model at zero temperature, driven by standing magnetic field wave, is studied by Monte Carlo simulation in two dimensions. The three different kinds of distribution of quenched random field are used here, uniform, bimodal and Gaussian. In all cases, three distinct dynamical phases were observed, namely, the pinned, structured and random. In the pinned phase no spin flip is observed. In the structured phase standing spin wave modes are observed. The random phase is shown no regular pattern. For a fixed value of the amplitude of the standing magnetic field wave, in the region of small quenched field, the system remains in a pinned phase. In the intermediate range of values of random field, a standing spin wave mode (structured phase) is observed. The regular pattern of this spin wave mode disappears for higher values of random field yielding a random phase. The comprehensive phase baundaries are drawn in all three cases. The boundary of pinned phase are analytically calculated for uniform and bimodal types of quenched random fields.


title: Temperature compensation via cooperative stability in protein degradation

abstract: Temperature compensation is a notable property of circadian oscillators that indicates the insensitivity of the oscillator system's period to temperature changes; the underlying mechanism, however, is still unclear. We investigated the influence of protein dimerization and cooperative stability in protein degradation on the temperature compensation ability of two oscillators. Here, cooperative stability means that high-order oligomers are more stable than their monomeric counterparts. The period of an oscillator is affected by the parameters of the dynamic system, which in turn are influenced by temperature. We adopted the Repressilator and the Atkinson oscillator to analyze the temperature sensitivity of their periods. Phase sensitivity analysis was employed to evaluate the period variations of different models induced by perturbations to the parameters. Furthermore, we used experimental data provided by other studies to determine the reasonable range of parameter temperature sensitivity. We then applied the linear programming method to the oscillatory systems to analyze the effects of protein dimerization and cooperative stability on the temperature sensitivity of their periods, which reflects the ability of temperature compensation in circadian rhythms. Our study explains the temperature compensation mechanism for circadian clocks. Compared with the no-dimer mathematical model and linear model for protein degradation, our theoretical results show that the nonlinear protein degradation caused by cooperative stability is more beneficial for realizing temperature compensation of the circadian clock.


title: Dynamics of flat actions on totally disconnected, locally compact groups

abstract: Let $G$ be a totally disconnected, locally compact group and let $H$ be a virtually flat (for example, polycyclic) group of automorphisms of $G$. We study the structure of, and relationships between, various subgroups of $G$ defined by the dynamics of $H$. In particular, we consider the following four subgroups: the intersection of all tidy subgroups for $H$ on $G$ (in the case that $H$ is flat); the intersection of all $H$-invariant open subgroups of $G$; the smallest closed $H$-invariant subgroup $D$ such that no $H$-orbit on $G/D$ accumulates at the trivial coset; and the group generated by the closures of contraction groups of elements of $H$ on $G$.


title: Some Results on the Regularization of LSQR for Large-Scale Discrete Ill-Posed Problems

abstract: LSQR, a Lanczos bidiagonalization based Krylov subspace iterative method, and its mathematically equivalent CGLS applied to normal equations system, are commonly used for large-scale discrete ill-posed problems. It is well known that LSQR and CGLS have regularizing effects, where the number of iterations plays the role of the regularization parameter. However, it has long been unknown whether the regularizing effects are good enough to find best possible regularized solutions. Here a best possible regularized solution means that it is at least as accurate as the best regularized solution obtained by the truncated singular value decomposition (TSVD) method. In this paper, we establish bounds for the distance between the $k$-dimensional Krylov subspace and the $k$-dimensional dominant right singular space. They show that the Krylov subspace captures the dominant right singular space better for severely and moderately ill-posed problems than for mildly ill-posed problems. Our general conclusions are that LSQR has better regularizing effects for the first two kinds of problems than for the third kind, and a hybrid LSQR with additional regularization is generally needed for mildly ill-posed problems. Exploiting the established bounds, we derive an estimate for the accuracy of the rank $k$ approximation generated by Lanczos bidiagonalization. Numerical experiments illustrate that the regularizing effects of LSQR are good enough to compute best possible regularized solutions for severely and moderately ill-posed problems, stronger than our theory predicts, but they are not for mildly ill-posed problems and additional regularization is needed.


title: Hyperbolic geometry in the work of Johann Heinrich Lambert

abstract: The memoir Theorie der Parallellinien (1766) by Johann Heinrich Lambert is one of the founding texts of hyperbolic geometry, even though its author's aim was, like many of his pre-decessors', to prove that such a geometry does not exist. In fact, Lambert developed his theory with the hope of finding a contradiction in a geometry where all the Euclidean axioms are kept except the parallel axiom and that the latter is replaced by its negation. In doing so, he obtained several fundamental results of hyperbolic geometry. This was sixty years before the first writings of Lobachevsky and Bolyai appeared in print. In the present paper, we present Lambert's main results and we comment on them. A French translation of the Theorie der Parallellinien, together with an extensive commentary, has just appeared in print (A. Papadopoulos and G. Th{\'e}ret, La th{\'e}orie des lignes parall{\`e}les de Johann Heinrich Lambert. Collection Sciences dans l'Histoire, Librairie Scientifique et Technique Albert Blanchard, Paris, 2014).