title: Task-Based Optimization of Computed Tomography Imaging Systems

abstract: The goal of this thesis is to provide a framework for the use of task-based metrics of image quality to aid in the design, implementation, and evaluation of CT image reconstruction algorithms and CT systems in general. We support the view that task-based metrics of image quality can be useful in guiding the algorithm design and implementation process in order to yield images of objectively superior quality and higher utility for a given task. Further, we believe that metrics such as the Hotelling observer (HO) SNR can be used as summary scalar metrics of image quality for the evaluation of images produced by novel reconstruction algorithms. In this work, we aim to construct a concise and versatile formalism for image reconstruction algorithm design, implementation, and assessment. The bulk of the work focuses on linear analytical algorithms, specifically the ubiquitous filtered back-projection (FBP) algorithm. However, due to the demonstrated importance of optimization-based algorithms in a wide variety of CT applications, we devote one chapter to the characterization of noise properties in TV-based iterative reconstruction, as the understanding of image statistics in optimization-based reconstruction is the limiting factor in applying HO metrics.


title: On Diffusion Limited Deposition

abstract: We propose a simple model of columnar growth through {\it diffusion limited aggregation} (DLA). Consider a graph $G_N\times\N$, where the basis has $N$ vertices $G_N:=\{1,\dots,N\}$, and two vertices $(x,h)$ and $(x',h')$ are adjacent if $|h-h'|\le 1$. Consider there a simple random walk {\it coming from infinity} which {\it deposits} on a growing cluster as follows: the cluster is a collection of columns, and the height of the column first hit by the walk immediately grows by one unit. Thus, columns do not grow laterally. We prove that there is a critical time scale $N/\log(N)$ for the maximal height of the piles, i.e., there exist constants $\alpha<\beta$ such that the maximal pile height at time $\alpha N/\log(N)$ is of order $\log(N)$, while at time $\beta N/\log(N)$ is larger than $N^\chi$. This suggests that a \emph{monopolistic regime} starts at such a time and only the highest pile goes on growing. If we rather consider a walk whose height-component goes down deterministically, the resulting \emph{ballistic deposition} has maximal height of order $\log(N)$ at time $N$. These two deposition models, diffusive and ballistic, are also compared with uniform random allocation and Polya's urn.


title: Higgs boson production in association with a jet at NNLO using jettiness subtraction

abstract: We use the recently proposed jettiness-subtraction scheme to provide the complete calculation of Higgs boson production in association with a jet in hadronic collisions through next-to-next-to-leading order in perturbative QCD. This method exploits the observation that the $N$-jettiness event-shape variable completely describes the singularity structure of QCD when final-state colored particles are present. Our results are in agreement with a recent computation of the $gg$ and $qg$ partonic initial states based on sector-improved residue subtraction. We present phenomenological results for both fiducial cross sections and distributions at the LHC.


title: The bi-Hamiltonian cohomology of a scalar Poisson pencil

abstract: We compute the bi-Hamiltonian cohomology of an arbitrary dispersionless Poisson pencil in a single dependent variable using a spectral sequence method. As in the KdV case, we obtain that $BH^p_d(\hat{F}, d_1,d_2)$ is isomorphic to $\mathbb{R}$ for $(p,d)=(0,0)$, to $C^\infty (\mathbb{R})$ for $(p,d)=(1,1)$, $(2,1)$, $(2,3)$, $(3,3)$, and vanishes otherwise.


title: Hologram of a pure state black hole

abstract: In this paper we extend the HKLL holographic smearing function method to reconstruct (quasi)local AdS bulk scalar observables in the background of a large AdS black hole formed by null shell collapse (a "pure state" black hole), from the dual CFT which is undergoing a sudden quench. In particular, we probe the near horizon and sub-horizon bulk locality. First we construct local bulk operators from the CFT in the leading semiclassical limit, $N\rightarrow\infty$. Then we look at effects due to the finiteness of $N$, where we propose a suitable coarse-graining prescription involving early and late time cut-offs to define semiclassical bulk observables which are approximately local; their departure from locality being non-perturbatively small in $N$. Our results have important implications on the black hole information problem.


title: A Morita theorem for modular finite W-algebras

abstract: We consider the Lie algebra $\mathfrak{g}$ of a simple, simply connected algebraic group over a field of large positive characteristic. For each nilpotent orbit $\mathcal{O} \subseteq \mathfrak{g}$ we choose a representative $e\in \mathcal{O}$ and attach a certain filtered, associative algebra $\widehat{U}(\mathfrak{g},e)$ known as a finite $W$-algebra, defined to be the opposite endomorphism ring of the generalised Gelfand-Graev module associated to $(\mathfrak{g}, e)$. This is shown to be Morita equivalent to a certain central reduction of the enveloping algebra of $U(\mathfrak{g})$. The result may be seen as a modular version of Skryabin's equivalence.


title: Closed-Form Bounds for the Rice $Ie$-Function

abstract: This work is devoted in the derivation of novel upper and lower bounds for the Rice $Ie$-function. These bounds are expressed in closed-form and are shown to be quite tight. This is particularly evident by the fact that for a certain range of parameter values, the derived lower bound virtually behaves as a remarkably accurate approximation. As a result, the offered expressions can be considered useful mathematical tools that can be efficiently employed in various analytical studies related to natural sciences and engineering. To this effect, they can be sufficiently applied in the area of digital communications over fading channels for the derivation of explicit representations for vital performance measures such as bit and symbol error probability, among others.


title: Pinball Loss Minimization for One-bit Compressive Sensing: Convex Models and Algorithms

abstract: The one-bit quantization is implemented by one single comparator that operates at low power and a high rate. Hence one-bit compressive sensing (1bit-CS) becomes attractive in signal processing. When measurements are corrupted by noise during signal acquisition and transmission, 1bit-CS is usually modeled as minimizing a loss function with a sparsity constraint. The one-sided $\ell_1$ loss and the linear loss are two popular loss functions for 1bit-CS. To improve the decoding performance on noisy data, we consider the pinball loss, which provides a bridge between the one-sided $\ell_1$ loss and the linear loss. Using the pinball loss, two convex models, an elastic-net pinball model and its modification with the $\ell_1$-norm constraint, are proposed. To efficiently solve them, the corresponding dual coordinate ascent algorithms are designed and their convergence is proved. The numerical experiments confirm the effectiveness of the proposed algorithms and the performance of the pinball loss minimization for 1bit-CS.


title: An Approach to Data Prefetching Using 2-Dimensional Selection Criteria

abstract: We propose an approach to data memory prefetching which augments the standard prefetch buffer with selection criteria based on performance and usage pattern of a given instruction. This approach is built on top of a pattern matching based prefetcher, specifically one which can choose between a stream, a stride, or a stream followed by a stride. We track the most recently called instructions to make a decision on the quantity of data to prefetch next. The decision is based on the frequency with which these instructions are called and the hit/miss rate of the prefetcher. In our approach, we separate the amount of data to prefetch into three categories: a high degree, a standard degree and a low degree. We ran tests on different values for the high prefetch degree, standard prefetch degree and low prefetch degree to determine that the most optimal combination was 1, 4, 8 lines respectively. The 2 dimensional selection criteria improved the performance of the prefetcher by up to 9.5% over the first data prefetching championship winner. Unfortunately performance also fell by as much as 14%, but remained similar on average across all of the benchmarks we tested.


title: Misrepair mechanism: a mechanism essential for individual adaptation, species adaptation and species evolution

abstract: In Misrepair-accumulation theory, we have proposed a Misrepair mechanism for interpreting aging. However, Misrepair mechanism is also important in biological adaptation. Misrepair is a strategy of repair for increasing the surviving chance of an organism when it suffers from severe injuries. As a surviving strategy, Misrepair mechanism plays also an important role in individual adaptation, species survival, and species differentiation. Firstly, Misrepair of an injury is one of the manners of individual adaptation; and Misrepair mechanism gives an organism a great potential in adapting to changeable and destructive environment. Secondly, Misrepair mechanism is important in maintaining and enlarging the diversity of genome DNAs of a species; and a large diversity of genome DNAs is essential for the adaptation and the differentiation of a species in different environments. On one hand, somatic Misrepairs are essential for maintaining the sufficient number of individuals in a species, which are the vectors of different genome DNAs. On the other hand, Misrepair of DNA is a source of DNA mutations and the DNA Misrepairs in germ cells may contribute to the diversity of genome DNAs in a species. In conclusion, Misrepair mechanism is a mechanism essential for individual adaptation, species adaptation, and species evolution.


title: Cosmological Particle Production at Strong Coupling

abstract: We study the dynamics of a strongly-coupled quantum field theory in a cosmological spacetime using the holographic AdS/CFT correspondence. Specifically we consider a confining gauge theory in an expanding FRW universe and track the evolution of the stress-energy tensor during a period of expansion, varying the initial temperature as well as the rate and amplitude of the expansion. At strong coupling, particle production is inseparable from entropy production. As a result, we find significant qualitative differences from the weak coupling results: at strong coupling the system rapidly loses memory of its initial state as the amplitude is increased. Furthermore, in the regime where the Hubble parameter is parametrically smaller than the initial temperature, the dynamics is well modelled as a plasma evolving hydrodynamically towards equilibrium.


title: Potential of longevity: hidden in structural complexity

abstract: In order to understand the phenomenon of longevity in biological world, the relationship between the potential of longevity and the structural complexity of an organism is analyzed. I. The potential of longevity is the maximum lifespan of an organism if it lives in an ideal environment. The longevity of an organism includes two parts: the time for development (development time) and the time for structure-maintenance (maintenance time). II. The development time of an organism depends on its structural complexity. The maintenance time is related to two factors: the degree of damage-exposure and the potential of functionality for structure-maintenance. The potential of functionality is built in structural complexity. Since both of development time and maintenance time are determined by structural complexity, the potential of longevity of an organism is hidden in structural complexity. III. The individuals of different species have different potentials of longevity because they have different structural complexities. An animal has limited longevity because it has limited structural complexity. Limited structural complexity and limited longevity are essential for the survival of a species. IV. Despite having the same potential of longevity, the individuals of a species can have different lifespans. The lifespan of an individual is more related to the degree of damage-exposure, which is determined by the living environment and the living habit of the individual. In conclusion, the potential of longevity of an organism is hidden in structural complexity, but the real lifespan of an organism is more related to the living environment.


title: Full quantum state reconstruction of symmetric two-mode squeezed thermal states via spectral homodyne detection

abstract: We suggest and demonstrate a scheme to reconstruct the symmetric two-mode squeezed thermal states of spectral sideband modes from an optical parametric oscillator. The method is based on a single homodyne detector and active stabilization of the cavity. The measurement scheme have been successfully tested on different two-mode squeezed thermal states, ranging from uncorrelated coherent states to entangled states.


title: Fibers of Cyclic Covering Fibrations of a Ruled Surface

abstract: We give an algorithm to classify singular fibers of finite cyclic covering fibrations of a ruled surface by using singularity diagrams. As the first application, we classify all fibers of 3-cyclic covering fibrations of genus 4 of a ruled surface and show that the signature of a complex surface with this fibration is non-positive by computing the local signature for any fiber. As the second application, we classify all fibers of hyperelliptic fibrations of genus 3 into 12 types according to the Horikawa index. We also prove that finite cyclic covering fibrations of a ruled surface have no multiple fibers if the degree of the covering is greater than 3.


title: Premature aging as a consequence of Mis-construction of tissues and organs during body development

abstract: Hutchinson-Gilford Progeria syndrome, Werner syndrome, and Cockayne syndrome are three genetic disorders, in which the children have premature aging features. To understand the phenomenon of premature aging, the similarity of aging features in these syndromes to that in normal aging is investigated. Although these three syndromes have different genetic backgrounds, all the patients have abnormal structures of tissues/organs like that in normal aging. Therefore, the abnormality in tissue structure is the common point in premature aging and normal aging. This abnormality links also a defective development and a defective repair, the Misrepair. Defective development is a result of Mis-construction of the structure of tissues and organs as consequence of genetic mutations. Aging is a result of Mis-reconstructions, the Misrepairs, for maintaining the structure of tissues/organs. Construction-reconstruction of the structure of an organism is thus the coupling point of development and aging. Mis- construction and Mis-reconstruction (Misrepair) are the essential processes for the development of aging-like feathers. In conclusion, premature aging is a result of Mis- construction of tissues and organs during body development as consequence of genetic disorders.


title: Training generative neural networks via Maximum Mean Discrepancy optimization

abstract: We consider training a deep neural network to generate samples from an unknown distribution given i.i.d. data. We frame learning as an optimization minimizing a two-sample test statistic---informally speaking, a good generator network produces samples that cause a two-sample test to fail to reject the null hypothesis. As our two-sample test statistic, we use an unbiased estimate of the maximum mean discrepancy, which is the centerpiece of the nonparametric kernel two-sample test proposed by Gretton et al. (2012). We compare to the adversarial nets framework introduced by Goodfellow et al. (2014), in which learning is a two-player game between a generator network and an adversarial discriminator network, both trained to outwit the other. From this perspective, the MMD statistic plays the role of the discriminator. In addition to empirical comparisons, we prove bounds on the generalization error incurred by optimizing the empirical MMD.


title: Velocity autocorrelation in liquid para-hydrogen by quantum simulations for first-principle computations of the neutron cross sections

abstract: Accurate knowledge of the single-molecule (self) translational dynamics of liquid para-H2 is an essential requirement for the calculation of the neutron scattering properties of this important quantum liquid. We show that, by using Centroid Molecular Dynamics (CMD) quantum simulations of the velocity autocorrelation function, calculations of the total neutron cross section (TCS) remarkably agree with experimental data at the thermal and epithermal incident neutron energies where para-H2 dynamics is actually dominated by the self contributions. This result shows that a proper account of the quantum nature of the fluid, as provided by CMD, is a necessary and very effective condition to obtain the correct absolute-scale cross section values directly from first-principle computations of the double differential cross section, and without the need of introducing any empirically adjusted quantity. At subthermal incident energies, appropriate modeling of the para-H2 intermolecular (distinct) dynamics also becomes crucial, but quantum simulations are not yet able to cope with it. Existing simple models which account for the distinct part provide an appropriate correction of self-only calculations and bring the computed results in reasonable accord with TCS experimental data available until very recently. However, if just published cross section measurements in the cold range are considered, the agreement turns out to be by far superior and very satisfactory. The possible origin of slight residual differences will be commented and suggest further computational and experimental efforts. Nonetheless, the ability to reproduce the total cross section in the wide range between 1 and 900 meV represents an encouraging and important validation step of the CMD method and of the present simple algorithm.


title: Efficient adaptive MCMC through precision estimation

abstract: A novel adaptive Markov chain Monte Carlo algorithm is presented. The algorithm utilizes sparsity in the partial correlation structure of a density to efficiently estimate the covariance matrix through the Cholesky factor of the precision matrix. The algorithm also utilizes the sparsity to sample efficiently from both MALA and Metropolis Hasting random walk proposals. Further, an algorithm that estimates the partial correlation structure of a density is proposed. Combining this with the Cholesky factor estimation algorithm results in an efficient black-box AMCMC method that can be used for general densities with unknown dependency structure. The method is compared with regular empirical covariance adaption for two examples. In both examples, the proposed method's covariance estimates converge faster to the true covariance matrix and the computational cost for each iteration is lower.


title: O$_2$ Emission Toward Orion H$_2$ Peak 1 and the Role of FUV-Illuminated C-Shocks

abstract: Molecular oxygen, O_2, has been the target of ground-based and space-borne searches for decades. Of the thousands of lines of sight surveyed, only those toward Rho Oph and Orion H_2 Peak 1 have yielded detections of any statistical significance. The detection of the O_2 N_J =3_3 -1_2 and 5_4 - 3_4 lines at 487.249 GHz and 773.840 GHz, respectively, toward Rho Ophiuchus has been attributed to a short-lived peak in the time-dependent, cold-cloud O_2 abundance, while the detection of the O_2 N_J =3_3 - 1_2, 5_4 - 3_4 lines, plus the 7_6 - 5_6 line at 1120.715 GHz, toward Orion has been ascribed to time-dependent preshock physical and chemical evolution and low-velocity (12 km/s) non-dissociative C-type shocks, both of which are fully shielded from far-ultraviolet (FUV) radiation, plus a postshock region that is exposed to a FUV field. We report a re-interpretation of the Orion O_2 detection based on new C-type shock models that fully incorporate the significant effects the presence of even a weak FUV field can have on the preshock gas, shock structure and postshock chemistry. In particular, we show that a family of solutions exists, depending on the FUV intensity, that reproduces both the observed O_2 intensities and O_2 line ratios. The solution in closest agreement with the shock parameters inferred for H_2 Peak 1 from other gas tracers assumes a 23 km/s shock impacting gas with a preshock density of 8x10^4 cm^-3 and G_0 =1, substantially different from that inferred for the fully-shielded shock case. As pointed out previously, the similarity between the LSR velocity of all three O_2 lines (~11 km/s) and recently measured H_2O 5_32 - 4_41 maser emission at 620.701 GHz toward H_2 Peak 1 suggests that the O_2 emission arises behind the same shocks responsible for the maser emission, though the O_2 emission is almost certainly more extended than the localized high density maser spots


title: Excitonic effects in two-dimensional semiconductors: Path integral Monte Carlo approach

abstract: One of the most striking features of novel 2D semiconductors (e.g., transition metal dichalcogenide monolayers or phosphorene) is a strong Coulomb interaction between charge carriers resulting in large excitonic effects. In particular, this leads to the formation of multi-carrier bound states upon photoexcitation (e.g., excitons, trions and biexcitons), which could remain stable at near-room temperatures and contribute significantly to optical properties of such materials. In the present work we have used the Path Integral Monte Carlo methodology to numerically study properties of multi-carrier bound states in 2D semiconductors. Specifically, we have accurately investigated and tabulated the dependence of single exciton, trion and biexciton binding energies on the strength of dielectric screening, including the limiting cases of very strong and very weak screening. The results of this work are potentially useful in the analysis of experimental data and benchmarking of theoretical and computational models.


title: Discretization effects in $N_c=2$ QCD and Random Matrix Theory

abstract: We summarize the analytical solution of the Chiral Perturbation Theory for the Hermitian Wilson Dirac operator of $N_c=2$ QCD with quarks in the fundamental representation. Results have been obtained for the quenched microscopic spectral density, the distribution of the chiralities over the real modes and the chiral condensate. The analytical results are compared with results from a Monte Carlo simulation of the corresponding Random Matrix Theory.


title: Interacting bosons in two-dimensional flat band systems

abstract: The Hubbard model of bosons on two dimensional lattices with a lowest flat band is discussed. In these systems there is a critical density, where the ground state is known exactly and can be represented as a charge density wave. Above this critical filling, depending on the lattice structure and the interaction strength, the additional particles are either delocalised and condensate in the ground state, or they form pairs. Pairs occur at strong interactions, e.g., on the chequerboard lattice. The general mechanism behind this phenomenon is discussed.


title: Finite Element Methods for the Stochastic Allen-Cahn Equation with Gradient-type Multiplicative Noises

abstract: This paper studies finite element approximations of the stochastic Allen-Cahn equation with gradient-type multiplicative noises that are white in time and correlated in space. The sharp interface limit as the parameter $\epsilon \rightarrow 0$ of the stochastic equation formally approximates a stochastic mean curvature flow which is described by a stochastically perturbed geometric law of the deterministic mean curvature flow. Both the stochastic Allen-Cahn equation and the stochastic mean curvature flow arise from materials science, fluid mechanics and cell biology applications. Two fully discrete finite element methods which are based on different time-stepping strategies for the nonlinear term are proposed. Strong convergence with sharp rates for both fully discrete finite element methods is proved. This is done with a crucial help of the H\"{o}lder continuity in time with respect to the spatial $L^2$-norm and $H^1$-seminorm for the strong solution of the stochastic Allen-Cahn equation, which are key technical lemmas proved in paper. It also relies on the fact that high moments of the strong solution are bounded in various spatial and temporal norms. Numerical experiments are provided to gauge the performance of the proposed fully discrete finite element methods and to study the interplay of the geometric evolution and gradient-type noises.


title: Estimation and Simulation of a COGARCH(p,q) model in the YUIMA project

abstract: In this paper we show how to simulate and estimate a COGARCH(p,q) model in the R package yuima. Several routines for simulation and estimation are available. Indeed for the generation of a COGARCH(p,q) trajectory, the user can choose between two alternative schemes. The first is based on the Euler discretization of the stochastic differential equations that identifies a COGARCH(p,q) model while the second one considers the explicit solution of the variance process. Estimation is based on the matching of the empirical with the theoretical autocorrelation function. In this case three different approaches are implemented: minimization of the mean square error, minimization of the absolute mean error and the generalized method of moments where the weighting matrix is continuously updated. Numerical examples are given in order to explain methods and classes used in the yuima package.


title: Gravitational self-force in non-vacuum spacetimes: an effective field theory derivation

abstract: In this paper we investigate the motion of small compact objects in non-vacuum spacetimes using methods from effective field theory in curved spacetime. Although a vacuum formulation is sufficient in many astrophysical contexts, there are applications such as the role of the self-force in enforcing cosmic-censorship in the context of the overcharging problem, which necessitate an extension into the non-vacuum regime. The defining feature of the self-force problem in non-vacuum spacetimes is the coupling between gravitational and non-gravitational field perturbations. The formulation of the self-force problem for non-vacuum spacetimes was recently provided in simultaneous papers by Zimmerman and Poisson [1] and Linz, Friedmann, Wiseman [2]. Here we distinguish ourselves by working with the effective action rather than the field equations. The formalism utilizes the multi-index notation developed by Zimmerman and Poisson [1] to accommodate the coupling between the different fields. Using dimensional regularization, we arrive at a finite expression for the local self-force expressed in terms of multi-index quantities evaluated in the background spacetime. We then apply the formalism to compute the coupled gravitational self-force in two explicit cases. First, we calculate the self-force on a massive particle possessing scalar charge and moving in an scalarvac spacetime. We then derive an expression for the self-force on an electrically charged, massive particle moving in an electrovac spacetime. In both cases, the force is expressed as a sum of local terms involving tensors defined in the background spacetime and evaluated at the current position of the particle, as well as tail integrals that depend on the past history of the particle.


title: Correcting systematic polarization effects in Keck LRISp spectropolarimetry to <0.05%

abstract: Spectropolarimetric measurements at moderate spectral resolutions are effective tracers of stellar magnetic fields and circumstellar environments when signal to noise ratios (SNRs) above 2000 can be achieved. The LRISp spectropolarimeter is capable of achieving these SNRs on faint targets with the 10m aperture of the Keck telescope, provided several instrumental artifacts can be suppressed. We describe here several methods to overcome instrumental error sources that are required to achieve these high SNRs on LRISp. We explore high SNR techniques such as defocusing and slit-stepping during integration with high spectral and spatial oversampling. We find that the instrument flexure and interference fringes introduced by the achromatic retarders create artificial signals at 0.5\% levels in the red channel which mimic real stellar signals and limit the sensitivity and calibration stability of LRISp. Careful spectral extraction and data filtering algorithms can remove these error sources. For faint targets and long exposures, cosmic ray hits are frequent and present a major limitation to the upgraded deep depletion red-channel CCD. These must be corrected to the same high SNR levels, requiring careful spectral extraction using iterative filtering algorithms. We demonstrate here characterization of these sources of instrumental polarization artifacts and present several methods used to successfully overcome these limitations. We have measured the linear to circular cross-talk and find it to be roughly 5\%, consistent with the known instrument limitations. We show spectropolarimetric signals on brown dwarfs are clearly detectable at 0.2\% amplitudes with sensitivities better than 0.05\% at full spectral sampling in atomic and molecular bands. Future LRISp users can perform high sensitivity observations with high quality calibration when following the described algorithms.


title: General Riemannian SOM

abstract: Kohonen's Self-Organizing Maps (SOMs) have proven to be a successful data-reduction method to identify the intrinsic lower-dimensional sub-manifold of a data set that is scattered in the higher-dimensional feature space. Motivated by the possibly non-Euclidian nature of the feature space and of the intrinsic geometry of the data set, we extend the definition of classic SOMs to obtain the General Riemannian SOM (GRiSOM). We additionally provide an implementation as a proof-of-concept for geometries with constant curvature. We furthermore perform the analytic and numerical analysis of the stability limits of certain (GRi)SOM configurations covering the different possible regular tessellation of the map space in each geometry. A deviation between the numerical and analytic stability limit has been observed for the square and hexagonal Euclidean maps for very small neighbourhoods in the map space as well as agreement in case of longer-ranged relations between the map nodes.


title: Quantum Process Tomography of an Optically-Controlled Kerr Non-linearity

abstract: Any optical quantum information processing machine would be comprised of fully-characterized constituent devices for both single state manipulations and tasks involving the interaction between multiple quantum optical states. Ideally for the latter, would be an apparatus capable of deterministic optical phase shifts that operate on input quantum states with the action mediated solely by auxiliary signal fields. Here we present the complete experimental characterization of a system designed for optically controlled phase shifts acting on single-photon level probe coherent states. Our setup is based on a warm vapor of rubidium atoms under the conditions of electromagnetically induced transparency with its dispersion properties modified through the use of an optically triggered N-type Kerr non-linearity. We fully characterize the performance of our device by sending in a set of input probe states and measuring the corresponding output via balanced homodyne tomography and subsequently performing the technique of coherent state quantum process tomography. This method provides us with the precise knowledge of how our optical phase shift will modify any arbitrary input quantum state engineered in the mode of the reconstruction.


title: Some applications of weighted norm inequalities to the error analysis of PDE constrained optimization problems

abstract: The purpose of this work is to illustrate how the theory of Muckenhoupt weights, Muckenhoupt weighted Sobolev spaces and the corresponding weighted norm inequalities can be used in the analysis and discretization of PDE constrained optimization problems. We consider: a linear quadratic constrained optimization problem where the state solves a nonuniformly elliptic equation; a problem where the cost involves pointwise observations of the state and one where the state has singular sources, e.g. point masses. For all three examples we propose and analyze numerical schemes and provide error estimates in two and three dimensions. While some of these problems might have been considered before in the literature, our approach allows for a simpler, Hilbert space-based, analysis and discretization and further generalizations.


title: Equivariant Poincar\'e series and topology of valuations

abstract: The equivariant with respect to a finite group action Poincar\'e series of a collection of $r$ valuations was defined earlier as a power series in $r$ variables with the coefficients from a modification of the Burnside ring of the group. Here we show that (modulo simple exceptions) the equivariant Poincar\'e series determines the equivariant topology of the collection of valuations.


title: Helium-3 and Helium-4 acceleration by high power laser pulses for hadron therapy

abstract: The laser driven acceleration of ions is considered a promising candidate for an ion source for hadron therapy of oncological diseases. Though proton and carbon ion sources are conventionally used for therapy, other light ions can also be utilized. Whereas carbon ions require 400 MeV per nucleon to reach the same penetration depth as 250 MeV protons, helium ions require only 250 MeV per nucleon, which is the lowest energy per nucleon among the light ions. This fact along with the larger biological damage to cancer cells achieved by helium ions, than that by protons, makes this species an interesting candidate for the laser driven ion source. Two mechanisms (Magnetic Vortex Acceleration and hole-boring Radiation Pressure Acceleration) of PW-class laser driven ion acceleration from liquid and gaseous helium targets are studied with the goal of producing 250 MeV per nucleon helium ion beams that meet the hadron therapy requirements. We show that He3 ions, having almost the same penetration depth as He4 with the same energy per nucleon, require less laser power to be accelerated to the required energy for the hadron therapy.


title: Sensitivity of full-sky experiments to large scale cosmic ray anisotropies

abstract: The two main advantages of space-based observation of extreme energy ($\gtrsim5\times10^{19}$ eV) cosmic rays (EECRs) over ground based observatories are the increased field of view and the full-sky coverage with nearly uniform systematics across the entire sky. The former guarantees increased statistics, whereas the latter enables a clean partitioning of the sky into spherical harmonics. The discovery of anisotropies would help to identify the long sought origin of EECRs. We begin an investigation of the reach of a full-sky space-based experiment such as EUSO to detect anisotropies in the extreme-energy cosmic-ray sky compared to ground based partial-sky experiments such as the Pierre Auger Observatory and Telescope Array. The technique is explained here, and simulations for a Universe with just two nonzero multipoles, monopole plus either dipole or quadrupole, are presented. These simulations quantify the advantages of space-based, all-sky coverage.


title: Spectral dimension and Bohr's formula for Schrodinger operators on unbounded fractal spaces

abstract: We establish an asymptotic formulas for the eigenvalue counting function of the Schr\"odinger operator $-\Delta +V$ for some unbounded potentials $V$ on several types of unbounded fractal spaces. We give sufficient conditions for Bohr's formula to hold on metric measure spaces which admit a cellular decomposition, and then verify these conditions for fractafolds and fractal fields based on nested fractals. In particular, we partially answer a question of Fan, Khandker, and Strichartz regarding the spectral asymptotics of the harmonic oscillator potential on the infinite blow-up of a Sierpinski gasket.

categories: math-ph math.FA math.MG math.MP math.PR math.SP


title: $k$-center Clustering under Perturbation Resilience

abstract: The $k$-center problem is a canonical and long-studied facility location and clustering problem with many applications in both its symmetric and asymmetric forms. Both versions of the problem have tight approximation factors on worst case instances. Therefore to improve on these ratios, one must go beyond the worst case. In this work, we take this approach and provide strong positive results both for the asymmetric and symmetric $k$-center problems under a natural input stability (promise) condition called $\alpha$-perturbation resilience [Bilu and Linia 2012], which states that the optimal solution does not change under any alpha-factor perturbation to the input distances. We provide algorithms that give strong guarantees simultaneously for stable and non-stable instances: our algorithms always inherit the worst-case guarantees of clustering approximation algorithms, and output the optimal solution if the input is $2$-perturbation resilient. Furthermore, we prove our result is tight by showing symmetric $k$-center under $(2-\epsilon)$-perturbation resilience is hard unless $NP=RP$. The impact of our results are multifaceted. This is the first tight result for any problem under perturbation resilience. Furthermore, our results illustrate a surprising relationship between symmetric and asymmetric $k$-center instances under perturbation resilience. Unlike approximation ratio, for which symmetric $k$-center is easily solved to a factor of 2 but asymmetric $k$-center cannot be approximated to any constant factor, both symmetric and asymmetric $k$-center can be solved optimally under resilience to 2-perturbations. Finally, our guarantees in the setting where only part of the data satisfies perturbation resilience makes these algorithms more applicable to real-life instances.


title: $B \rightarrow D l \nu$ Form Factors at Non-Zero Recoil and Extraction of $|V_{cb}|$

abstract: We present a lattice QCD calculation of the $B \rightarrow D l \nu$ semileptonic decay form factors $f_+(q^2)$ and $f_0(q^2)$ for the entire physical $q^2$ range. Non-relativistic QCD (NRQCD) bottom quarks and Highly Improved Staggered Quark (HISQ) charm and light quarks are employed together with $N_f = 2+1$ MILC gauge configurations. A joint fit to our lattice and BaBar experimental data allows an extraction of the CKM matrix element $|V_{cb}|$. We also determine the phenomenologically interesting ratio $R(D) = {\cal B}(B \rightarrow D \tau \nu_\tau) / {\cal B}(B \rightarrow D l \nu_l)$ ($l = e, \mu$). We find $|V_{cb}|_{excl.}^{B \rightarrow D} = 0.0402(17)(13)$, where the first error consists of the lattice simulation errors and the experimental statistical error and the second error is the experimental systematic error. For the branching fraction ratio we find $R(D) = 0.300(8)$.


title: Spin-current probe for phase transition in an insulator

abstract: Spin fluctuation and transition have always been one of central topics of magnetism and condensed matter science. Experimentally, the spin fluctuation is found transcribed onto scattering intensity in the neutron scattering process, which is represented by dynamical magnetic susceptibility and maximized at phase transitions. Importantly, a neutron carries spin without electric charge, and it can bring spin into a sample without being disturbed by electric energy, although large facilities such as a nuclear reactor is necessary. Here we show that spin pumping, frequently used in nanoscale spintronic devices, provides a desktop micro probe for spin transition; spin current is a flux of spin without an electric charge and its transport reflects spin excitation. We demonstrate detection of antiferromagnetic transition in ultra-thin CoO films via frequency dependent spin-current transmission measurements, which provides a versatile probe for phase transition in an electric manner in minute devices.


title: Momentum-space calculation of ${}^4$He triatomic system with realistic potential

abstract: The system of three ${}^4\mathrm{He}$ atoms with realistic interactions is studied in the momentum-space framework. It is demonstrated that short and long-range difficulties encountered in the coordinate-space calculations can be reliably resolved in the momentum-space calculations. Well-converged and accurate results are obtained for the ground and excited trimer energies, atom-dimer scattering length, phase shifts, inelasticity parameters, and elastic and breakup cross sections. A significant correction to previously published results is found for the elastic cross section at very low energy.


title: Hypoelliptic heat kernel on nilpotent Lie groups

abstract: The starting point of our analysis is an old idea of writing an eigenfunction expansion for a heat kernel considered in the case of a hypoelliptic heat kernel on a nilpotent Lie group $G$. One of the ingredients of this approach is the generalized Fourier transform. The formula one gets using this approach is explicit as long as we can find all unitary irreducible representations of $G$. In the current paper we consider an $n$-step nilpotent Lie group $G_{n}$ as an illustration of this technique. First we apply Kirillov's orbit method to describe these representations for $G_{n}$. This allows us to write the corresponding hypoelliptic heat kernel using an integral formula over a Euclidean space. As an application, we describe a short-time behavior of the hypoelliptic heat kernel in our case.


title: Surfaces of general type with q=2 are rigidified

abstract: Let $S$ be a minimal smooth projective surface of general type with irregularity $q=2$. We show that, if $S$ has a nontrivial holomorphic automorphism acting trivially on the cohomology with rational coefficients, then it is a surface isogenous to a product. As a consequence of this geometric characterization, one infers that no nontrivial automorphism of surfaces of general type with $q=2$ (which are not necessarily minimal) can be homotopic to the identity. In particular, such surfaces are rigidified in the sense of Fabrizio Catanese.


title: Determination of the relative phase in \psi' and J/\psi\ decays into baryon and antibaryon

abstract: With the recent measurements of psi' and J/psi decay into octet-baryon pairs, we study the relative phase between the strong and the electromagnetic amplitudes, and find a large phase by fitting the present data. The fits take into account the details of experimental effects, including energy spread and initial state radiation. We also predict some branching fractions of psi' decays and the continuum production rates at the J/psi mass based on the relative phase and absolute amplitudes obtained from the fits.


title: Robust Biomolecular Finite Automata

abstract: We present a uniform method for translating an arbitrary nondeterministic finite automaton (NFA) into a deterministic mass action input/output chemical reaction network (I/O CRN) that simulates it. The I/O CRN receives its input as a continuous time signal consisting of concentrations of chemical species that vary to represent the NFA's input string in a natural way. The I/O CRN exploits the inherent parallelism of chemical kinetics to simulate the NFA in real time with a number of chemical species that is linear in the size of the NFA. We prove that the simulation is correct and that it is robust with respect to perturbations of the input signal, the initial concentrations of species, the output (decision), and the rate constants of the reactions of the I/O CRN.


title: Using Ensemble Models in the Histological Examination of Tissue Abnormalities

abstract: Classification models for the automatic detection of abnormalities on histological samples do exists, with an active debate on the cost associated with false negative diagnosis (underdiagnosis) and false positive diagnosis (overdiagnosis). Current models tend to underdiagnose, failing to recognize a potentially fatal disease. The objective of this study is to investigate the possibility of automatically identifying abnormalities in tissue samples through the use of an ensemble model on data generated by histological examination and to minimize the number of false negative cases.


title: Directed Flow of Charged Kaons in Au+Au Collisions from the BES Program at RHIC

abstract: We report the measurement of the directed flow ($v_{1}$) for charged kaons in Au+Au collisions at $\sqrt{s_{NN}}$ =7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV as a function of rapidity and compare these results for pions, protons and antiprotons. These new kaon results may help to constrain the medium properties and collision dynamics including the in-medium kaon potential and baryon number transport in these collisions.


title: Textual Spatial Cosine Similarity

abstract: When dealing with document similarity many methods exist today, like cosine similarity. More complex methods are also available based on the semantic analysis of textual information, which are computationally expensive and rarely used in the real time feeding of content as in enterprise-wide search environments. To address these real-time constraints, we developed a new measure of document similarity called Textual Spatial Cosine Similarity, which is able to detect similitude at the semantic level using word placement information contained in the document. We will see in this paper that two degenerate cases exist for this model, which coincide with Cosine Similarity on one side and with a paraphrasing detection model to the other.


title: Bound states of charges on top of graphene in magnetic field

abstract: We show theoretically that in the external magnetic field like charges on top of graphene monolayer may be mutually attracted to form macro-molecules. For this to happen graphene needs to be in Quantum Hall plateau state with local chemical potential being between the Landau levels. Graphene electron(s) gets localized in the middle between charges and provides overscreening of Coulomb repulsion between the charges. The size of the resulting macro-molecules is of the order of the magnetic length ($\sim 10$ nm for magnetic field 10 T). The possible stable macro-molecules that unit charges can form on graphene in magnetic field are classified. The binding survives significant temperatures, exceeding mobility barriers for many ionically bond impurities. The influence of possible lattice-scale effects of valley-mixing are discussed. Tuning the doping of graphene or the magnetic field, the binding of impurities can be turned on and off and the macro-molecule size may be tuned. This opens the perspective to nanoscopic manipulation of ions on graphene by using magnetic field and gating.


title: Combined Fast Reversible Liquid-like Elastic Deformation with Topological Phase Transition in Na$_3$Bi

abstract: By means of first-principles calculations, we identified the structural phase transition of Na$_3$Bi from hexagonal ground state to cubic $cF$16 phase above 0.8 GPa, in agreement with the experimental findings. Upon the releasing of pressure, \emph{cF}16 phase of Na$_3$Bi is mechanically stable at ambient condition. The calculations revealed that the $cF$16 phase is topological semimetal, in similarity to well-known HgTe and it even exhibits an unusually low $C^\prime$ modulus (only about 1.9 GPa) and a huge anisotropy, $A^u$ of as high as 11, the third highest value among all known cubic crystals in their elastic behaviors. These facts render \emph{cF}16-type Na$_3$Bi very soft with a liquid-like elastic deformation in the (110)$<$1$\overline{1}$0$>$ slip system. Importantly, as accompanied with this deformation, Na$_3$Bi shows a topological phase transition from a topological semimetal state at its strain-free cubic phase to a topological insulating state at its distorted phase. Because the $C^\prime$ elastic deformation almost costs no energy in a reversible and liquid-like soft manner, \emph{cF}16-type Na$_3$Bi would potentially provide a fast on/off switching way between topological insulator and topological semimetal, which would be beneficial to the quantum electronic devices for practical applications.


title: Experimental demonstration of PT-symmetric stripe lasers

abstract: Recently, the coexistence of parity-time (PT) symmetric laser and absorber has gained tremendous research attention. While the PT symmetric absorber has been observed in microwave metamaterials, the experimental demonstration of PT symmetric laser is still absent. Here we experimentally study PT-symmetric laser absorber in stripe waveguide. Using the concept of PT symmetry to exploit the light amplification and absorption, PT-symmetric laser absorbers have been successfully obtained. Different from the single-mode PT symmetric lasers, the PT-symmetric stripe lasers have been experimentally confirmed by comparing the relative wavelength positions and mode spacing under different pumping conditions. When the waveguide is half pumped, the mode spacing is doubled and the lasing wavelengths shift to the center of every two initial lasing modes. All these observations are consistent with the theoretical predictions and confirm the PT-symmetry breaking well.


title: SPDEs with two reflecting walls and two singular drifts

abstract: We study SPDEs with two reflecting walls $\Lambda^1$, $\Lambda^2$ and two singular drifts $\frac{c_1}{(X-\Lambda^1)^{\vartheta}}$, $\frac{c_2}{(\Lambda^2-X)^{\vartheta}}$, driven by space-time white noise. First, we establish the existence and uniqueness of the solutions $X$ for $\vartheta\geq 0$. Second, we obtain the following pathwise properties of the solutions $X$. If $\vartheta>3$, then a.s. $\Lambda^1<X<\Lambda^2$ for all $t\geq0$; If $0<\vartheta<3$, then $X$ hits $\Lambda^1$ or $\Lambda^2$ with positive probability in finite time. Thus $\vartheta=3$ is the critical parameter for $X$ to hit reflecting walls.


title: Hulls of Surfaces

abstract: In this paper it is shown that every compact two-dimensional manifold $S$, with or without boundary, can be embedded in $\mathbb C^3$ as a smooth submanifold $\Sigma$ in such a way that the polynomially convex hull of $\Sigma$, though strictly larger than $\Sigma$, contains no analytic disc.


title: Exponential mixing for stochastic model of two-dimensional second grade fluids

abstract: In this paper, we establish the exponential mixing property of stochastic models for the incompressible second grade fluid. The general criterion established by Cyril Odasso plays an important role.


title: Mismatched Estimation in Large Linear Systems

abstract: We study the excess mean square error (EMSE) above the minimum mean square error (MMSE) in large linear systems where the posterior mean estimator (PME) is evaluated with a postulated prior that differs from the true prior of the input signal. We focus on large linear systems where the measurements are acquired via an independent and identically distributed random matrix, and are corrupted by additive white Gaussian noise (AWGN). The relationship between the EMSE in large linear systems and EMSE in scalar channels is derived, and closed form approximations are provided. Our analysis is based on the decoupling principle, which links scalar channels to large linear system analyses. Numerical examples demonstrate that our closed form approximations are accurate.


title: Structures in 9Be, 10Be and 10B studied with tensor-optimized shell model

abstract: We investigate the structures of $^{9,10}$Be and $^{10}$B with the tensor-optimized shell model (TOSM) using the effective interaction based on the bare nucleon-nucleon interaction AV8$^\prime$. The tensor correlation is treated in TOSM with the full optimization of 2p2h configurations including high-momentum components. The short-range correlation is described in the unitary correlation operator method (UCOM). It is found that the level orders of the low-lying states of $^{9,10}$Be and $^{10}$B are entirely reproduced. For $^9$Be, ground band states are located relatively in higher energy than the experiments, which indicates the missing $\alpha$ clustering correlation in these states as seen in the case of $^8$Be with TOSM. In addition, the tensor force gives the larger attraction for $T$=1/2 states than for $T$=3/2 ones for $^9$Be. For $^{10}$Be, the tensor contribution of $0^+_2$ shows the largest value among the $0^+$ states. This can be related to the $\alpha$ clustering correlation in this state. It is also found that the level order of three nuclei depends on the tensor force in comparison with the results obtained with the Minnesota interaction without the tensor force.


title: Towards Noncommutative Linking Numbers Via the Seiberg-Witten Map

abstract: In the present work some geometric and topological implications of noncommutative Wilson loops are explored via the Seiberg-Witten map. In the abelian Chern-Simons theory on a three dimensional manifold, it is shown that the effect of noncommutativity is the appearance of $6^n$ new knots at the $n$-th order of the Seiberg-Witten expansion. These knots are trivial homology cycles which are Poincar\'e dual to the high-order Seiberg-Witten potentials. Moreover the linking number of a standard 1-cycle with the Poincar\'e dual of the gauge field is shown to be written as an expansion of the linking number of this 1-cycle with the Poincar\'e dual of the Seiberg-Witten gauge fields. In the process we explicitly compute the noncommutative 'Jones-Witten' invariants up to first order in the noncommutative parameter. Finally in order to exhibit a physical example, we apply these ideas explicitly to the Aharonov-Bohm effect. It is explicitly displayed at first order in the noncommutative parameter, we also show the relation to the noncommutative Landau levels.


title: Quantum interference and structure-dependent orbital-filling effects on the thermoelectric properties of quantum dot molecules

abstract: The quantum interference and orbital filling effects on the thermoelectric (TE) properties of quantum dot molecules with high figure of merit are illustrated via the full solution to the Hubbard- Anderson model in the Coulomb blockade regime. It is found that under certain condition in the triangular QD molecule (TQDM), destructive quantum interference (QI) can occur, which leads to vanishing small electrical conductance, while the Seebeck coefficient is modified dramatically. When TQDM is in the charge localization state due to QI, the Seebeck coefficient is seriously suppressed at low temperature, but highly enhanced at high temperature. Meanwhile, the behavior of Lorenz number reveals that it is easier to block charge transport via destructive QI than the electron heat transport at high temperatures. The maximum power factor (PF) in TQDM occurs at full-filling condition. Nevertheless, low-filling condition is preferred for getting maximum PF in serially coupled triple QDs in general. In double QDs, the maximum PF can be achieved either with orbital-depletion or orbital-filling as a result of electron-hole symmetry. Our theoretical work provides a useful guideline for advancing the nanoscale TE technology.


title: High-energy magnetic excitations in overdoped La$_{2-x}$Sr$_{x}$CuO$_{4}$ studied by neutron and resonant inelastic X-ray scattering

abstract: We have performed neutron inelastic scattering and resonant inelastic X-ray scattering (RIXS) at the Cu-$L_3$ edge to study high-energy magnetic excitations at energy transfers of more than 100 meV for overdoped La$_{2-x}$Sr$_{x}$CuO$_{4}$ with $x=0.25$ ($T_c=15$ K) and $x=0.30$ (non-superconducting) using identical single crystal samples for the two techniques. From constant-energy slices of neutron scattering cross-sections, we have identified magnetic excitations up to ~250 meV for $x=0.25$. Although the width in the momentum direction is large, the peak positions along the (pi, pi) direction agree with the dispersion relation of the spin-wave in the non-doped La$_{2}$CuO$_{4}$ (LCO), which is consistent with the previous RIXS results of cuprate superconductors. Using RIXS at the Cu-$L_3$ edge, we have measured the dispersion relations of the so-called paramagnon mode along both (pi, pi) and (pi, 0) directions. Although in both directions the neutron and RIXS data connect with each other and the paramagnon along (pi, 0) agrees well with the LCO spin-wave dispersion, the paramagnon in the (pi, pi) direction probed by RIXS appears to be less dispersive and the excitation energy is lower than the spin-wave of LCO near (pi/2, pi/2). Thus, our results indicate consistency between neutron inelastic scattering and RIXS, and elucidate the entire magnetic excitation in the (pi, pi) direction by the complementary use of two probes. The polarization dependence of the RIXS profiles indicates that appreciable charge excitations exist in the same energy range of magnetic excitations, reflecting the itinerant character of the overdoped sample. A possible anisotropy in the charge excitation intensity might explain the apparent differences in the paramagnon dispersion in the (pi, pi) direction as detected by the X-ray scattering.


title: Block Markov Superposition Transmission of RUN Codes

abstract: In this paper, we propose a simple procedure to construct (decodable) good codes with any given alphabet (of moderate size) for any given (rational) code rate to achieve any given target error performance (of interest) over additive white Gaussian noise (AWGN) channels. We start with constructing codes over groups for any given code rates. This can be done in an extremely simple way if we ignore the error performance requirement for the time being. Actually, this can be satisfied by repetition (R) codes and uncoded (UN) transmission along with time-sharing technique. The resulting codes are simply referred to as RUN codes for convenience. The encoding/decoding algorithms for RUN codes are almost trivial. In addition, the performance can be easily analyzed. It is not difficult to imagine that a RUN code usually performs far away from the corresponding Shannon limit. Fortunately, the performance can be improved as required by spatially coupling the RUN codes via block Markov superposition transmission (BMST), resulting in the BMST-RUN codes. Simulation results show that the BMST-RUN codes perform well (within one dB away from Shannon limits) for a wide range of code rates and outperform the BMST with bit-interleaved coded modulation (BMST-BICM) scheme.


title: Nano-Stitching of Graphene Bilayers: A First-Principles Study

abstract: A nano-stitching method is proposed and investigated to modify graphene bilayers. Based this method, four types of low energy carbon allotropes, "wormhole graphene" allotropoes, are obtained and their structures, stabilities and electronic properties are investigated using first principles methods. We find that all of these wormhole graphene allotropoes are more favorable than graphdiyne and dynamically stable. Similar to carbon nanotubes and fullerences, these graphene allotropes are expected to act as two-dimensional periodic nano-capsules for encapsulating magnetic atoms or functional clusters for a variety of applications.


title: Electronic and magnetic influences of a stacking fault in cobalt nanoscale islands on the Ag(111) surface

abstract: Utilizing spin-polarized scanning tunneling microscopy and spectroscopy, we found coexistence of perpendicularly and in-plane magnetized cobalt nanoscale islands on the Ag(111) surface, and the relationship between the moire corrugation amplitude and the magnetization direction of the islands; the islands with the stronger moire corrugation show the perpendicular magnetization, and the ones with the weaker moire corrugation do the in-plane. Density functional theory calculations reproduce the relationship and explain the differences between the two types of the islands with an fcc stacking fault in the intrinsic hcp stacking of cobalt.


title: The Dependence of the $A_V$ Prior for SN\,Ia on Host Mass and Disk Inclination

abstract: Supernovae type Ia (SNIa) are used as "standard candles" for cosmological distance scales. To fit their light curve shape -- absolute luminosity relation, one needs to assume an intrinsic color and a likelihood of host galaxy extinction or a convolution of these, a color distribution prior. The host galaxy extinction prior is typically assumed to be an exponential drop-off for the current supernova programs ($P(A_V) \propto e^{-A_V/\tau_0}$). We explore the validity of this prior using the distribution of extinction values inferred when two galaxies accidentally overlap (an occulting galaxy pair). We correct the supernova luminosity distances from the SDSS-III Supernova projects (SDSS-SN) by matching the host galaxies to one of three templates from occulting galaxy pairs based on the host galaxy mass and the $A_V$-bias - prior-scale ($\tau_0$) relation from Jha et al. (2007). We find that introducing an $A_V$ prior that depends on host mass results in lowered luminosity distances for the SDSS-SN on average but it does not reduce the scatter in individual measurements. This points, in our view, to the need for many more occulting galaxy templates to match to SNIa host galaxies to rule out this possible source of scatter in the SNIa distance measurements. We match occulting galaxy templates based on both mass and projected radius and we find that one should match by stellar mass first with radius as a secondary consideration. We discuss the caveats of the current approach and our aim is to convince the reader that a library of occulting galaxy pairs observed with HST will provide sufficient priors to improve (optical) SNIa measurements to the next required accuracy in Cosmology.


title: Logics of Strong Noncontingency

abstract: Inspired by Hintikka's treatment of question embedding verbs in [8] and the variations of noncontingency operator, we propose a logic with strong noncontingency operator $\blacktriangle$ as the only primitive modality. A proposition is strongly noncontingent, if no matter whether it is true or false, it does it necessarily; otherwise, it is weakly contingent. This logic is not a normal modal logic, since $\blacktriangle(\phi\to\psi)\to(\blacktriangle\phi\to\blacktriangle\psi)$ is invalid. We compare the relative expressivity of this logic and other logics, such as standard modal logic, noncontingency logic, and logic of essence and accident, and investigate its frame definability. Apart from these results, we also propose a suitable notion of bisimulation for the logic of strong noncontingency, based on which we characterize this logic within modal logic and within first-order logic. We also axiomatize the logic of strong noncontingency over various frame classes. Our work is also related to the treatment of agreement operator in [10].


title: Logarithmic Corrections to the Entanglement Entropy

abstract: In a $d$-dimensional conformal field theory, it has been known that a relevant deformation operator with the conformal dimension, $\Delta=\frac{d+2}{2}$, generates a logarithmic correction to the entanglement entropy. In the large 't Hooft coupling limit, we can investigate such a logarithmic correction holographically by deforming an AdS space with a massive scalar field dual to the operator with $\Delta=\frac{d+2}{2}$. There are two sources generating the logarithmic correction. One is the metric deformation and the other is the minimal surface deformation. In this work, we investigate the change of the entanglement entropy caused by the minimal surface deformation and find that the second order minimal surface deformation leads to an additional logarithmic correction.


title: A non-commutative Beurling's theorem with respect to unitarily invariant norms

abstract: In 1967, Arveson invented a non-commutative generalization of classical $H^{\infty},$ known as finite maximal subdiagonal subalgebras, for a finite von Neumann algebra $\mathcal M$ with a faithful normal tracial state $\tau$. In 2008, Blecher and Labuschagne proved a version of Beurling's theorem on $H^\infty$-right invariant subspaces in a non-commutative $L^{p}(\mathcal M,\tau)$ space for $1\le p\le \infty$. In the present paper, we define and study a class of norms ${\mathcal{N}}_{c}(\mathcal M, \tau)$ on $\mathcal{M},$ called normalized, unitarily invariant, $\Vert \cdot \Vert_{1}$-dominating, continuous norms, which properly contains the class $\{ \Vert \cdot \Vert_{p}:1\leq p< \infty \}.$ For $\alpha \in \mathcal{N}_{c}(\mathcal M, \tau),$ we define a non-commutative $L^{\alpha }({\mathcal{M}},\tau)$ space and a non-commutative $H^{\alpha}$ space. Then we obtain a version of the Blecher-Labuschagne-Beurling invariant subspace theorem on $H^\infty$-right invariant subspaces in a non-commutative $L^{\alpha }({\mathcal{M}},\tau)$ space. Key ingredients in the proof of our main result include a characterization theorem of $H^\alpha$ and a density theorem for $L^\alpha(\mathcal M,\tau)$.


title: A Theory of Formal Synthesis via Inductive Learning

abstract: Formal synthesis is the process of generating a program satisfying a high-level formal specification. In recent times, effective formal synthesis methods have been proposed based on the use of inductive learning. We refer to this class of methods that learn programs from examples as formal inductive synthesis. In this paper, we present a theoretical framework for formal inductive synthesis. We discuss how formal inductive synthesis differs from traditional machine learning. We then describe oracle-guided inductive synthesis (OGIS), a framework that captures a family of synthesizers that operate by iteratively querying an oracle. An instance of OGIS that has had much practical impact is counterexample-guided inductive synthesis (CEGIS). We present a theoretical characterization of CEGIS for learning any program that computes a recursive language. In particular, we analyze the relative power of CEGIS variants where the types of counterexamples generated by the oracle varies. We also consider the impact of bounded versus unbounded memory available to the learning algorithm. In the special case where the universe of candidate programs is finite, we relate the speed of convergence to the notion of teaching dimension studied in machine learning theory. Altogether, the results of the paper take a first step towards a theoretical foundation for the emerging field of formal inductive synthesis.


title: Formation of primordial supermassive stars by burst accretion

abstract: A promising formation channel of SMBHs at redshift 6 is the so-called DC model, which posits that a massive seed BH forms through gravitational collapse of a $\sim 10^5~M_\odot$ SMS. We study the evolution of such a SMS growing by rapid mass accretion. In particular, we examine the impact of time-dependent mass accretion of repeating burst and quiescent phases that are expected to occur with a self-gravitating circumstellar disk. We show that the stellar evolution with such episodic accretion differs qualitatively from that expected with a constant accretion rate, even if the mean accretion rate is the same. Unlike the case of constant mass accretion, whereby the star expands roughly following $R_* \simeq 2.6 \times 10^3 R_\odot (M_*/100~M_\odot)^{1/2}$, the protostar can substantially contract during the quiescent phases between accretion bursts. The stellar effective temperature and ionizing photon emissivity increase accordingly as the star contracts, which can cause strong ionizing feedback and halt the mass accretion onto the star. With a fixed duration of the quiescent phase $\Delta t_{\rm q}$, such contraction occurs in early evolutionary phases, i.e. for $M_* \lesssim 10^3~M_\odot$ with $\Delta t_{\rm q} \simeq 10^3$ yr. For later epochs and larger masses but the same $\Delta t_{\rm q}$, contraction is negligible even during quiescent phases. With larger quiescent times $\Delta t_{\rm q}$, however, the star continues to contract during quiescent phases even for the higher stellar masses. We show that such behavior is well understood by comparing the interval time and the thermal relaxation time for a bloated surface layer. We conclude that the UV radiative feedback becomes effective if the quiescent phase associated by the burst accretion is longer than $\sim 10^3$ yr, which is possible in an accretion disk forming in the direct collapse model.


title: On the equation-of-motion versus in-in approach in cosmological perturbation theory

abstract: In this paper, we study several issues in the linear equation-of-motion (EoM) and in-in approaches of computing the two-point correlation functions in multi-field inflation. We prove the equivalence between this EoM approach and the first-principle in-in formalism. We check this equivalence using several explicit examples, including cases with scale-invariant corrections and scale-dependent features. Motivated by the explicit proof, we show that the usual procedures in these approaches can be extended and applied to some interesting model categories beyond what has been studied in the literature so far. These include the density perturbations with strong couplings and correlated multi-field initial states.


title: From 0D to 1D spatial models using OCMat

abstract: We show that the standard class of optimal control models in OCMat can be used to analyze 1D spatial distributed systems. This approach is an intermediate step on the way to the FEM discretization approach presented in Grass and Uecker (2015). Therefore, the spatial distributed model is transformed into a standard model by a finite difference discretization. This (high dimensional) standard model is then analyzed using OCMAT. As an example we apply this method to the spatial distributed shallow lake model formulated in Brock and Xepapadeas (2008). The results are then compared with those of the FEM discretization in GRass and Uecker (2015)


title: On some extensions of the Ailon-Rudnick Theorem

abstract: In this paper we present some extensions of the Ailon-Rudnick Theorem, which says that if $f,g\in{\mathbb C}[T]$, then $\gcd(f^n-1,g^m-1)$ is bounded for all $n,m\ge 1$. More precisely, using a uniform bound for the number of torsion points on curves and results on the intersection of curves with algebraic subgroups of codimension at least $2$, we present two such extensions in the univariate case. We also give two multivariate analogues of the Ailon-Rudnick Theorem based on Hilbert's irreducibility theorem and a result of Granville and Rudnick about torsion points on hypersurfaces.


title: Vacuum polarization around a three-dimensional black hole

abstract: We calculate the Euclidean propagator for a conformally coupled massless scalar field in the background of the three-dimensional black hole. The expectation value $\langle\varphi^2\rangle$ in the Hartle-Hawking state is obtained in the spacetime.


title: Direct detection of pure spin-current by x-ray pump-probe measurements

abstract: By synchronizing a microwave waveform with the synchrotron x-ray pulses, we use the ferromagnetic resonance (FMR) of the Py (Ni81Fe19) layer in a Py/Cu/Cu75Mn25/Cu/Co multilayer to pump a pure spin current into the Cu75Mn25 spacer layer, and then directly probe the spin current in the Cu75Mn25 layer by a time-resolved x-ray magnetic circular dichroism (XMCD). This element-specific pump-probe measurement unambiguously identifies the AC spin current in the Cu75Mn25 layer. In addition, phase resolved x-ray measurements reveal a characteristic bipolar phase behavior of the Co spins that is a fingerprint of spin-current driven spin precession.


title: Tangents, rectifiability, and corkscrew domains

abstract: In a recent paper, Cs\"ornyei and Wilson prove that curves in Euclidean space of $\sigma$-finite length have tangents on a set of positive $\mathscr{H}^{1}$-measure. They also show that a higher dimensional analogue of this result is not possible without some additional assumptions. In this note, we show that if $\Sigma\subseteq \mathbb{R}^{d+1}$ has the property that each ball centered on $\Sigma$ contains two large balls in different components of $\Sigma^{c}$ and $\Sigma$ has $\sigma$-finite $\mathscr{H}^{d}$-measure, then it has $d$-dimensional tangent points in a set of positive $\mathscr{H}^{d}$-measure. We also give shorter proofs that Semmes surfaces are uniformly rectifiable and, if $\Omega\subseteq \mathbb{R}^{d+1}$ is an exterior corkscrew domain whose boundary has locally finite $\mathscr{H}^{d}$-measure, one can find a Lipschitz subdomain intersecting a large portion of the boundary.


title: Computationally efficient formulation of relay operator for Preisach hysteresis modeling

abstract: An algebraic expression for the Preisach hysteron, which is a non-ideal (delayed) relay operator, is formulated for a computationally efficient real-time implementation. This allows representing the classical scalar Preisach hysteresis model as a summation of a large number of weighted hysterons which computation can be accomplished in parallel. The latter makes possible an efficient FPGA or ASIC realization of the scalar Preisach hysteresis model that can be useful for multiple applications. The signal flow which specifies the model implementation is provided in form of the block diagram. The proposed computation of Preisach hysterons, aggregated to the entire Preisach hysteresis model, is evaluated numerically and on a real-time hardware platform.


title: Deciding 4-colorability of planar triangulations

abstract: We show, without using the Four Color Theorem, that for each planar triangulation, the number of its proper vertex colorings by 4 colors is a determinant and thus can be calculated in a polynomial time. In particular, we can efficiently decide if the number is non-zero.


title: Algebraic bounds for heterogeneous site percolation on directed and undirected graphs

abstract: We analyze site percolation on directed and undirected graphs with site-dependent open-site probabilities. We construct upper bounds on cluster susceptibilities, vertex connectivity functions, and the expected number of simple open cycles through a chosen arc; separate bounds are given on finite and infinite (di)graphs. These produce lower bounds for percolation and uniqueness transitions in infinite (di)graphs, and for the formation of a giant component in finite (di)graphs. The bounds are formulated in terms of appropriately weighted adjacency and non-backtracking (Hashimoto) matrices. It turns out to be the uniqueness criterion that is most closely associated with an asymptotically vanishing probability of forming a giant strongly-connected component on a large finite (di)graph.


title: Algebraic identification of the effective connectivity of constrained geometric network models of neural signaling

abstract: Cellular neural circuit and networks consisting of interconnected neurons and glia are ulti- mately responsible for the information processing associated with information processing in the brain. While there are major efforts aimed at mapping the structural and (electro)physiological connectivity of brain networks, such as the White House BRAIN Initiative aimed at the devel- opment of neurotechnologies capable of high density neural recordings, theoretical and compu- tational methods for analyzing and making sense of all this data seem to be further behind. Here, we propose and provide a summary of an approach for calculating effective connectivity from experimental observations of neuronal network activity. The proposed method operates on network-level data, makes use of all relevant prior knowledge, such as dynamical models of individual cells in the network and the physical structural connectivity of the network, and is broadly applicable to large classes of biological and non-biological networks.


title: Quantum decoherence: a logical perspective

abstract: The so-called classical limit of quantum mechanics is generally studied in terms of the decoherence of the state operator that characterizes a system. This is not the only possible approach to decoherence. In previous works we have presented the possibility of studying the classical limit in terms of the decoherence of relevant observables of the system. On the basis of this approach, in this paper we introduce the classical limit from a logical perspective, by studying the way in which the logical structure of quantum properties corresponding to relevant observables acquires Boolean characteristics.


title: Phase Diagram of the Interacting Majorana Chain Model

abstract: The Hubbard chain and spinless fermion chain are paradigms of strongly correlated systems, very well understood using Bethe ansatz, Density Matrix Renormalization Group (DMRG) and field theory/renormalization group (RG) methods. They have been applied to one-dimensional materials and have provided important insights for understanding higher dimensional cases. Recently, a new interacting fermion model has been introduced, with possible applications to topological materials. It has a single Majorana fermion operator on each lattice site and interactions with the shortest possible range that involve 4 sites. We present a thorough analysis of the phase diagram of this model in one dimension using field theory/RG and DMRG methods. It includes a gapped supersymmetric region and a novel gapless phase with coexisting Luttinger liquid and Ising degrees of freedom. In addition to a first order transition, three critical points occur: tricritical Ising, Lifshitz and a novel generalization of the commensurate-incommensurate transition. We also survey various gapped phases of the system that arise when the translation symmetry is broken by dimerization and find both trivial and topological phases with 0, 1 and 2 Majorana zero modes bound to the edges of the chain with open boundary conditions.


title: Efficient computation of the Grunwald-Letnikov fractional diffusion derivative using adaptive time step memory

abstract: Computing numerical solutions to fractional differential equations can be computationally intensive due to the effect of non-local derivatives in which all previous time points contribute to the current iteration. In general, numerical approaches that depend on truncating part of the system history while efficient, can suffer from high degrees of error and inaccuracy. Here we present an adaptive time step memory method for smooth functions applied to the Grunwald-Letnikov fractional diffusion derivative. This method is computationally efficient and results in smaller errors during numerical simulations. Sampled points along the system history at progressively longer intervals are assumed to reflect the values of neighboring time points. By including progressively fewer points backward in time, a temporally weighted history is computed that includes contributions from the entire past of the system, maintaining accuracy, but with fewer points actually calculated, greatly improving computational efficiency.


title: Escape dynamics and fractal basin boundaries in Seyfert galaxies

abstract: The escape dynamics in a simple analytical gravitational model which describes the motion of stars in a Seyfert galaxy is investigated in detail. We conduct a thorough numerical analysis distinguishing between regular and chaotic orbits as well as between trapped and escaping orbits, considering only unbounded motion for several energy levels. In order to distinguish safely and with certainty between ordered and chaotic motion, we apply the Smaller ALingment Index (SALI) method. It is of particular interest to locate the escape basins through the openings around the collinear Lagrangian points $L_1$ and $L_2$ and relate them with the corresponding spatial distribution of the escape times of the orbits. Our exploration takes place both in the physical $(x,y)$ and in the phase $(x,\dot{x})$ space in order to elucidate the escape process as well as the overall orbital properties of the galactic system. Our numerical analysis reveals the strong dependence of the properties of the considered escape basins with the total orbital energy, with a remarkable presence of fractal basin boundaries along all the escape regimes. It was also observed, that for energy levels close to the critical escape energy the escape rates of orbits are large, while for much higher values of energy most of the orbits have low escape periods or they escape immediately to infinity. We also present evidence obtained through numerical simulations that our model can describe the formation and the evolution of the observed spiral structure in Seyfert galaxies. We hope our outcomes to be useful for a further understanding of the escape mechanism in galaxies with active nuclei.


title: On the Existence of Drifting Orbits for Non-Convex Hamiltonian Systems

abstract: We show the existence of drifting orbits for certain perturbations of non-convex Hamiltonian systems with several degrees of freedom. These orbits remain in the vicinity of resonant surfaces where the action variables can undergo changes $O(1)$ infinitely often although the size of perturbations $O(\epsilon)$ can be arbitrarily small. The first drifts occur in a period of time $O(1/\epsilon)$ and then reoccur with frequencies independent of $\epsilon$. We also perform numerical simulations to compare the effects of two conditions for instability in two four-dimensional examples with random parameters.


title: $C^1$-triangulations of semialgebraic sets

abstract: We show that every semialgebraic set admits a semialgebraic triangulation such that each closed simplex is $C^1$ differentiable. As an application, we give a straightforward definition of the integration $\int_X \omega$ over a compact semialgebraic subset $X$ of a differential form $\omega$ on an ambient algebraic manifold, that provides a significant simplification of the theory of semialgebraic singular chains and integrations. Our results hold over every (possibly non-archimedian) real closed field.


title: On the origin of the 11.3 micron unidentified infrared emission feature

abstract: The 11.3 $\mu$m emission feature is a prominent member of the family of unidentified infrared emission (UIE) bands and is frequently attributed to out-of-plane bending modes of polycyclic aromatic hydrocarbon (PAH) molecules. We have performed quantum mechanical calculations of 60 neutral PAH molecules and found that it is difficult to reconcile the observed astronomical feature with any or a mix of these PAH molecules. We have further analyzed the fitting of spectra of several astronomical objects by the NASA PAH database program and found that reasonable fittings to the observed spectra are only possible by including significant contributions from oxygen and/or magnesium containing molecules in the mix. A mixed of pure PAH molecules, even including units of different sizes, geometry and charged states, is unable to fit the astronomical spectra. Preliminary theoretical results on the vibrational spectra of simple molecules with mixed aromatic/aliphatic structures show that these structures have consistent bundles of vibrational modes and could be viable carriers of the UIE bands.


title: Electric-field guiding of magnetic skyrmions

abstract: We theoretically study equilibrium and dynamic properties of nanosized magnetic skyrmions in thin magnetic films with broken inversion symmetry, where electric field couples to magnetization via spin-orbit coupling. Based on a symmetry-based phenomenology and micromagnetic simulations we show that this electric-field coupling, via renormalizing the micromagnetic energy, modifies the equilibrium properties of the skyrmion. This change, in turn, results in a significant alteration of the current-induced skyrmion motion. Particularly, speed and direction of the skyrmion can be manipulated by designing a desired energy landscape electrically, which we describe within Thiele's analytical model and demonstrate in micromagnetic simulations including electric-field-controlled magnetic anisotropy. We additionally use this electric-field control to construct gates for controlling skyrmion motion exhibiting a transistor-like and multiplexer-like function. The proposed electric-field effect can thus provide a low energy electrical knob to extend the reach of information processing with skyrmions.


title: Generating admissible space-time meshes for moving domains in $d+1$-dimensions

abstract: In this paper we present a discontinuous Galerkin finite element method for the solution of the transient Stokes equations on moving domains. For the discretization we use an interior penalty Galerkin approach in space, and an upwind technique in time. The method is based on a decomposition of the space-time cylinder into finite elements. Our focus lies on three-dimensional moving geometries, thus we need to triangulate four dimensional objects. For this we will present an algorithm to generate $d+1$-dimensional simplex space-time meshes and we show under natural assumptions that the resulting space-time meshes are admissible. Further we will show how one can generate a four-dimensional object resolving the domain movement. First numerical results for the transient Stokes equations on triangulations generated with the newly developed meshing algorithm are presented.


title: On the Bhabha scattering for $z=2$ Lifshitz QED

abstract: In this paper we compute and discuss the differential cross section of the Bhabha scattering in the framework of the $z=2$ Lifshitz QED. We start by constructing the classical solutions for the fermionic fields, in particular the completeness relations, and also deriving the theory's propagators. Afterwards, we compute the photon exchange and pair annihilation contributions for the Bhabha's process, and upon the results we establish the magnitude of the theory's free parameter by looking for small deviations of the QED tree results.


title: Impurity Effects on Bound States in the Vortex Core of Topological S-wave Superconductor

abstract: We study the impurity effects on the Caroli-de Gennes-Matricon (CdGM) states, particularly on the level spacings in a vortex core in topological s-wave superconductor (SC) by two means, numerically and analytically. The topological s-wave SC belongs to the same class as a chiral p-wave SC and thus there are two inequivalent vortices in terms of any symmetry operation. We take into account this inequivalence and numerically calculate the scattering rates based on an improved version of Kopnin-Kravtsov (iKK) scheme, which enables us to treat the discrete levels in the presence of white-noise disorder. We also construct the Andreev equation for the topological s-wave SC and obtain the Andreev bound states analytically. We use a correspondence between the wave functions for the Bogoliubov-de Gennes equation and the Andreev equation in the iKK scheme and deduce the formula of scattering rates described by the wave function for the Andreev equation. With this formula, we discuss the origin of impurity scattering rates for CdGM states of topological s-wave SC and the dependence on the types of vortices related to the inequivalence.


title: Complete integrability of geodesic motion in Sasaki-Einstein toric $Y^{p,q}$ spaces

abstract: We construct explicitly the constants of motion for geodesics in the $5$-dimensional Sasaki-Einstein spaces $Y^{p,q}$. To carry out this task we use the knowledge of the complete set of Killing vectors and Killing-Yano tensors on these spaces. In spite of the fact that we generate a multitude of constants of motion, only five of them are functionally independent implying the complete integrability of geodesic flow on $Y^{p,q}$ spaces. In the particular case of the homogeneous Sasaki-Einstein manifold $T^{1,1}$ the integrals of motion have simpler forms and the relations between them are described in detail.


title: A wide diversity of 3D surfaces Generator using a new implicit function

abstract: We present in this paper a new family of implicit function for synthesizing a wide variety of 3D surfaces. The basis of this family consists of the usual functions that are: the function rectangular pulses, the function saw-tooth pulses, the function of triangular pulses, the staircase function and the power function. By combining these common functions, named constituent functions, in one implicit function and by varying some parameters of this function we can synthesize a wide variety of 3D surfaces with the possibility to set their deformations.


title: On the $\eta-\mu$/gamma and the $\lambda-\mu$/gamma Composite Distributions

abstract: This work is devoted to the formulation and derivation of the $\eta{-}\mu{/}$gamma and $\lambda{-}\mu{/}$gamma distributions which correspond to physical fading models. These distributions are composite and are based on the $\eta-\mu$ and $\lambda-\mu$ generalized multipath models, respectively, and the gamma shadowing model. Novel analytic expressions are derived for the corresponding envelope probability density functions. Importantly, the proposed models provide accurate characterisation of the simultaneous occurrence of multipath fading and shadowing effects which is achieved thanks to the remarkable flexibility offered by their parameters that render them capable of providing good fittings to experimental data associated with realistic communication scenarios. This is additionally justified by the fact that they include as special cases the widely known fading models such as Hoyt/gamma, Nakagami-m/gamma and Rayleigh/gamma. As a result, they can be meaningfully utilized in various analytical studies related to the performance evaluation of digital communications over composite multipath/shadowing fading channels.


title: Branching within branching I: The extinction problem

abstract: We consider a discrete-time host-parasite model for a population of cells which are colonized by proliferating parasites. The cell population grows like an ordinary Galton-Watson process, but in reflection of real biological settings the multiplication mechanisms of cells and parasites are allowed to obey some dependence structure. More precisely, the number of offspring produced by a mother cell determines the reproduction law of a parasite living in this cell and also the way the parasite offspring is shared into the daughter cells. In this article, we provide a formal introduction of this branching-within-branching model and then focus on the property of parasite extinction. We establish equivalent conditions for almost sure extinction of parasites, and find a strong relation of this event to the behavior of parasite multiplication along a randomly chosen cell line through the cell tree, which forms a branching process in random environment. In a second paper, the case when parasites survive is studied by proving limit results.


title: Optimal Dividend Strategies for Two Collaborating Insurance Companies

abstract: We consider a two-dimensional optimal dividend problem in the context of two insurance companies with compound Poisson surplus processes, who collaborate by paying each other's deficit when possible. We solve the stochastic control problem of maximizing the weighted sum of expected discounted dividend payments (among all admissible dividend strategies) until ruin of both companies, by extending results of univariate optimal control theory. In the case that the dividends paid by the two companies are equally weighted, the value function of this problem compares favorably with the one of merging the two companies completely. We identify this optimal value function as the smallest viscosity supersolution of the respective Hamilton-Jacobi-Bellman equation and provide an iterative approach to approximate it numerically. Curve strategies are identified as the natural analogue of barrier strategies in this two-dimensional context. A numerical example is given for which such a curve strategy is indeed optimal among all admissible dividend strategies, and for which this collaboration mechanism also outperforms the suitably weighted optimal dividend strategies of the two stand-alone companies.


title: Branching within branching II: Limit theorems

abstract: This continues work started in part I on a general branching-within-branching model for host-parasite co-evolution. Here we focus on asymptotic results for relevant processes in the case when parasites survive. In particular, limit theorems for the processes of contaminated cells and of parasites are established by using martingale theory and the technique of size-biasing. The results for both processes are of Kesten-Stigum type by including equivalent integrability conditions for the martingale limits to be positive with positive probability. The case when these conditions fail is also studied. For the process of contaminated cells, we show that a proper Heyde-Seneta norming exists such that the limit is nondegenerate.


title: Spatial adiabatic passage via interaction-induced band separation

abstract: The development of advanced quantum technologies and the quest for a deeper understanding of many-particle quantum mechanics requires control over the quantum state of interacting particles to a high degree of fidelity. However, the quickly increasing density of the spectrum, together with the appearance of crossings in time-dependent processes, makes any effort to control the system hard and resource intensive. Here we show that in trapped systems regimes can exist, in which isolated energy bands appear that allow to easily generalize known single-particle techniques. We demonstrate this for the well-known spatial adiabatic passage effect, which can control the center-of-mass state of atoms with high fidelity.


title: Global integration of the Schr\"odinger equation: a short iterative scheme within the wave operator formalism using discrete Fourier transforms

abstract: A global solution of the Schr\"odinger equation for explicitly time-dependent Hamiltonians is derived by integrating the non-linear differential equation associated with the time-dependent wave operator. A fast iterative solution method is proposed in which, however, numerous integrals over time have to be evaluated. This internal work is done using a numerical integrator based on Fast Fourier Transforms (FFT). The case of a transition between two potential wells of a model molecule driven by intense laser pulses is used as an illustrative example. This application reveals some interesting features of the integration technique. Each iteration provides a global approximate solution on grid points regularly distributed over the full time propagation interval. Inside the convergence radius, the complete integration is competitive with standard algorithms, especially when high accuracy is required.


title: Location Prediction of Social Images via Generative Model

abstract: The vast amount of geo-tagged social images has attracted great attention in research of predicting location using the plentiful content of images, such as visual content and textual description. Most of the existing researches use the text-based or vision-based method to predict location. There still exists a problem: how to effectively exploit the correlation between different types of content as well as their geographical distributions for location prediction. In this paper, we propose to predict image location by learning the latent relation between geographical location and multiple types of image content. In particularly, we propose a geographical topic model GTMI (geographical topic model of social image) to integrate multiple types of image content as well as the geographical distributions, In GTMI, image topic is modeled on both text vocabulary and visual feature. Each region has its own distribution over topics and hence has its own language model and vision pattern. The location of a new image is estimated based on the joint probability of image content and similarity measure on topic distribution between images. Experiment results demonstrate the performance of location prediction based on GTMI.


title: Towards the efficiency limits of silicon solar cells: how thin is too thin?

abstract: It is currently possible to fabricate crystalline silicon solar cells with the absorber thickness ranging from a few hundreds of micrometers (conventional wafer-based cells) to devices as thin as $1\,\mu\mathrm{m}$. In this work, we use a model single-junction solar cell to calculate the limits of energy conversion efficiency and estimate the optimal absorber thickness. The limiting efficiency for cells in the thickness range between 40 and $500\,\mu\mathrm{m}$ is very similar and close to 29%. In this regard, we argue that decreasing the thickness below around $40\,\mu\mathrm{m}$ is counter-productive, as it significantly reduces the maximum achievable efficiency, even when optimal light trapping is implemented. We analyse the roles of incomplete light trapping and extrinsic (bulk and surface) recombination mechanisms. For a reasonably high material quality, consistent with present-day fabrication techniques, the optimal thickness is always higher than a few tens of micrometers. We identify incomplete light trapping and parasitic losses as a major roadblock in improving the efficiency upon the current record of 25.6% for silicon solar cells. Finally, considering the main parameters that impact solar cell performance, we quantify the constraints and requirements for achieving a specified energy conversion efficiency, which is important for a proper design strategy of high efficiency silicon solar cells.


title: Hausdorff dimension of the graphs of the classical Weierstrass functions

abstract: We show that the graph of the classical Weierstrass function $\sum_{n=0}^\infty \lambda^n \cos (2\pi b^n x)$ has Hausdorff dimension $2+\log\lambda/\log b$, for every integer $b\ge 2$ and every $\lambda\in (1/b,1)$. Replacing $\cos(2\pi x)$ by a general non-constant $C^2$ periodic function, we obtain the same result under a further assumption that $\lambda b$ is close to $1$.