title: Spatially Resolved Magnetic Field Structure in the Disk of a T Tauri Star

abstract: Magnetic fields in accretion disks play a dominant role during the star formation process but have hitherto been observationally poorly constrained. Field strengths have been inferred on T Tauri stars themselves and possibly in the innermost part of the accretion disk, but the strength and morphology of the field in the bulk of the disk have not been observed. Unresolved measurements of polarized emission (arising from elongated dust grains aligned perpendicular to the field) imply average fields aligned with the disks. Theoretically, the fields are expected to be largely toroidal, poloidal, or a mixture of the two, which imply different mechanisms for transporting angular momentum in the disks of actively accreting young stars such as HL Tau. Here we report resolved measurements of the polarized 1.25 mm continuum emission from HL Tau's disk. The magnetic field on a scale of 80 AU is coincident with the major axis (~210 AU diameter) of the disk. From this we conclude that the magnetic field inside the disk at this scale cannot be dominated by a vertical component, though a purely toroidal field does not fit the data well either. The unexpected morphology suggests that the magnetic field's role for the accretion of a T Tauri star is more complex than the current theoretical understanding.


title: Planetary population synthesis coupled with atmospheric escape: a statistical view of evaporation

abstract: We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on a core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 $R_{\oplus}$. We find that evaporation can lead to a bimodal distribution of planetary sizes (Owen & Wu 2013) and to an "evaporation valley" running diagonally downwards in the orbital distance - planetary radius plane, separating bare cores from low-mass planet that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100\% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of $Kepler$ candidates.


title: Directly imaging damped Ly-alpha galaxies at z>2. III: The star formation rates of neutral gas reservoirs at z~2.7

abstract: We present results from a survey designed to probe the star formation properties of 32 damped Ly-alpha systems (DLAs) at z~2.7. By using the "double-DLA" technique that eliminates the glare of the bright background quasars, we directly measure the rest-frame FUV flux from DLAs and their neighbouring galaxies. At the position of the absorbing gas, we place stringent constraints on the unobscured star formation rates (SFRs) of DLAs to 2-sigma limits of <0.09-0.27 M/yr, corresponding to SFR surface densities <10^(-2.6)-10^(-1.5) M/yr/kpc^2. The implications of these limits for the star formation law, metal enrichment, and cooling rates of DLAs are examined. By studying the distribution of impact parameters as a function of SFRs for all the galaxies detected around these DLAs, we place new direct constraints on the bright end of the UV luminosity function of DLA hosts. We find that <13% of the hosts have SFR>2 M/yr at impact parameters b < (SFR/(M/yr))^(0.8)+6 kpc, differently from current samples of confirmed DLA galaxies. Our observations also disfavor a scenario in which the majority of DLAs arise from bright LBGs at distances 20 < b < 100 kpc. These new findings corroborate a picture in which DLAs do not originate from highly star forming systems that are coincident with the absorbers, and instead suggest that DLAs are associated with faint, possibly isolated, star-forming galaxies. Potential shortcomings of this scenario and future strategies for further investigation are discussed.


title: Pre-LHB Evolution of the Earth's Obliquity

abstract: The Earth's obliquity is stabilized by the Moon, which facilitates a rapid precession of the Earth's spin-axis, de-tuning the system away from resonance with orbital modulation. It is however, likely that the architecture of the Solar System underwent a dynamical instability-driven transformation, where the primordial configuration was more compact. Hence, the characteristic frequencies associated with orbital perturbations were likely faster in the past, potentially allowing for secular resonant encounters. In this work we examine if at any point in the Earth's evolutionary history, the obliquity varied significantly. Our calculations suggest that even though the orbital perturbations were different, the system nevertheless avoided resonant encounters throughout its evolution. This indicates that the Earth obtained its current obliquity during the formation of the Moon.


title: Quantum distillation and confinement of vacancies in a doublon sea

abstract: Ultracold atomic gases have revolutionized the study of non-equilibrium dynamics in quantum many-body systems. Many counterintuitive non-equilibrium effects have been observed, such as suppressed thermalization in a one-dimensional (1D) gas, the formation of repulsive self-bound dimers, and identical behaviors for attractive and repulsive interactions. Here, we observe the expansion of a bundle of ultracold 1D Bose gases in a flat-bottomed optical lattice potential. By combining in situ measurements with photoassociation, we follow the spatial dynamics of singly, doubly, and triply occupied lattice sites. The system sheds interaction energy by dissolving some doublons and triplons. Some singlons quantum distill out of the doublon center, while others remain confined. Our Gutzwiller mean-field model captures these experimental features in a physically clear way. These experiments might be used to study thermalization in systems with particle losses or the evolution of quantum entanglement, or if applied to fermions, to prepare very low entropy states.


title: A Uniform History for Galaxy Evolution

abstract: Recent observations indicate a remarkable similarity in the properties of evolving galaxies at fixed mass and redshift, prompting us to consider the possibility that most galaxies may evolve with a common history encompassing star formation, quasar accretion, and eventual quiescence. We quantify this by defining a "synchronization timescale" for galaxies as a function of mass and redshift that characterizes the extent to which different galaxies of a common mass are evolving in the same manner at various cosmic epochs. We measure this synchronization timescale using 9 different star-forming galaxy observations from the literature and SDSS quasar observations spanning $0 < z \lesssim 6$. Surprisingly, this synchronization timescale is a constant, approximately 1.5 Gyr for all combinations of mass and time. We also find that the ratio between the stellar mass of galaxies turning off star formation and black hole mass of turnoff quasars is approximately 30:1, much lower than the 500:1 for quiescent galaxies at low redshift. As a result, we propose a model in which the star-forming "main sequence", analogous quasar behavior, and other observations form a galactic evolution "main sequence", in which star formation occurs earliest, followed by supermassive black hole accretion, and feedback between the two are dominated by deterministic rather than stochastic processes.


title: Entropy mode loops and cosmological correlations during perturbative reheating

abstract: Recently, it has been shown that during preheating the entropy modes circulating in the loops, which correspond to the inflaton decay products, meaningfully modify the cosmological correlation functions at superhorizon scales. In this paper, we determine the significance of the same effect when reheating occurs in the perturbative regime. In a typical two scalar field model, the magnitude of the loop corrections are shown to depend on several parameters like the background inflaton amplitude in the beginning of reheating, the inflaton decay rate and the inflaton mass. Although the loop contributions turn out to be small as compared to the preheating case, they still come out larger than the loop effects during inflation.


title: Crowded Cluster Cores: Algorithms for Deblending in Dark Energy Survey Images

abstract: Deep optical images are often crowded with overlapping objects. This is especially true in the cores of galaxy clusters, where images of dozens of galaxies may lie atop one another. Accurate measurements of cluster properties require deblending algorithms designed to automatically extract a list of individual objects and decide what fraction of the light in each pixel comes from each object. We present new software called the Gradient And INterpolation based deblender (GAIN) as a secondary deblender to improve deblending the images of cluster cores. This software relies on using image intensity gradient and using an image interpolation technique usually used to correct flawed terrestrial digital images. We test this software on Dark Energy Survey coadd images. GAIN helps extracting unbiased photometry measurement for blended sources. It also helps improving detection completeness while introducing only a modest amount of spurious detections. For example, when applied to deep images simulated with high level of deblending difficulties, this software improves detection completeness from 91% to 97% for sources above the 10? limiting magnitude at 25.3 mag. We expect this software to be a useful tool for cluster population measurements.


title: On the Frequency of Potential Venus Analogs from Kepler Data

abstract: The field of exoplanetary science has seen a dramatic improvement in sensitivity to terrestrial planets over recent years. Such discoveries have been a key feature of results from the {\it Kepler} mission which utilizes the transit method to determine the size of the planet. These discoveries have resulted in a corresponding interest in the topic of the Habitable Zone (HZ) and the search for potential Earth analogs. Within the Solar System, there is a clear dichotomy between Venus and Earth in terms of atmospheric evolution, likely the result of the large difference ($\sim$ factor of two) in incident flux from the Sun. Since Venus is 95\% of the Earth's radius in size, it is impossible to distinguish between these two planets based only on size. In this paper we discuss planetary insolation in the context of atmospheric erosion and runaway greenhouse limits for planets similar to Venus. We define a ``Venus Zone'' (VZ) in which the planet is more likely to be a Venus analog rather than an Earth analog. We identify 43 potential Venus analogs with an occurrence rate ($\eta_{\venus}$) of $0.32^{+0.05}_{-0.07}$ and $0.45^{+0.06}_{-0.09}$ for M dwarfs and GK dwarfs respectively.


title: The diversity of quasars unified by accretion and orientation

abstract: Quasars are rapidly accreting supermassive black holes at the center of massive galaxies. They display a broad range of properties across all wavelengths, reflecting the diversity in the physical conditions of the regions close to the central engine. These properties, however, are not random, but form well-defined trends. The dominant trend is known as Eigenvector 1, where many properties correlate with the strength of optical iron and [OIII] emission. The main physical driver of Eigenvector 1 has long been suspected to be the quasar luminosity normalized by the mass of the hole (the Eddington ratio), an important quantity of the black hole accretion process. But a definitive proof has been missing. Here we report an analysis of archival data that reveals that Eddington ratio indeed drives Eigenvector 1. We also find that orientation plays a significant role in determining the observed kinematics of the gas, implying a flattened, disklike geometry for the fast-moving clouds close to the hole. Our results show that most of the diversity of quasar phenomenology can be unified with two simple quantities, Eddington ratio and orientation.


title: Extended M{\o}ller-Plesset perturbation theory for dynamical and static correlations

abstract: We present a novel method that appropriately handles both dynamical and static electron correlation in a balanced manner, using a perturbation theory on a spin-extended Hartree-Fock (EHF) wave function reference. While EHF is a suitable candidate for degenerate systems where static correlation is ubiquitous, it is known that most of dynamical correlation is neglected in EHF. In this work, we derive a perturbative correction to a fully spin-projected self-consistent wave function based on second-order M{\o}ller-Plesset perturbation theory (MP2). The proposed method efficiently captures the ability of EHF to describe static correlation in degeneracy, combined with MP2's ability to treat dynamical correlation effects. We demonstrate drastic improvements on molecular ground state and excited state potential energy curves and singlet-triplet splitting energies over both EHF and MP2 with similar computational effort to the latter.


title: Time-dependent potential barriers and superarrivals

abstract: Scattering of a Gaussian wavepacket from rectangular potential barriers with increasing widths or heights is studied numerically. It is seen that during a certain time interval the time-evolving transmission probability increases compared to the corresponding unperturbed cases. In the literature this effect is known as superarrival in transmission probability. We present a trajectory-based explanation for this effect by using the concept of quantum potential energy and computing a selection of Bohmian trajectories. Relevant parameters in superarrivals are determined for the case that the barrier width increases linearly during the dispersion of the wavepacket. Nonlinear in time perturbation is also considered.


title: Fixing the Reference Frame for PPMXL Proper Motions Using Extragalactic Sources

abstract: We quantify and correct systematic errors in PPMXL proper motions using extragalactic sources from the first two LAMOST data releases and the Veron-Cetty & Veron Catalog of Quasars. Although the majority of the sources are from the Veron catalog, LAMOST makes important contributions in regions that are not well-sampled by previous catalogs, particularly at low Galactic latitudes and in the south Galactic cap. We show that quasars in PPMXL have measureable and significant proper motions, which reflect the systematic zero-point offsets present in the catalog. We confirm the global proper motion shifts seen by Wu, Ma, & Zhou (2011), and additionally find smaller-scale fluctuations of the QSO-derived corrections to an absolute frame. We average the proper motions of 158,106 extragalactic objects in bins of 3x3 degrees and present a table of proper motion corrections.


title: Weak Values and Modular Variables From a Quantum Phase Space Perspective

abstract: We address two major conceptual developments introduced by Aharonov and collaborators through a \textit{quantum phase space} approach: the concept of \textit{modular variables} devised to explain the phenomena of quantum dynamical non-locality and the \textit{two-state formalism} for Quantum Mechanics which is a retrocausal time-symmetric interpretation of quantum physics which led to the discovery of \textit{weak values.} We propose that a quantum phase space structure underlies these profound physical insights in a unifying manner. For this, we briefly review the Weyl-Wigner and the coherent state formalisms as well as the inherent symplectic structures of quantum projective spaces in order to gain a deeper understanding of the weak value concept. We also review Schwinger's finite quantum kinematics so that we may apply this discrete formalism to understand Aharonov's modular variable concept in a different manner that has been proposed before in the literature. We discuss why we believe that this\ is indeed the correct kinematic framework for the modular variable concept and how this may shine some light on the physical distinction between quantum dynamical non-locality and the kinematic non-locality, generally associated with entangled quantum systems.


title: Storage and retrieval of ultrafast single photons using a room-temperature diamond quantum memory

abstract: We report the storage and retrieval of single photons, via a quantum memory, in the optical phonons of room-temperature bulk diamond. The THz-bandwidth heralded photons are generated by spontaneous parametric downconversion and mapped to phonons via a Raman transition, stored for a variable delay, and released on demand. The second-order correlation of the memory output is $g^{(2)}(0) = 0.65 \pm 0.07$, demonstrating preservation of non-classical photon statistics throughout storage and retrieval. The memory is low-noise, high-speed and broadly tunable; it therefore promises to be a versatile light-matter interface for local quantum processing applications.


title: Simplified Models for Dark Matter and Missing Energy Searches at the LHC

abstract: The study of collision events with missing energy as searches for the dark matter (DM) component of the Universe are an essential part of the extensive program looking for new physics at the LHC. Given the unknown nature of DM, the interpretation of such searches should be made broad and inclusive. This report reviews the usage of simplified models in the interpretation of missing energy searches. We begin with a brief discussion of the utility and limitation of the effective field theory approach to this problem. The bulk of the report is then devoted to several different simplified models and their signatures, including s-channel and t-channel processes. A common feature of simplified models for DM is the presence of additional particles that mediate the interactions between the Standard Model and the particle that makes up DM. We consider these in detail and emphasize the importance of their inclusion as final states in any coherent interpretation. We also review some of the experimental progress in the field, new signatures, and other aspects of the searches themselves. We conclude with comments and recommendations regarding the use of simplified models in Run-II of the LHC.


title: Transparent subdiffraction optics: nanoscale light confinement without metal

abstract: The integration of nanoscale electronics with conventional optical devices is restricted by the diffraction limit of light. Metals can confine light at the subwavelength scales needed, but they are lossy, while dielectric materials do not confine evanescent waves outside a waveguide or resonator, leading to cross talk between components. We introduce a paradigm shift in light confinement strategy and show that light can be confined below the diffraction limit using completely transparent artificial media. Our approach relies on controlling the optical momentum of evanescent waves, an important electromagnetic property overlooked in photonic devices. For practical applications, we propose a class of waveguides using this approach that outperforms the cross talk performance by 1 order of magnitude as compared to any existing photonic structure. Our work overcomes a critical stumbling block for nanophotonics by completely averting the use of metals and can impact electromagnetic devices from the visible to microwave frequency ranges.


title: Blondel et les oscillations auto-entretenues

abstract: In 1893, the "physicist-engineer" Andr\'e Blondel invents the oscilloscope for displaying voltage and current variables. With this powerful means of investigation, he first studies the phenomena of the arc then used for the coastal and urban lighting and then, the singing arc used as a transmitter of radio waves in wireless telegraphy. In 1905, he highlights a new type of non-sinusoidal oscillations in the singing arc. Twenty years later, Balthasar van der Pol will recognize that such oscillations were in fact "relaxation oscillations". To explain this phenomenon, he uses a representation in the phase plane and shows that its evolution takes the form of small cycles. During World War I the triode gradually replaces the singing arc in transmission systems. At the end of the war, using analogy, Blondel transposes to the triode most of the results he had obtained for the singing arc. In April 1919, he publishes a long memoir in which he introduces the terminology "self-sustained oscillations" and proposes to illustrate this concept starting from the example of the Tantalus cup which is then picked up by Van der Pol and Philippe Le Corbeiller. He then provides the definition of a self sustained system which is quite similar to that given later by Aleksandr Andronov and Van der Pol. To study the stability of oscillations sustained by the triode and by the singing arc he uses, this time, a representation in the complex plane and he expresses the amplitude in polar coordinates. He then justifies the maintaining of oscillations by the existence cycles which nearly present all the features of Poincar\'e's limit cycles. Finally, in November 1919, Blondel performs, a year before Van der Pol, the setting in equation of the triode oscillations. In March 1926, Blondel establishes the differential equation characterizing the oscillations of the singing arc, completely similar to that obtained concomitantly by Van der Pol for the triode. Thus, throughout his career, Blondel, has made fundamental and relatively unknown contributions to the development of the theory of nonlinear oscillations. The purpose of this article is to analyze his main work in this area and to measure their importance or influence by placing them in the perspective of the development of this theory.


title: On the evolution of rotating accreting white dwarfs and type Ia supernovae

abstract: The potential importance of the angular momentum which is gained by accreting white dwarfs (WDs) has been increasingly recognized in the context of type Ia supernova (SN Ia) single-degenerate model. The expectation that the spin of the WD can delay the explosion should help the single-degenerate model to be consistent with the observed properties of most SNe Ia, in particular by avoiding hydrogen contamination. In this article, we attempt to study the most prominent single-degenerate supersoft (WD + MS) channel when the rotation of accreting WDs is considered. We present a detailed binary population synthesis study to examine the predicted population of SNe Ia for this channel. For our standard model, we find that 77% of these SNe Ia explode with WD masses which are low enough to be supported by solid-body rotation (1.378-1.5Msun); this is a substantially higher proportion than found by previous work. Only 2% have WD explosion masses >2.0Msun; these require the initial WD mass to be larger than 1.0 Msun. We further discuss the possible origin of the diversity of SNe Ia from the pre- and post- accretion properties of the WDs in this population. We also suggest that some SN Ia progenitors with substantial circumstellar hydrogen, including some apparent type IIn SNe, might be related to WDs which required support from differential rotation to avoid explosion, since these can still be accreting from hydrogen-rich donors with a relatively high mass-transfer rate at the time of the SN explosion.


title: Co-adaptation in a Handwriting Recognition System

abstract: Handwriting is a natural and versatile method for human-computer interaction, especially on small mobile devices such as smart phones. However, as handwriting varies significantly from person to person, it is difficult to design handwriting recognizers that perform well for all users. A natural solution is to use machine learning to adapt the recognizer to the user. One complicating factor is that, as the computer adapts to the user, the user also adapts to the computer and probably changes their handwriting. This paper investigates the dynamics of co-adaptation, a process in which both the computer and the user are adapting their behaviors in order to improve the speed and accuracy of the communication through handwriting. We devised an information-theoretic framework for quantifying the efficiency of a handwriting system where the system includes both the user and the computer. Using this framework, we analyzed data collected from an adaptive handwriting recognition system and characterized the impact of machine adaptation and of human adaptation. We found that both machine adaptation and human adaptation have significant impact on the input rate and must be considered together in order to improve the efficiency of the system as a whole.


title: One-loop corrections to neutralino-stop coannihilation revisited

abstract: We discuss the ${\cal O}(\alpha_s)$ supersymmetric QCD corrections to neutralino-stop coannihilation into a top quark and a gluon in the Minimal Supersymmetric Standard Model (MSSM). This particular channel can be numerically important in wide ranges of the MSSM parameter space with rather light stops. We discuss technical details such as the renormalization scheme and the phase-space slicing method with two cutoffs. We also comment on improvements with respect to earlier works on the given process. Further, we study for the first time the phenomenologically very interesting interplay of neutralino-stop coannihilation with neutralino-pair annihilation into quark pairs taking the full next-to-leading order SUSY-QCD corrections into account. We demonstrate that the numerical impact of these corrections on the total (co)annihilation cross section and finally on the theoretically predicted neutralino relic density is significant.


title: Strong General Position

abstract: We say that a finite set S of points in R^d is in "strong general position" if for any collection {F_1,..., F_r} of r pairwise disjoint subsets of S (1 <= r <= |S|) we have: d-dim (the intersection of aff F_1,aff F_2,...,aff F_r) = min{d+1, (d-dim aff F_1)+(d-dim aff F_2)+...+(d-dim aff F_r)}. In this paper we reduce the set of conditions that one has to check in order to determine if S is in "strong general position".


title: Characterizing SL2S galaxy groups using the Einstein radius

abstract: We analyzed the Einstein radius, $\theta_E$, in our sample of SL2S galaxy groups, and compared it with $R_A$ (the distance from the arcs to the center of the lens), using three different approaches: 1.- the velocity dispersion obtained from weak lensing assuming a Singular Isothermal Sphere profile ($\theta_{E,I}$), 2.- a strong lensing analytical method ($\theta_{E,II}$) combined with a velocity dispersion-concentration relation derived from numerical simulations designed to mimic our group sample, 3.- strong lensing modeling ($\theta_{E,III}$) of eleven groups (with four new models presented in this work) using HST and CFHT images. Finally, $R_A$ was analyzed as a function of redshift $z$ to investigate possible correlations with L, N, and the richness-to-luminosity ratio (N/L). We found a correlation between $\theta_{E}$ and $R_A$, but with large scatter. We estimate $\theta_{E,I}$ = (2.2 $\pm$ 0.9) + (0.7 $\pm$ 0.2)$R_A$, $\theta_{E,II}$ = (0.4 $\pm$ 1.5) + (1.1 $\pm$ 0.4)$R_A$, and $\theta_{E,III}$ = (0.4 $\pm$ 1.5) + (0.9 $\pm$ 0.3)$R_A$ for each method respectively. We found a weak evidence of anti-correlation between $R_A$ and $z$, with Log$R_A$ = (0.58$\pm$0.06) - (0.04$\pm$0.1)$z$, suggesting a possible evolution of the Einstein radius with $z$, as reported previously by other authors. Our results also show that $R_A$ is correlated with L and N (more luminous and richer groups have greater $R_A$), and a possible correlation between $R_A$ and the N/L ratio. Our analysis indicates that $R_A$ is correlated with $\theta_E$ in our sample, making $R_A$ useful to characterize properties like L and N (and possible N/L) in galaxy groups. Additionally, we present evidence suggesting that the Einstein radius evolves with $z$.


title: Effects of self-consistency and plasmon-pole models on GW calculations for closed-shell molecules

abstract: We present theoretical calculations of quasiparticle energies in closed-shell molecules using the GW method. We compare three different approaches: a full-frequency $G_0W_0$ (FF-$G_0W_0$) method with density functional theory (DFT-PBE) used as a starting mean field; a full-frequency $GW_0$ (FF-$GW_0$) method where the interacting Green's function is approximated by replacing the DFT energies with self-consistent quasiparticle energies or Hartree-Fock energies; and a $G_0W_0$ method with a Hybertsen-Louie generalized plasmon-pole model (HL GPP-$G_0W_0$). While the latter two methods lead to good agreement with experimental ionization potentials and electron affinities for methane, ozone, and beryllium oxide molecules, FF-$G_0W_0$ results can differ by more than one electron volt from experiment. We trace this failure of the FF-$G_0W_0$ method to the occurrence of incorrect self-energy poles describing shake-up processes in the vicinity of the quasiparticle energies.


title: The Hildreth's Algorithm with Applications to Soft Constraints for User Interface Layout

abstract: The Hildreth's algorithm is a row action method for solving large systems of inequalities. This algorithm is efficient for problems with sparse matrices, as opposed to direct methods such as Gaussian elimination or QR-factorization. We apply the Hildreth's algorithm, as well as a randomized version, along with prioritized selection of the inequalities, to efficiently detect the highest priority feasible subsystem of equations. We prove convergence results and feasibility criteria for both cyclic and randomized Hildreth's algorithm, as well as a mixed algorithm which uses Hildreth's algorithm for inequalities and Kaczmarz algorithm for equalities. These prioritized, sparse systems of inequalities commonly appear in constraint-based user interface (UI) layout specifications. The performance and convergence of these proposed algorithms are evaluated empirically using randomly generated UI layout specifications of various sizes. The results show that these methods offer improvements in performance over standard methods like Matlab's LINPROG, a well-known efficient linear programming solver, and the recent developed Kaczmarz algorithm with prioritized IIS detection.


title: A structured matrix factorization framework for large scale calcium imaging data analysis

abstract: We present a structured matrix factorization approach to analyzing calcium imaging recordings of large neuronal ensembles. Our goal is to simultaneously identify the locations of the neurons, demix spatially overlapping components, and denoise and deconvolve the spiking activity of each neuron from the slow dynamics of the calcium indicator. The matrix factorization approach relies on the observation that the spatiotemporal fluorescence activity can be expressed as a product of two matrices: a spatial matrix that encodes the location of each neuron in the optical field and a temporal matrix that characterizes the calcium concentration of each neuron over time. We present a simple approach for estimating the dynamics of the calcium indicator as well as the observation noise statistics from the observed data. These parameters are then used to set up the matrix factorization problem in a constrained form that requires no further parameter tuning. We discuss initialization and post-processing techniques that enhance the performance of our method, along with efficient and largely parallelizable algorithms. We apply our method to {\it in vivo} large scale multi-neuronal imaging data and also demonstrate how similar methods can be used for the analysis of {\it in vivo} dendritic imaging data.


title: New analytic solution for the heat flow through a general harmonic network

abstract: We present a new analytic expression for the heat current through a general harmonic network coupled with Ohmic reservoirs. We use a new method that enables us to express the stationary state of the network in terms of the eigenvectors and eigenvalues of a generalized cubic eigenvalue problem. In this way, we obtain exact formulae for the heat current and the local temperature inside the network. Our method does not rely in the usual assumptions of weak coupling to the environments or on the existence of an infinite cutoff in the environmental spectral densities. We use this method to study non-equilibrium processes without the weak coupling and Markovian approximations. As a first application of our method, we revisit the problem of heat conduction in 2D and 3D crystals with binary mass disorder. We complement previous results showing that for small systems the scaling of the heat current with the system size greatly depends on the strength of the interaction between system and reservoirs. This somewhat counter-intuitive result seems not to have been noticed before.


title: Non-Convex Boosting Overcomes Random Label Noise

abstract: The sensitivity of Adaboost to random label noise is a well-studied problem. LogitBoost, BrownBoost and RobustBoost are boosting algorithms claimed to be less sensitive to noise than AdaBoost. We present the results of experiments evaluating these algorithms on both synthetic and real datasets. We compare the performance on each of datasets when the labels are corrupted by different levels of independent label noise. In presence of random label noise, we found that BrownBoost and RobustBoost perform significantly better than AdaBoost and LogitBoost, while the difference between each pair of algorithms is insignificant. We provide an explanation for the difference based on the margin distributions of the algorithms.


title: Multicomponent solution solidification with arrested phase separation model for liquid-to-glass transition

abstract: Due to nonuniform aggregation in liquid state, from the thermodynamic point of view any glass-forming liquid in the vicinity of the liquid-to-solid phase transition temperature, irrespective of its actual chemical composition, shall be described in terms of a complex multicomponent solution whose comprised of the same chemical elements components have characteristic atomic arrangement deviating to various extent from the thermodynamic ground state with respect to the size, shape, density, structure, and stoichiometry. Therefore, glass transition appears to be a process of non-equilibrium solidification of multicomponent solution upon its rapid cooling. The essential feature of this process is that the attempts of the liquid and solid phases of the solidifying solution to separate out are largely arrested due to quenching. Thus, the solidification occurs in the absence of solid-liquid interface, so the substance in the liquid-to-glass transition region is observed behaving like fluid with rapidly growing viscosity that reflects the formation of mechanically rigid and stable bound configurations. It is shown that glass transition shall be classified as phase transition in multicomponent solutions and not a standalone phenomenon.


title: Gradual Diffusive Capture: Slow Death by Many Mosquito Bites

abstract: We study the dynamics of a single diffusing particle (a "man") with diffusivity $D_M$ that is attacked by another diffusing particle (a "mosquito") with fixed diffusivity $D_m$. Each time the mosquito meets and bites the man, the diffusivity of the man is reduced by a fixed amount, while the diffusivity of the mosquito is unchanged. The mosquito is also displaced by a small distance $\pm a$ with respect to the man after each encounter. The man is defined as dead when $D_M$ reaches zero. At the moment when the man dies, his probability distribution of displacements $x$ is given by a Cauchy form, which asymptotically decays as $x^{-2}$, while the distribution of times $t$ when the man dies asymptotically decays as $t^{-3/2}$, which has the same form as the one-dimensional first-passage probability.


title: A Framework for Practical Parallel Fast Matrix Multiplication

abstract: Matrix multiplication is a fundamental computation in many scientific disciplines. In this paper, we show that novel fast matrix multiplication algorithms can significantly outperform vendor implementations of the classical algorithm and Strassen's fast algorithm on modest problem sizes and shapes. Furthermore, we show that the best choice of fast algorithm depends not only on the size of the matrices but also the shape. We develop a code generation tool to automatically implement multiple sequential and shared-memory parallel variants of each fast algorithm, including our novel parallelization scheme. This allows us to rapidly benchmark over 20 fast algorithms on several problem sizes. Furthermore, we discuss a number of practical implementation issues for these algorithms on shared-memory machines that can direct further research on making fast algorithms practical.


title: On the Reciprocal of the Binary Generating Function for the Sum of Divisors

abstract: If \(A \) is a set of natural numbers containing \(0 \), then there is a unique nonempty "reciprocal" set \(B \) of natural numbers (containing \(0 \)) such that every positive integer can be written in the form \(a + b \), where \(a \in A \) and \(b \in B \), in an even number of ways. Furthermore, the generating functions for \(A \) and \(B \) over \(\FF_2 \) are reciprocals in \(\FF_2 [[q]] \). We consider the reciprocal set \(B \) for the set \(A \) containing \(0 \) and all integers such that \(\sigma(n) \) is odd, where \(\sigma(n) \) is the sum of all the positive divisors of \(n \). This problem is motivated by Euler's "Pentagonal Number Theorem", a corollary of which is that the set of natural numbers \(n \) so that the number \(p(n) \) of partitions of an integer \(n \) is odd is the reciprocal of the set of generalized pentagonal numbers (integers of the form \(k(3k\pm1)/2 \), where \(k \) is a natural number). An old (1967) conjecture of Parkin and Shanks is that the density of integers \(n \) so that \(p(n) \) is odd (equivalently, even) is \(1/2 \). Euler also found that \(\sigma(n) \) satisfies an almost identical recurrence as that given by the Pentagonal Number Theorem, so we hope to shed light on the Parkin-Shanks conjecture by computing the density of the reciprocal of the set containing the natural numbers with \(\sigma(n) \) odd (\(\sigma(0)=1 \) by convention). We conjecture this particular density is \(1/32 \) and prove that it lies between \(0 \) and \(1/16 \). We finish with a few surprising connections between certain Beatty sequences and the sequence of integers \(n \) for which \(\sigma(n) \) is odd.


title: Fullerene photoemission time delay explores molecular cavity in attoseconds

abstract: Photoelectron spectroscopy earlier probed oscillations in C60 valence emissions, producing series of minima whose energy separation depends on the molecular cavity. We show here that the quantum phase at these cavity minima exhibits variations from strong electron correlations in C60, causing rich structures in the emission time delay. Hence, these minima offer unique spectral zones to directly explore multielectron forces via attosecond RABITT interferometry not only in fullerenes, but also in clusters and nanostructures for which such minima are likely abundant.


title: The Hirzebruch $\chi_y$-genus and Poincare polynomial revisited

abstract: The Hirzebruch $\chi_y$-genus and Poincare polynomial share some similar features. In this article we investigate two of their similar features simultaneously. Through this process we shall derive several new results as well as reprove and improve some known results.


title: -1-Phenomena for the pluri $\chi_y$-genus and elliptic genus

abstract: Several independent articles have observed that the Hirzebruch $\chi_y$-genus has an important feature, which the author calls -1-phenomenon and tells us that the coefficients of the Taylor expansion of the $\chi_y$-genus at $y=-1$ have explicit expressions. Hirzebruch's original $\chi_y$-genus can be extended towards two directions: the pluri-case and the case of elliptic genus. This paper contains two parts in which we investigate the -1-phenomena in these two generalized cases respectively and show that in each case there exists a -1-phenomenon in a suitable sense. Our main results in the first part have an application, which states that all characteristic numbers (Chern numbers and Pontrjagin numbers) on manifolds can be expressed, in a very explicit way, in terms of some rationally linear combination of indices of some elliptic operators. This gives an analytic interpretation of characteristic numbers and affirmatively answers a question posed by the author several years ago. The second part contains our attempt to generalize this -1-phenomenon to elliptic genus, a modern version of the $\chi_y$-genus. We first extend the elliptic genus of an almost-complex manifold to a twisted version where an extra complex vector bundle is involved, and show that it is a weak Jacobi form under some assumptions. A suitable manipulation on the theory of Jacobi form will produce new modular forms from this weak Jacobi form and thus much arithmetic information related to the underlying manifold can be obtained, in which the -1-phenomenon of the original $\chi_y$-genus is hidden.


title: Efficient Algorithms for Discrepancy Minimization in Convex Sets

abstract: A result of Spencer states that every collection of $n$ sets over a universe of size $n$ has a coloring of the ground set with $\{-1,+1\}$ of discrepancy $O(\sqrt{n})$. A geometric generalization of this result was given by Gluskin (see also Giannopoulos) who showed that every symmetric convex body $K\subseteq R^n$ with Gaussian measure at least $e^{-\epsilon n}$, for a small $\epsilon>0$, contains a point $y\in K$ where a constant fraction of coordinates of $y$ are in $\{-1,1\}$. This is often called a partial coloring result. While both these results were inherently non-algorithmic, recently Bansal (see also Lovett-Meka) gave a polynomial time algorithm for Spencer's setting and Rothvo\ss gave a randomized polynomial time algorithm obtaining the same guarantee as the result of Gluskin and Giannopoulos. This paper has several related results. First we prove another constructive version of the result of Gluskin and Giannopoulos via an optimization of a linear function. This implies a linear programming based algorithm for combinatorial discrepancy obtaining the same result as Spencer. Our second result gives a new approach to obtains partial colorings and shows that every convex body $K\subseteq R^n$, possibly non-symmetric, with Gaussian measure at least $e^{-\epsilon n}$, for a small $\epsilon>0$, contains a point $y\in K$ where a constant fraction of coordinates of $y$ are in $\{-1,1\}$. Finally, we give a simple proof that shows that for any $\delta >0$ there exists a constant $c>0$ such that given a body $K$ with $\gamma_n(K)\geq \delta$, a uniformly random $x$ from $\{-1,1\}^n$ is in $cK$ with constant probability. This gives an algorithmic version of a special case of the result of Banaszczyk.


title: Eigenvalues of vector fields, Bott's residue formula and integral invariants

abstract: Given a compatible vector field on a compact connected almost-complex manifold, we show in this article that the multiplicities of eigenvalues among the zero point set of this vector field have intimate relations. We highlight a special case of our result and reinterpret it as a vanishing-type result in the framework of the celebrated Atiyah-Bott-Singer localization formula. This new point of view, via the Chern-Weil theory and a strengthened version of Bott's residue formula observed by Futaki and Morita, can lead to an obstruction to Killing real holomorphic vector fields on compact Hermitian manifolds in terms of a curvature integral.


title: Two dimensional kicked quantum Ising model: dynamical phase transitions

abstract: Using an efficient one and two qubit gate simulator, operating on graphical processing units, we investigate ergodic properties of a quantum Ising spin 1/2 model on a two dimensional lattice, which is periodically driven by a $\delta$-pulsed transverse magnetic field. We consider three different dynamical properties: (i) level density and (ii) level spacing distribution of the Floquet quasienergy spectrum, as well as (iii) time-averaged autocorrelation function of components of the magnetization. Varying the parameters of the model, we found transitions between ordered (non ergodic) and quantum chaotic (ergodic) phases, but the transitions between flat and non-flat spectral density {\em do not} correspond to transitions between ergodic and non-ergodic local observables. Even more surprisingly, we found nice agreement of level spacing distribution with the Wigner surmise of random matrix theory for almost all values of parameters except where the model is essentially noninteracting, even in the regions where local observables are not ergodic or where spectral density is non-flat. These findings put in question the versatility of the interpretation of level spacing distribution in many-body systems and stress the importance of the concept of locality.


title: Non-additive simple potentials for pre-programmed self-assembly

abstract: A major goal in nanoscience and nanotechnology is the self-assembly of any desired complex structure with a system of particles interacting through simple potentials. To achieve this objective, intense experimental and theoretical efforts are currently concentrated in the development of the so called "patchy" particles. Here we follow a completely different approach and introduce a very accessible model to produce a large variety of pre-programmed two-dimensional (2D) complex structures. Our model consists of a binary mixture of particles that interact through isotropic in plane interactions that is able to self-assemble into targeted lattices by the appropriate choice of a small number of geometrical parameters and interaction strengths. We study the system using Monte Carlo computer simulations and, despite its simplicity, we are able to self assemble potentially useful structures such as chains, stripes, Kagome, twisted Kagome, honeycomb, square, Archimedean and quasicrystaline tilings. Our model is designed such that it may be imediately implemented in experiments using existing techniques to build particles with different shapes and interactions. Thus, it represents a promising strategy for bottom-up nano-fabrication.


title: Valence-shell photoionization of Ag-like Xe$^{7+}$ ions : experiment and theory

abstract: We report on experimental and theoretical results for the photoionization of Ag-like xenon ions, Xe$^{7+}$, in the photon energy range 95 to 145~eV. The measurements were carried out at the Advanced Light Source at an energy resolution of $\Delta$E = 65 meV with additional measurements made at $\Delta$E = 28 meV and 39 meV. Small resonance features below the ground-state ionization threshold, at about 106 eV, are due to the presence of metastable Xe$^{7+} (4d^{10} 4f~^2{\rm F}^{\circ}_{5/2,7/2})$ ions in the ion beam. On the basis of the accompanying theoretical calculations using the Dirac Atomic R-matrix Codes (DARC), an admixture of only a few percent of metastable ions in the parent ion beam is inferred, with almost 100\% of the parent ions in the $(4d^{10}5s ~^2{\rm S_{1/2}})$ ground level. The cross-section is dominated by a very strong resonance associated with $4d \rightarrow 5f$ excitation and subsequent autoionization. This prominent feature in the measured spectrum is the $4d^95s5f ~^2{\rm P}^{\circ}$ resonance located at (122.139 $\pm$ 0.01)~eV. An absolute peak cross-section of 1.2 Gigabarns was measured at 38 meV energy resolution. The experimental natural width $\Gamma$ = 76 $\pm$ 3 meV of this resonance compares well with the theoretical estimate of 88 meV obtained from the DARC calculation with 249 target states. Given the complexity of the system, overall satisfactory agreement between theory and experiment is obtained for the photon energy region investigated.


title: Quantum Edge Detection for Image Segmentation in Optical Environments

abstract: A quantum edge detector for image segmentation in optical environments is presented in this work. A Boolean version of the same detector is presented too. The quantum version of the new edge detector works with computational basis states, exclusively. This way, we can easily avoid the problem of quantum measurement retrieving the result of applying the new detector on the image. Besides, a new criterion and logic based on projections onto vertical axis of Bloch's Sphere exclusively are presented too. This approach will allow us: 1) a simpler development of logic quantum operations, where they will closer to those used in the classical logic operations, 2) building simple and robust classical-to-quantum and quantum-to-classical interfaces. Said so far is extended to quantum algorithms outside image processing too. In a special section on metric and simulations, a new metric based on the comparison between the classical and quantum versions algorithms for edge detection of images is presented. Notable differences between the results of classical and quantum versions of such algorithms (outside and inside of quantum computer, respectively) show the existence of implementation problems involved in the experiment, and that they have not been properly modeled for optical environments. However, although they are different, the quantum results are equally valid. The latter is clearly seen in the computer simulations


title: Intrinsic scales for high-dimensional LEVY-driven models with non-Markovian synchronizing updates

abstract: We propose stochastic $N$-component synchronization models $(x_{1}(t),...,x_{N}(t))$, $x_{j}\in\mathbb{R}^{d}$, $t\in\mathbb{R}_{+}$, whose dynamics is described by Levy processes and synchronizing jumps. We prove that symmetric models reach synchronization in a stochastic sense: differences between components $d_{kj}^{(N)}(t)=x_{k}(t)-x_{j}(t)$ have limits in distribution as $t\rightarrow\infty$. We give conditions of existence of natural (intrinsic) space scales for large synchronized systems, i.e., we are looking for such sequences $\{b_{N}\}$ that distribution of $d_{kj}^{(N)}(\infty)/b_{N}$ converges to some limit as $N\rightarrow\infty$. It appears that such sequence exists if the Levy process enters a domain of attraction of some stable law. For Markovian synchronization models based on $\alpha$-stable Levy processes this results holds for any finite $N$ in the precise form with $b_{N}=(N-1)^{1/\alpha}$. For non-Markovian models similar results hold only in the asymptotic sense. The class of limiting laws includes the Linnik distributions. We also discuss generalizations of these theorems to the case of non-uniform matrix-based intrinsic scales. The central point of our proofs is a representation of characteristic functions of $d_{kj}^{(N)}(t)$ via probability distribution of a superposition of $N$ independent renewal processes.


title: Crystal structure on rigged configurations and the filling map

abstract: In this paper, we extend work of the first author on a crystal structure on rigged configurations of simply-laced type to all non-exceptional affine types using the technology of virtual rigged configurations and crystals. Under the bijection between rigged configurations and tensor products of Kirillov-Reshetikhin crystals specialized to a single tensor factor, we obtain a new tableaux model for Kirillov-Reshetikhin crystals. This is related to the model in terms of Kashiwara-Nakashima tableaux via a filling map, generalizing the recently discovered filling map in type $D_n^{(1)}$.


title: Local Hamiltonians for quantitative Green's function embedding methods

abstract: Embedding calculations that find approximate solutions to the Schr\"{o}dinger equation for large molecules and realistic solids are performed commonly in a three step procedure involving (i) construction of a model system with effective interactions approximating the low energy physics of the initial realistic system, (ii) mapping the model system onto an impurity Hamiltonian, and (iii) solving the impurity problem. We have developed a novel procedure for parametrizing the impurity Hamiltonian that avoids the mathematically uncontrolled step of constructing the low energy model system. Instead, the impurity Hamiltonian is immediately parametrized to recover the self-energy of the realistic system in the limit of high frequencies or short time. The effective interactions parametrizing the fictitious impurity Hamiltonian are local to the embedded regions, and include all the non-local interactions present in the original realistic Hamiltonian in an implicit way. We show that this impurity Hamiltonian can lead to excellent total energies and self-energies that approximate the quantities of the initial realistic system very well. Moreover, we show that as long as the effective impurity Hamiltonian parametrization is designed to recover the self-energy of the initial realistic system for high frequencies, we can expect a good total energy and self-energy. Finally, we propose two practical ways of evaluating effective integrals for parametrizing impurity models.


title: Trees, ladders and graphs

abstract: We introduce a new method to construct uncountably chromatic graphs from non special trees and ladder systems. Answering a question of P. Erd\H{o}s and A. Hajnal from 1985, we construct graphs of chromatic number $\omega_1$ without uncountable $\omega$-connected subgraphs. Second, we build triangle free graphs of chromatic number $\omega_1$ without subgraphs isomorphic to $H_{\omega,\omega+2}$.


title: A Cascadic Multigrid Method for Eigenvalue Problem

abstract: A cascadic multigrid method is proposed for eigenvalue problems based on the multilevel correction scheme. With this new scheme, an eigenvalue problem on the finest space can be solved by smoothing steps on a series of multilevel finite element spaces and eigenvalue problem solving on the coarsest finite element space. Choosing the appropriate sequence of finite element spaces and the number of smoothing steps, the optimal convergence rate with the optimal computational work can be arrived. Some numerical experiments are presented to validate our theoretical analysis.


title: SAD effects on grantsmanship

abstract: Seasonal affective disorder (SAD) is a state of depression induced by a lack of sufficient sunlight that occurs at high latitudes during the fall and winter. One effect of SAD is that causes people to be more risk-adverse, an effect that should be considered by granting agencies of high latitude countries. Funding agencies often have programmes aimed at high-risk, innovative research. However, the time of the year during which these purposefully high-risk proposals are evaluated usually does not take into consideration the effects of SAD. In high-latitude countries (e.g., Canada, UK, Nordic and Baltic countries), evaluating proposals for high-risk programmes during the late fall might significantly detract from the very purpose of such programmes. At this time of the year, grant evaluators might be in a darkness-induced state of mild depression. As such, evaluators might be more likely to opt for safe investments, more of the same, the well established, which is the antithesis of innovative research.


title: Game Theoretical Methods in Nonlinear PDEs

abstract: Nonlinear PDEs, mean value properties, and stochastic differential games are intrinsically connected. In this short expository note, we will describe how the solutions to certain PDEs (of $p$-Laplacian type) can be interpreted as limits of values of a specific Tug-of-War game, when the step-size $\epsilon$ determining the allowed length of move of a token, decreases to $0$.


title: A new time-frequency method to reveal quantum dynamics of atomic hydrogen in intense laser pulses: Synchrosqueezing Transform

abstract: This study introduces a new adaptive time-frequency (TF) analysis technique, synchrosqueezing transform (SST), to explore the dynamics of a laser-driven hydrogen atom at an {\it ab initio} level, upon which we have demonstrated its versatility as a new viable venue for further exploring quantum dynamics. For a signal composed of oscillatory components which can be characterized by instantaneous frequency, the SST enables rendering the decomposed signal based on the phase information inherited in the linear TF representation with mathematical support. Compared with the classical type TF methods, the SST clearly depicts several intrinsic quantum dynamical processes such as selection rules, AC Stark effects, and high harmonic generation.


title: Modeling effective FRW cosmologies with perfect fluids from states of the hybrid quantum Gowdy model

abstract: We employ recently developed approximation methods in the hybrid quantization of the Gowdy $T^3$ model with linear polarization and a massless scalar field to obtain physically interesting solutions of this inhomogeneous cosmology. More specifically, we propose approximate solutions of the quantum Gowdy model constructed in such a way that, for the Hamiltonian constraint, they effectively behave as those corresponding to a flat homogeneous and isotropic universe filled with a perfect fluid, even though these quantum states are far from being homogeneous and isotropic. We analyze how one can get different perfect fluid effective behaviors, including the cases of dust, radiation, and cosmological constant.


title: Path algebra algorithm for finding longest increasing subsequence

abstract: New algorithm for finding longest increasing subsequence is discussed. This algorithm is based on the ideas of idempotent mathematics and uses Max-Plus idempotent semiring. Problem of finding longest increasing sub- sequence is reformulated in a matrix form and solved with linear algebra.


title: T-shaped Single-photon Router

abstract: We study the transport properties of a single photon scattered by a two-level system (TLS) in a T-shaped waveguide, which is made of two coupled-resonator waveguides (CRWs)--- an infinite CRW and a semi-infinite CRW. The spontaneous emission of the TLS directs single photons from one CRW to the other. Although the transfer rate is different for the wave incident from different CRWs, due to the boundary breaking the translational symmetry, the boundary can enhance the transfer rate found in Phys. Rev. Lett. 111, 103604 (2013) and Phys. Rev. A 89, 013805 (2014), as the transfer rate could be unity for the wave incident from the semi-infinite CRW.


title: Sleeping Giants? - X-ray search for low-luminosity AGN candidates in nearby optically bright galaxies

abstract: In this Master's project, the X-ray nuclear properties of a sample of bright nearby galaxies are explored. This is done by matching their comprehensive optical spectroscopic classification to the latest available XMM-Newton catalogue - 3XMM-DR4. The good coverage (approx. 38 per cent) ensures that a statistically representative sample is investigated. All nuclear and morphological subsets found within the original sample of 486 galaxies are encompassed, but early-type galaxies and galaxies with optical features characteristic for active galactic nuclei (AGN) are favoured. The results from the investigation of the properties of our cross-matched sample are overall consistent with the presence of a large fraction of X-ray - detected low-luminosity AGN (LLAGN). The X-ray - detected galaxies within our HII and transition-LINER subsets are of particular interest, as they could harbour LLAGN missed by optical spectroscopic selection. The properties of these nuclei are explored by X-ray spectral fitting of available XMM-Newton observations. The majority of the examined spectra are found to be consistent with a thermal plasma model and an underlying power-law continuum, and the relative contributions of a potential LLAGN and of thermal emission to the source luminosity are presented. In the case that these spectra are not highly contaminated, as confirmed for some of the nuclei by examining available higher-resolution Chandra images, we classify approx. 43 per cent of the HII nuclei and 40 per cent of the transition-LINER nuclei as optically undetected LLAGN-candidates. This classification is found to be reasonably successful, given the project limitations, when compared to X-ray and multi-wavelength results found in literature for many of the nuclei.


title: Elliptic K3 surfaces associated with the product of two elliptic curves: Mordell-Weil lattices and their fields of definition

abstract: To a pair of elliptic curves, one can naturally attach two K3 surfaces: the Kummer surface of their product and a double cover of it, called the Inose surface. They have prominently featured in many interesting constructions in algebraic geometry and number theory. There are several more associated elliptic K3 surfaces, obtained through base change of the Inose surface; these have been previously studied by Kuwata. We give an explicit description of the geometric Mordell-Weil groups of each of these elliptic surfaces in the generic case (when the elliptic curves are non-isogenous). In the non-generic case, we describe a method to calculate explicitly a finite index subgroup of the Mordell-Weil group, which may be saturated to give the full group. Our methods rely on several interesting group actions, the use of rational elliptic surfaces, as well as connections to the geometry of low degree curves on cubic and quartic surfaces. We apply our techniques to compute the full Mordell-Weil group in several examples of arithmetic interest, arising from isogenous elliptic curves with complex multiplication, for which these K3 surfaces are singular.


title: Viscoelastic modulus reconstruction using time harmonic vibrations

abstract: This paper presents a new iterative reconstruction method to provide high-resolution images of shear modulus and viscosity via the internal measurement of displacement fields in tissues. To solve the inverse problem, we compute the Fr\'echet derivatives of the least-squares discrepancy functional with respect to the shear modulus and shear viscosity. The proposed iterative reconstruction method using this Fr\'echet derivative does not require any differentiation of the displacement data for the full isotropic linearly viscoelastic model, whereas the standard reconstruction methods require at least double differentiation. Because the minimization problem is ill-posed and highly nonlinear, this adjoint-based optimization method needs a very well-matched initial guess. We find a good initial guess. For a well-matched initial guess, numerical experiments show that the proposed method considerably improves the quality of the reconstructed viscoelastic images.


title: Real-time Information, Uncertainty and Quantum Feedback Control

abstract: Feedback is the core concept in cybernetics and its effective use has made great success in but not limited to the fields of engineering, biology, and computer science. When feedback is used to quantum systems, two major types of feedback control protocols including coherent feedback control (CFC) and measurement-based feedback control (MFC) have been developed. In this paper, we compare the two types of quantum feedback control protocols by focusing on the real-time information used in the feedback loop and the capability in dealing with parameter uncertainty. An equivalent relationship is established between quantum CFC and non-selective quantum MFC in the form of operator-sum representation. Using several examples of quantum feedback control, we show that quantum MFC can theoretically achieve better performance than quantum CFC in stabilizing a quantum state and dealing with Hamiltonian parameter uncertainty. The results enrich understanding of the relative advantages between quantum MFC and quantum CFC, and can provide useful information in choosing suitable feedback protocols for quantum systems.


title: Fuzzy bags, Polyakov loop and gauge/string duality

abstract: Confinement in SU($N$) gauge theory is due to the linear potential between colored objects. At short distances, the linear contribution could be considered as the quadratic correction to the leading Coulomb term. Recent lattice data show that such quadratic corrections also appear in the deconfined phase, in both the thermal quantities and the Polyakov loop. These contributions are studied systematically employing the gauge/string duality. "Confinement" in ${\mathcal N}=4$ SU($N$) Super Yang-Mills (SYM) theory could be achieved kinematically when the theory is defined on a compact space manifold. In the large-$N$ limit, deconfinement of ${\mathcal N}=4$ SYM on $\mathbb{S}^3$ at strong coupling is dual to the Hawking-Page phase transition in the global Anti-de Sitter spacetime. Meantime, all the thermal quantities and the Polyakov loop achieve significant quadratic contributions. Similar results can also be obtained at weak coupling. However, when confinement is induced dynamically through the local dilaton field in the gravity-dilaton system, these contributions can not be generated consistently. This is in accordance with the fact that there is no dimension-2 gauge-invariant operator in the boundary gauge theory. Based on these results, we suspect that quadratic corrections, and also confinement, should be due to global or non-local effects in the bulk spacetime.


title: Nonlinear optical magnetometry with accessible in situ optical squeezing

abstract: We demonstrate compact and accessible squeezed-light magnetometry using four-wave mixing in a single hot rubidium vapor cell. The strong intrinsic coherence of the four wave mixing process results in nonlinear magneto-optical rotation (NMOR) on each mode of a two mode relative-intensity squeezed state. This framework enables 4.7 dB of quantum noise reduction while the opposing polarization rotation signals of the probe and conjugate fields add to increase the total signal to noise ratio.


title: Adaptive Robust Optimization with Dynamic Uncertainty Sets for Multi-Period Economic Dispatch under Significant Wind

abstract: The exceptional benefits of wind power as an environmentally responsible renewable energy resource have led to an increasing penetration of wind energy in today's power systems. This trend has started to reshape the paradigms of power system operations, as dealing with uncertainty caused by the highly intermittent and uncertain wind power becomes a significant issue. Motivated by this, we present a new framework using adaptive robust optimization for the economic dispatch of power systems with high level of wind penetration. In particular, we propose an adaptive robust optimization model for multi-period economic dispatch, and introduce the concept of dynamic uncertainty sets and methods to construct such sets to model temporal and spatial correlations of uncertainty. We also develop a simulation platform which combines the proposed robust economic dispatch model with statistical prediction tools in a rolling horizon framework. We have conducted extensive computational experiments on this platform using real wind data. The results are promising and demonstrate the benefits of our approach in terms of cost and reliability over existing robust optimization models as well as recent look-ahead dispatch models.


title: Multiplicative properties of the number of $k$-regular partitions

abstract: In a previous paper of the second author with K. Ono, surprising multiplicative properties of the partition function were presented. Here, we deal with $k$-regular partitions. Extending the generating function for $k$-regular partitions multiplicatively to a function on $k$-regular partitions, we show that it takes its maximum at an explicitly described small set of partitions, and can thus easily be computed. The basis for this is an extension of a classical result of Lehmer, from which an inequality for the generating function for $k$-regular partitions is deduced which seems not to have been noticed before.


title: Gauge Invariance and QCD Twist-3 Factorization for Single Spin Asymmetries

abstract: The collinear factorization at twist-3 for Drell-Yan processes is studied with the motivation to solve the discrepancy in literature about the single spin asymmetry in the lepton angular distribution, and to show how QCD gauge invariance is realized in the hadronic tensor. The obtained result here agrees with our early result derived with a totally different approach. In addition to the asymmetry we can construct another two observables to identify the spin effect. We show that the gauge invariance of different contributions in the hadronic tensor is made in different ways by summing the effects of gluon exchanges. More interestingly is that we can show that the virtual correction to one structure function of the hadronic tensor, hence to some weighted SSA observables, is completely determined by the quark form factor. This will simplify the calculation of higher order corrections. The corresponding result in semi-inclusive DIS is also given for the comparison with Drell-Yan processes.


title: On Dark Matter Selected High-Scale Supersymmetry

abstract: The prediction for the Higgs mass in the dark matter selected high-scale SUSY is explored. We show the bounds on SUSY-breaking scale in models of SM $+\tilde{w}$ and SM $+\tilde{h}/\tilde{s}$ due to the observed Higgs mass at the LHC. We propose that effective theory below scale $\tilde{m}$ described by SM $+\tilde{w}$ is possibly realized in gauge mediation with multiple spurion fields that exhibit significant mass hierarchy, and that by SM $+\tilde{h}/\tilde{s}$ can be realized with direct singlet-messenger-messenger coupling for singlet Yukawa coupling $\lambda\sim(v/\tilde{m})^{1/2}g_{\text{SM}}$. Finally, the constraint on high-scale SUSY is investigated in the light of inflation physics if these two subjects are directly related.


title: Measurement-Based Noiseless Linear Amplification for Quantum Communication

abstract: Entanglement distillation is an indispensable ingredient in extended quantum communication networks. Distillation protocols are necessarily non-deterministic and require advanced experimental techniques such as noiseless amplification. Recently it was shown that the benefits of noiseless amplification could be extracted by performing a post-selective filtering of the measurement record to improve the performance of quantum key distribution. We apply this protocol to entanglement degraded by transmission loss of up to the equivalent of 100km of optical fibre. We measure an effective entangled resource stronger than that achievable by even a maximally entangled resource passively transmitted through the same channel. We also provide a proof-of-principle demonstration of secret key extraction from an otherwise insecure regime. The measurement-based noiseless linear amplifier offers two advantages over its physical counterpart: ease of implementation and near optimal probability of success. It should provide an effective and versatile tool for a broad class of entanglement-based quantum communication protocols.


title: Generation of electron beams from a laser-based advanced accelerator at Shanghai Jiao Tong University

abstract: At Shanghai Jiao Tong University, we have established a research laboratory for advanced acceleration research based on high-power lasers and plasma technologies. In a primary experiment based on the laser wakefield acceleration (LWFA) scheme, multi-hundred MeV electron beams having a reasonable quality are generated using 20-40 TW, 30 femtosecond laser pulses interacting independently with helium, neon, nitrogen and argon gas jet targets. The laser-plasma interaction conditions are optimized for stabilizing the electron beam generation from each type of gas. The electron beam pointing angle stability and divergence angle as well as the energy spectra from each gas jet are measured and compared.


title: A role for recurrent processing in object completion: neurophysiological, psychophysical and computational"evidence

abstract: Recognition of objects from partial information presents a significant challenge for theories of vision because it requires spatial integration and extrapolation from prior knowledge. We combined neurophysiological recordings in human cortex with psychophysical measurements and computational modeling to investigate the mechanisms involved in object completion. We recorded intracranial field potentials from 1,699 electrodes in 18 epilepsy patients to measure the timing and selectivity of responses along human visual cortex to whole and partial objects. Responses along the ventral visual stream remained selective despite showing only 9-25% of the object. However, these visually selective signals emerged ~100 ms later for partial versus whole objects. The processing delays were particularly pronounced in higher visual areas within the ventral stream, suggesting the involvement of additional recurrent processing. In separate psychophysics experiments, disrupting this recurrent computation with a backward mask at ~75ms significantly impaired recognition of partial, but not whole, objects. Additionally, computational modeling shows that the performance of a purely bottom-up architecture is impaired by heavy occlusion and that this effect can be partially rescued via the incorporation of top-down connections. These results provide spatiotemporal constraints on theories of object recognition that involve recurrent processing to recognize objects from partial information.


title: First Look Murchison Widefield Array observations of Abell 3667

abstract: The Murchison Widefield Array (MWA) is a new low frequency interferometric radio telescope, operating in the remote Murchison Radio Observatory in Western Australia. In this paper we present the first MWA observations of the well known radio relics in Abell 3667 (A3667) between 120 and 226 MHz. We clearly detect the radio relics in A3667 and present flux estimates and spectral indices for these features. The average spectral index of the north-west (NW) and south-east (SE) relics is -0.9 +/- 0.1 between 120 and 1400 MHz. We are able to resolve spatial variation in the spectral index of the NW relic from -1.7 to -0.4, which is consistent with results found at higher frequencies.


title: Collaborative Deep Learning for Recommender Systems

abstract: Collaborative filtering (CF) is a successful approach commonly used by many recommender systems. Conventional CF-based methods use the ratings given to items by users as the sole source of information for learning to make recommendation. However, the ratings are often very sparse in many applications, causing CF-based methods to degrade significantly in their recommendation performance. To address this sparsity problem, auxiliary information such as item content information may be utilized. Collaborative topic regression (CTR) is an appealing recent method taking this approach which tightly couples the two components that learn from two different sources of information. Nevertheless, the latent representation learned by CTR may not be very effective when the auxiliary information is very sparse. To address this problem, we generalize recent advances in deep learning from i.i.d. input to non-i.i.d. (CF-based) input and propose in this paper a hierarchical Bayesian model called collaborative deep learning (CDL), which jointly performs deep representation learning for the content information and collaborative filtering for the ratings (feedback) matrix. Extensive experiments on three real-world datasets from different domains show that CDL can significantly advance the state of the art.


title: Formation of Fermi surfaces and the appearance of liquid phases in holographic theories with hyperscaling violation

abstract: We consider a holographic fermionic system in which the fermions are interacting with a U(1) gauge field in the presence of a dilaton field in a gravity bulk of a charged black hole with hyperscaling violation. Using both analytical and numerical methods, we investigate the properties of the infrared and ultaviolet Green's functions of the holographic fermionic system. Studying the spectral functions of the system, we find that as the hyperscaling violation exponent is varied, the fermionic system possesses Fermi, non-Fermi, marginal-Fermi and log-oscillating liquid phases. Various liquid phases of the fermionic system with hyperscaling violation are also generated with the variation of the fermionic mass. We also explore the properties of the flat band and the Fermi surface of the non-relativistic fermionic fixed point dual to the hyperscaling violation gravity.


title: Near 3:2 and 2:1 mean motion resonances formation in the systems observed by Kepler

abstract: The Kepler mission has released ~4229 transiting planet candidates. There are approximately 222 candidate systems with three planets. Among them, the period ratios of planet pairs near 1.5 and 2.0 reveal that two peaks exist for which the proportions of the candidate systems are ~7.0% and 18.0%, respectively. In this work, we study the formation of mean motion resonance (MMR) systems, particularly for the planetary configurations near 3:2 and 2:1 MMRs, and we concentrate on the interplay between the resonant configuration and the combination of stellar accretion rate, stellar magnetic field, speed of migration and additional planets. We perform more than 1000 runs by assuming a system with a solar-like star and three surrounding planets. From the statistical results, we find that under the formation scenario, the proportions near 1.5 and 2.0 can reach 14.5% and 26.0%, respectively. In addition, $\dot M=0.1\times 10^{-8} ~M_\odot ~{\rm yr^{-1}}$ is propitious toward the formation of 3:2 resonance, whereas $\dot M=2\times 10^{-8} ~M_\odot ~{\rm yr^{-1}}$ contributes to the formation of 2:1 resonance. The speed-reduction factor of type I migration $f_1\geq 0.3$ facilitates 3:2 MMRs, whereas $f_1\geq 0.1$ facilitates 2:1 MMRs. If additional planets are present in orbits within the innermost or beyond the outermost planet in a three-planet system, 3:2:1 MMRs can be formed, but the original systems trapped in 4:2:1 MMRs are not affected by the supposed planets. In summary, we conclude that this formation scenario will provide a likely explanation for Kepler candidates involved in 2:1 and 3:2 MMRs.


title: Sunspot waves and flare energy release

abstract: We address a possibility of the flare process initiation and further maintenance of its energy release due to a transformation of sunspot longitudinal waves into transverse magnetic loop oscillations with initiation of reconnection. This leads to heating maintaining after the energy release peak and formation of a flat stage on the X-ray profile. We applied the time-distance plots and pixel wavelet filtration (PWF) methods to obtain spatio-temporal distribution of wave power variations in SDO/AIA data. To find magnetic waveguides, we used magnetic field extrapolation of SDO/HMI magnetograms. The propagation velocity of wave fronts was measured from their spatial locations at specific times. In correlation curves of the 17 GHz (NoRH) radio emission we found a monotonous energy amplification of 3-min waves in the sunspot umbra before the 2012 June 7 flare. This dynamics agrees with an increase in the wave-train length in coronal loops (SDO/AIA, 171 {\AA}) reaching the maximum 30 minutes prior to the flare onset. A peculiarity of this flare time profile in soft X-rays (RHESSI, 3-25 keV) is maintaining the constant level of the flare emission for 10 minutes after the short impulse phase, which indicates at the energy release continuation. Throughout this time, we found 30-sec period transverse oscillations of the flare loop in the radio-frequency range (NoRH, 17 GHz). This periodicity is apparently related to the transformation of propagating longitudinal 3-min waves from the sunspot into the loop transverse oscillations. The magnetic field extrapolation showed the existence of the magnetic waveguide (loop) connecting the sunspot with the energy release region. A flare loop heating can be caused by the interaction (reconnections) of this transversally oscillating waveguide with the underlying twisted loops.


title: Right $n$-angulated categories arising from covariantly finite subcategories

abstract: We define the notion of right $n$-angulated category, which generalizes the notion of right triangulated category. Let $\mathcal{C}$ be an additive category or $n$-angulated category and $\mathcal{X}$ a covariantly finite subcategory, we show that under certain conditions the quotient $\mathcal{C}/\mathcal{X}$ is a right $n$-angulated category. This result generalizes some previous work.


title: Curvaton with Nonminimal Derivative Coupling to Gravity II: Full Perturbation Analysis

abstract: In our previous work \cite{Feng:2013pba}, we have shown a curvaton model where the curvaton has a nonminimal derivative coupling to gravity. Such a coupling could bring us scale-invariance of the perturbations for wide range constant values of the equation-of-state of the cosmic background at the early time. In this paper, we continue our study by fully analyzing its perturbations up to the third order. Apart from the usual 2-point correlation function that has already been calculated in \cite{Feng:2013pba}, we have also taken into account the 3-point correlation functions including pure scalar part, pure tensor part, as well as the cross-correlations between scalar and tensor perturbation modes. We find that for pure scalar part, the 3-point correlation functions can generate non-Gaussianities that fits the PLANCK data very well. For pure tensor and mixed parts, the shape functions have peaks at squeezed and equilateral limits respectively, responsible for sizable $f_{NL}^{sqz}$ and $f_{NL}^{eql}$, which could be tested by the future observatioanl data.


title: Bounds on the origin of extragalactic ultrahigh energy cosmic rays from the IceCube neutrino observations

abstract: We study general implications of the IceCube observations in the energy range from $10^{6}$ GeV to $10^{10}$ GeV for the origin of extragalactic ultrahigh energy cosmic rays assuming that high energy neutrinos are generated by the photomeson production of protons in the extragalactic universe. The PeV-energy neutrino flux observed by IceCube gives strong bounds on the photomeson-production optical depth of protons in their sources and the intensity of the proton component of extragalactic cosmic rays. The neutrino flux implies that extragalactic cosmic-ray sources should have the optical depth greater than $\sim 0.01$ and contribute to more than a few percent of the observed bulk of cosmic rays at 10 PeV. If the spectrum of cosmic rays from these extragalactic sources extends well beyond 1 EeV, the neutrino flux indicates that extragalactic cosmic rays are dominant in the observed total cosmic-ray flux at 1 EeV and above, favoring the dip transition model of cosmic rays. The cosmic-ray sources are also required to be efficient neutrino emitters with the optical depth close to unity in this case. The highest energy cosmic-ray ($\sim 10^{11}$ GeV) sources should not be strongly evolved with redshift to account for the IceCube observations, suggesting that any cosmic-ray radiation scenarios involving distant powerful astronomical objects with strong cosmological evolution are strongly disfavored. These considerations conclude that none of the known extragalactic astronomical objects can be simultaneously a source of both PeV and trans-EeV energy cosmic rays. We also discuss a possible effect of cosmic-ray propagation in magnetized intergalactic space to the connection between the observed total cosmic-ray flux and neutrino flux.


title: Ages of Type Ia Supernovae Over Cosmic Time

abstract: We derive empirical models for galaxy mass assembly histories, and convolve these with theoretical delay time distribution (DTD) models for Type Ia supernovae (SNe Ia) to derive the distribution of progenitor ages for all SNe Ia occurring at a given epoch of cosmic time. In actively star-forming galaxies, the progression of the star formation rate is shallower than a $t^{-1}$ SN Ia DTD, so mean SN Ia ages peak at the DTD peak in all star-forming galaxies at all epochs of cosmic history. In passive galaxies which have ceased star formation through some quenching process, the SN Ia age distribution peaks at the quenching epoch, which in passive galaxies evolves in redshift to track the past epoch of major star formation. Our models reproduce the SN Ia rate evolution in redshift, the relationship between SN Ia stretch and host mass, and the distribution of SN Ia host masses in a manner qualitatively consistent with observations. Our model naturally predicts that low-mass galaxies tend to be actively star-forming while massive galaxies are generally passive, consistent with observations of galaxy "downsizing". Consequently, the mean ages of SNe Ia undergo a sharp transition from young ages at low host mass to old ages at high host mass, qualitatively similar to the transition of mean SN Ia Hubble residuals with host mass. The age discrepancy evolves with redshift in a manner currently not accounted for in SN Ia cosmology analyses. We thus suggest that SNe Ia selected only from actively star-forming galaxies will yield the most cosmologically uniform sample, due to the homogeneity of young SN Ia progenitor ages at all cosmological epochs.


title: Kinematic matrix theory and universalities in self-propellers and active swimmers

abstract: We describe an efficient and parsimonious matrix-based theory for studying the ensemble behavior of self-propellers and active swimmers, such as nanomotors or motile bacteria, that are typically studied by differential-equation-based Langevin or Fokker-Planck formalisms. The kinematic effects for elementary processes of motion are incorporated into a matrix, called the "kinematrix", from which we immediately obtain correlators and the mean and variance of angular and position variables (and thus effective diffusivity) by simple matrix algebra. The kinematrix formalism enables us recast the behaviors of a diverse range of self-propellers into a unified form, revealing universalities in their ensemble behavior in terms of new emergent time scales. Active fluctuations and hydrodynamic interactions can be expressed as an additive composition of separate self-propellers.


title: Global smooth axisymmetric solutions of 3-D Inhomogenenous incompressible Navier-Stokes system

abstract: In this paper, we investigate the global regularity to 3-D inhomogeneous incompressible Navier-Stokes system with axisymmetric initial data which does not have swirl component for the initial velocity. We first prove that the $L^\infty$ norm to the quotient of the inhomogeneity by $r,$ namely $a/r\eqdefa\bigl(1/\r-1\bigr)\bigl/r,$ controls the regularity of the solutions. Then we prove the global regularity of such solutions provided that the $L^\infty$ norm of $a_0/r$ is sufficiently small. Finally, with additional assumption that the initial velocity belongs to $L^p$ for some $p\in [1,2),$ we prove that the velocity field decays to zero with exactly the same rate as the classical Navier-Stokes system.


title: Sun-to-Earth Characteristics of Two Coronal Mass Ejections Interacting near 1 AU: Formation of a Complex Ejecta and Generation of a Two-Step Geomagnetic Storm

abstract: On 2012 September 30 - October 1 the Earth underwent a two-step geomagnetic storm. We examine the Sun-to-Earth characteristics of the coronal mass ejections (CMEs) responsible for the geomagnetic storm with combined heliospheric imaging and in situ observations. The first CME, which occurred on 2012 September 25, is a slow event and shows an acceleration followed by a nearly invariant speed in the whole Sun-Earth space. The second event, launched from the Sun on 2012 September 27, exhibits a quick acceleration, then a rapid deceleration and finally a nearly constant speed, a typical Sun-to-Earth propagation profile for fast CMEs \citep{liu13}. These two CMEs interacted near 1 AU as predicted by the heliospheric imaging observations and formed a complex ejecta observed at Wind, with a shock inside that enhanced the pre-existing southward magnetic field. Reconstruction of the complex ejecta with the in situ data indicates an overall left-handed flux rope-like configuration, with an embedded concave-outward shock front, a maximum magnetic field strength deviating from the flux rope axis and convex-outward field lines ahead of the shock. While the reconstruction results are consistent with the picture of CME-CME interactions, a magnetic cloud-like structure without clear signs of CME interactions \citep{lugaz14} is anticipated when the merging process is finished.


title: Gaussian Memory in Kinematic Matrix Theory for Self-Propellers

abstract: We extend the kinematic matrix ("kinematrix") formalism [Phys. Rev. E 89, 062304 (2014)], which via simple matrix algebra accesses ensemble properties of self-propellers influenced by uncorrelated noise, to treat Gaussian correlated noises. This extension brings into reach many real-world biological and biomimetic self-propellers for which inertia is significant. Applying the formalism, we analyze in detail ensemble behaviors of a 2D self-propeller with velocity fluctuations and orientation evolution driven by an Ornstein-Uhlenbeck process. On the basis of exact results, a variety of dynamical regimes determined by the inertial, speed-fluctuation, orientational diffusion, and emergent disorientation time scales are delineated and discussed.


title: Experimental Verification of Overlimiting Current by Surface Conduction and Electro-osmotic Flow in Microchannels

abstract: Possible mechanisms of overlimiting current in unsupported electrolytes, exceeding diffusion limitation, have been intensely studied for their fundamental significance and applications to desalination, separations, sensing, and energy storage. In bulk membrane systems, the primary physical mechanism is electro-convection, driven by electro-osmotic instability on the membrane surface. It has recently been predicted that confinement by charged surfaces in microchannels or porous media favors two new mechanisms, electro-osmotic flow (EOF) and surface conduction (SC), driven by large electric fields in the depleted region acting on the electric double layers on the sidewalls. Here, we provide the first direct evidence for the transition from SC to EOF above a critical channel height, using in situ particle tracking and current-voltage measurements in a micro/nanofluidic device. The dependence of the over-limiting conductance on channel depth (d) is consistent with theoretical predictions, scaling as d^-1 for SC and d^4/5 for EOF with a transition around d=8um. This complete picture of surface-driven over-limiting current can guide engineering applications of ion concentration polarization phenomena in microfluidics and porous media.


title: Topology and inference for multi-type Yule trees

abstract: We introduce two models for multi-type random trees motivated by studies of trait dependence in the evolution of species. Our discrete time model, the multi-type ERM tree, is a generalization of Markov propagation models on a random tree generated by a binary search or `equal rates Markov' mechanism. Our continuous time model, the multi-type Yule tree with mutations, is a multi-type generalization of the tree generated by a pure birth or Yule process. We study type dependent topological properties of these two random tree models. We derive asymptotic results that allow one to infer model parameters from data on types at the leaves and at branch-points that are one step away from the leaves.


title: Terahertz radiation by optical rectification in a hydrogen-bonded organic molecular ferroelectric crystal, 2-phenylmalondialdehyde

abstract: Terahertz radiation by optical rectification has been observed at room temperature in a hydrogen-bonded organic molecular ferroelectric crystal, 2-phenyl malondialdehyde (PhMDA). The radiated electromagnetic wave consisted of a single-cycle terahertz pulse with a temporal width of $\sim$ 0.5 ps. The terahertz radiation amplitude divided by the sample thickness in PhMDA was nearly equivalent to that in a typical terahertz wave emitter ZnTe. This is attributable to a long coherence length in the range of 130 $\sim$ 800 $\mu$m for the terahertz radiation from PhMDA. We also discussed the possibility of PhMDA as a terahertz wave emitter in terms of the phase-matching condition.


title: Anthropic considerations in nuclear physics

abstract: In this short review, I discuss the sensitivity of the generation of the light and the life-relevant elements like carbon and oxygen under changes of the parameters of the Standard Model pertinent to nuclear physics. Chiral effective field theory allows for a systematic and precise description of the forces between two, three, and four nucleons. In this framework, variations under the light quark masses and the electromagnetic fine structure constant can also be consistently calculated. Combining chiral nuclear effective field theory with Monte Carlo simulations allows to further calculate the properties of nuclei, in particular of the Hoyle state in carbon, that plays a crucial role in the generation of the life-relevant elements in hot, old stars. The dependence of the triple-alpha process on the fundamental constants of Nature is calculated and some implications for our anthropic view of the Universe are discussed.


title: Visualization of ferroelectric domains in boracite using emission of terahertz radiation

abstract: We report on the emission of terahertz radiation by irradiation of femtosecond laser pulses in non-centrosymmetric paraelectric and ferroelectric phases of Co$_3$B$_7$O$_{13}$I boracite. The Generation of the terahertz waves in both phases is caused by optical rectification via a second-order nonlinear optical effect. In the ferroelectric phase, we successfully visualized ferroelectric domains by analyzing the polarization state of the terahertz wave radiated from the crystal. In a large area of the crystal ($\sim$ 500 $\times$ 500 $\mu$m$^2$), the observed polarization vector of the radiated terahertz wave was tilted from directions of spontaneous polarization, i.e., [100]$_{\rm cub}$, [010]$_{\rm cub}$, and [001]$_{\rm cub}$ in cubic setting, which can be explained by the presence of a ferroelectric 90$^\circ$ domain wall of the (101)$_{\rm cub}$ plane.


title: Formation of ultracold $^{7}$Li$^{85}$Rb molecules in the lowest triplet electronic state by photoassociation and their detection by ionization spectroscopy

abstract: We report the formation of ultracold $^{7}$Li$^{85}$Rb molecules in the $a^{3}\Sigma^{+}$ electronic state by photoassociation (PA) and their detection via resonantly enhanced multiphoton ionization (REMPI). With our dual-species Li and Rb magneto-optical trap (MOT) apparatus, we detect PA resonances with binding energies up to ~62 cm$^{-1}$ below the $^{7}$Li 2s $^{2}S_{1/2}$ $+$ $^{85}$Rb 5p $^{2}P_{1/2}$ asymptote. In addition, we use REMPI spectroscopy to probe the $a^{3}\Sigma^{+}$ state and excited electronic $3^{3} \Pi$ and $4^{3} \Sigma^{+}$ states, and identify $a^{3}\Sigma^{+} (v" = 7 - 13)$, $3^{3} \Pi (v'_{\Pi} = 0 - 10)$ and $4^{3} \Sigma^{+} (v'_{\Sigma} = 0 - 5)$ vibrational levels. Our line assignments agree well with ab initio calculations. These preliminary spectroscopic studies on previously unobserved electronic states are crucial to discovering transition pathways for transferring ultracold LiRb molecules created via PA to deeply bound rovibrational levels of the electronic ground state.


title: Nonlinear generalized sections of vector bundles

abstract: We present an extension of J. F. Colombeau's theory of nonlinear generalized functions to spaces of generalized sections of vector bundles. Our construction builds on classical functional analytic notions, which is the key to having a canonical geometric embedding of vector bundle valued distributions into spaces of generalized sections. This permits to have tensor products, invariance under diffeomorphisms, covariant derivatives and the sheaf property. While retaining as much compatibility to L. Schwartz' theory of distributions as possible, our theory provides the basis for a rigorous and general treatment of singular pseudo-Riemannian geometry in the setting of Colombeau nonlinear generalized functions.


title: Spectral energy distributions of the BL Lac PKS 2155-304 from XMM-Newton

abstract: We have used all 20 archival XMM-Newton observations of PKS 2155-304 with simultaneous X-ray and UV/optical data to study its long term flux and spectral variability. We find significant variations, in all bands, on time-scales of years with an rms amplitude of ~ 35-45 per cent, through the optical/UV variations are not correlated with those in the X-ray. We constructed SEDs that span more than three orders of magnitude in frequency and we first fitted them with a log-parabolic model; such models have been applied many times in the past for this, and other, blazars. These fits were poor, so we then examined combined power-law and log-parabolic fits that are improvements. These models indicate that the optical/UV and X-ray flux variations are mainly driven by model normalization variations, but the X-ray band flux is also affected by spectral variations, as parametrized with the model 'curvature' parameter, b. Overall, the energy at which the emitted power is maximum correlates positively with the total flux. As the spectrum shifts to higher frequencies, the spectral 'curvature' increase, in contrast to what is expected if a single log-parabolic model were an acceptable representation of the broad band SEDs. Our results suggest that the optical/UV and X-ray emissions in this source may arise from different lepton populations.