title: Non-supersymmetric Conifold

abstract: We find a new family of non-supersymmetric numerical solutions of IIB supergravity which are dual to the N=1 cascading "conifold" theory perturbed by certain combinations of relevant single trace and marginal double trace operators with non infinitesimal couplings. The SUSY is broken but the resulting ground states, and their gravity duals, remain stable, at least perturbatively.Despite the complicated field theory dynamics the gravity solutions have a simple structure. They feature the Ricci-flat non-Kahler metric on the deformed conifold and the imaginary self-dual three-form flux accompanied by a constant dilaton.


title: Neutrino-Nucleus Interactions around 1 GeV

abstract: The extraction of neutrino oscillation parameters often requires the unambiguous identification of the reaction mechanism to be used for the energy determination. The use of nuclear targets and of beams with broad energy distributions complicates this task because it requires the reliable description of many different elementary processes and the use of event generators to isolate and identify them. In this talk I briefly summarize our present understanding of this task.


title: NLO results for five, six and seven jets in electron-positron annihilation

abstract: We present next-to-leading order corrections in the leading colour approximation for jet rates in electron-positron annihilation up to seven jets. The results for the two-, three-, and four-jet rates agree with known results. The NLO jet rates have been known previously only up to five jets. The results for the six- and seven-jet rate are new. The results are obtained by a new and efficient method based on subtraction and numerical integration.


title: Failure of protection of Majorana based qubits against decoherence

abstract: Qubit realizations based on Majorana bound states have been considered promising candidates for quantum information processing which is inherently inert to decoherence. We put the underlying general arguments leading to this conjecture to the test from an open quantum system perspective. It turns out that, from a fundamental point of view, the Majorana qubit is as susceptible to decoherence as any local paradigm of a qubit.


title: HST-COS Observations of Hydrogen, Helium, Carbon and Nitrogen Emission from the SN 1987A Reverse Shock

abstract: We present the most sensitive ultraviolet observations of Supernova 1987A to date. Imaging spectroscopy from the Hubble Space Telescope-Cosmic Origins Spectrograph shows many narrow (dv \sim 300 km/s) emission lines from the circumstellar ring, broad (dv \sim 10 -- 20 x 10^3 km/s) emission lines from the reverse shock, and ultraviolet continuum emission. The high signal-to-noise (> 40 per resolution element) broad LyA emission is excited by soft X-ray and EUV heating of mostly neutral gas in the circumstellar ring and outer supernova debris. The ultraviolet continuum at \lambda > 1350A can be explained by HI 2-photon emission from the same region. We confirm our earlier, tentative detection of NV \lambda 1240 emission from the reverse shock and we present the first detections of broad HeII \lambda1640, CIV \lambda1550, and NIV] \lambda1486 emission lines from the reverse shock. The helium abundance in the high-velocity material is He/H = 0.14 +/- 0.06. The NV/H-alpha line ratio requires partial ion-electron equilibration (T_{e}/T_{p} \approx 0.14 - 0.35). We find that the N/C abundance ratio in the gas crossing the reverse shock is significantly higher than that in the circumstellar ring, a result that may be attributed to chemical stratification in the outer envelope of the supernova progenitor. The N/C abundance ratio may have been stratified prior to the ring expulsion, or this result may indicate continued CNO processing in the progenitor subsequent to the expulsion of the circumstellar ring.


title: Bars rejuvenating bulges? Evidence from stellar population analysis

abstract: We obtained stellar ages and metallicities via spectrum fitting for a sample of 575 bulges with spectra available from the Sloan Digital Sky Survey. The structural properties of the galaxies have been studied in detail in Gadotti (2009b) and the sample contains 251 bulges in galaxies with bars. Using the whole sample, where galaxy stellar mass distributions for barred and unbarred galaxies are similar, we find that bulges in barred and unbarred galaxies occupy similar loci in the age vs. metallicity plane. However, the distribution of bulge ages in barred galaxies shows an excess of populations younger than ~ 4 Gyr, when compared to bulges in unbarred galaxies. Kolmogorov-Smirnov statistics confirm that the age distributions are different with a significance of 99.94%. If we select sub-samples for which the bulge stellar mass distributions are similar for barred and unbarred galaxies, this excess vanishes for galaxies with bulge mass log M < 10.1 M_Sun while for more massive galaxies we find a bimodal bulge age distribution for barred galaxies only, corresponding to two normal distributions with mean ages of 10.4 and 4.7 Gyr. We also find twice as much AGN among barred galaxies, as compared to unbarred galaxies, for low-mass bulges. By combining a large sample of high quality data with sophisticated image and spectral analysis, we are able to find evidence that the presence of bars affect the mean stellar ages of bulges. This lends strong support to models in which bars trigger star formation activity in the centers of galaxies.


title: Prediction of Orientational Phase Transition in Boron Carbide

abstract: The assessed binary phase diagram of boron-carbon exhibits a single alloy phase designated "B4C" with rhombohedral symmetry occupying a broad composition range that falls just short of the nominal carbon content of 20%. As this composition range is nearly temperature independent, the phase diagram suggests a violation of the third law of thermodynamics, which typically requires compounds to achieve a definite stoichiometry at low temperatures. By means of first principles total energy calculations we predict the existence of two stoichiometric phases at T=0K: one of composition B4C with monoclinic symmetry; the other of composition B13C2 with rhombohedral symmetry. Using statistical mechanics to extend to finite temperatures, we demonstrate that the monoclinic phase reverts to the rhombohedral phase above T=600K, along with a slight reduction on carbon content.


title: Quantization via Empirical Divergence Maximization

abstract: Empirical divergence maximization (EDM) refers to a recently proposed strategy for estimating f-divergences and likelihood ratio functions. This paper extends the idea to empirical vector quantization where one seeks to empirically derive quantization rules that maximize the Kullback-Leibler divergence between two statistical hypotheses. We analyze the estimator's error convergence rate leveraging Tsybakov's margin condition and show that rates as fast as 1/n are possible, where n equals the number of training samples. We also show that the Flynn and Gray algorithm can be used to efficiently compute EDM estimates and show that they can be efficiently and accurately represented by recursive dyadic partitions. The EDM formulation have several advantages. First, the formulation gives access to the tools and results of empirical process theory that quantify the estimator's error convergence rate. Second, the formulation provides a previously unknown derivation for the Flynn and Gray algorithm. Third, the flexibility it affords allows one to avoid a small-cell assumption common in other approaches. Finally, we illustrate the potential use of the method through an example.


title: Adam Adamandy Kocha\'nski's approximations of \pi: reconstruction of the algorithm

abstract: In his 1685 paper "Observationes cyclometricae" published in Acta Eruditorum, Adam Adamandy Kocha\'nski presented an approximate ruler-and-compass construction for rectification of the circle. It is not generally known that the first part of this paper included an interesting sequence of rational approximations of \pi. Kocha\'nski gave only a partial explanation of the algorithm used to produce these approximations, while promising to publish details at a later time, which has never happened. We reconstruct the complete algorithm and discuss some of its properties. We also argue that Kocha\'nski was very close to discovery of continued fractions and convergents of \pi.


title: Transverse Momentum Dependent Fragmentation and Quark Distribution Functions from the NJL-jet Model

abstract: Using the model of Nambu and Jona-Lasinio to provide a microscopic description of both the structure of the nucleon and of the quark to hadron elementary fragmentation functions, we investigate the transverse momentum dependence of the unpolarized quark distributions in the nucleon and of the quark to pion and kaon fragmentation functions. The transverse momentum dependence of the fragmentation functions is determined within a Monte Carlo framework, with the notable result that the average $P_\perp^2$ of the produced kaons is significantly larger than that of the pions. We also find that $<P_\perp^2>$ has a sizable $z$ dependence, in contrast with the naive Gaussian ansatz for the fragmentation functions. Diquark correlations in the nucleon give rise to a non-trivial flavor dependence in the unpolarized transverse momentum dependent quark distribution functions. The $<k_T^2>$ of the quarks in the nucleon are also found to have a sizable $x$ dependence. Finally, these results are used as input to a Monte Carlo event generator for semi-inclusive deep inelastic scattering (SIDIS), which is used to determine the average transverse momentum squared of the produced hadrons measured in SIDIS, namely $<P_T^2>$. Again we find that the average $P_T^2$ of the produced kaons in SIDIS is significantly larger than that of the pions and in each case $\la P_T^2 \ra$ has a sizable $z$ dependence.


title: Light curve modelling for mutual transits

abstract: In this paper we describe an algorithm and deduce the related mathematical formulae that allows the computation of observed fluxes in stellar and planetary systems with arbitrary number of bodies being part of a transit or occultation event. The presented method does not have any limits or constraints for the geometry and can be applied for almost all of the available limb darkening models as well. As a demonstration, we apply this scheme to gather information for the orbital inclinations from multiple transiting planetary systems in cases when mutual transits occur. We also show here that these mutual events constrain inclinations unambiguously, yielding a complete picture for the whole system.


title: Spontaneous chiral symmetry breaking in the Tayler instability

abstract: The chiral symmetry breaking properties of the Tayler instability are discussed. Effective amplitude equations are determined in one case. This model has three free parameters that are determined numerically. Comparison with chiral symmetry breaking in biochemistry is made.


title: Higher Derivative Gravity from the Universal Renormalization Group Machine

abstract: We study the renormalization group flow of higher derivative gravity, utilizing the functional renormalization group equation for the average action. Employing a recently proposed algorithm, termed the universal renormalization group machine, for solving the flow equation, all the universal features of the one-loop beta-functions are recovered. While the universal part of the beta-functions admits two fixed points, we explicitly show that the existence of one of them depends on the choice of regularization scheme, indicating that it is most probably unphysical.


title: Matrix Models for Supersymmetric Chern-Simons Theories with an ADE Classification

abstract: We consider N=3 supersymmetric Chern-Simons (CS) theories that contain product U(N) gauge groups and bifundamental matter fields. Using the matrix model of Kapustin, Willett and Yaakov, we examine the Euclidean partition function of these theories on an S^3 in the large N limit. We show that the only such CS theories for which the long range forces between the eigenvalues cancel have quivers which are in one-to-one correspondence with the simply laced affine Dynkin diagrams. As the A_n series was studied in detail before, in this paper we compute the partition function for the D_4 quiver. The D_4 example gives further evidence for a conjecture that the saddle point eigenvalue distribution is determined by the distribution of gauge invariant chiral operators. We also see that the partition function is invariant under a generalized Seiberg duality for CS theories.


title: Introduction to stability conditions

abstract: These are notes of a course given at the 'school on moduli spaces' at the Newton Institute in January 2011. The abstract theory of stability conditions (due to Bridgeland and Douglas) on abelian and triangulated categories is developed via tilting and t-structures. Special emphasis is put on the bounded derived category of coherent sheaves on smooth projective varieties (in particular for curves and K3 surfaces). The lectures were targeted at an audience with little prior knowledge of triangulated categories and stability conditions but with an interest in vector bundles on curves.


title: The necessity of entanglement and the equivalency of Bell's theorem with the second law of thermodynamics

abstract: We demonstrate that both Wigner's form of Bell's inequalities as well as a form of the second law of thermodynamics, as manifest in Carath\'{e}odory's principle, can be derived from the same simple experimental and statistical mechanical assumptions combined with the trivial behavior of integers. This suggests that Bell's theorem is merely a well-disguised statement of the second law. It also suggests that entanglement is necessary for quantum theory to be in full accord with the second law and thus builds on the results of Wiesniak, Vedral, and Brukner \cite{Marcin-Wiesniak:2008fv} who showed it was necessary for consistency with the third law.


title: Global structure and kinematics of stellar haloes in cosmological hydrodynamic simulations

abstract: We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamical simulations to study the global structure and kinematics of stellar spheroids of Milky Way mass disc galaxies. Font et al. have recently demonstrated that these simulations are able to successfully reproduce the satellite luminosity functions and the metallicity and surface brightness profiles of the spheroids of the Milky Way and M31. A key to the success of the simulations is a significant contribution to the spheroid from stars that formed in situ. While the outer halo is dominated by accreted stars, stars formed in the main progenitor of the galaxy dominate at r < ~30 kpc. In the present study we show that this component was primarily formed in a proto-disc at high redshift and was subsequently liberated from the disc by dynamical heating associated with mass accretion. As a consequence of its origin, the in situ component of the spheroid has different kinematics (namely net prograde rotation with respect to the disc) than that of the spheroid component built from the disruption of satellites. In addition, the in situ component has a flattened distribution, that is due in part to its rotation. We make comparisons with measurements of the shape and kinematics of local galaxies, including the Milky Way and M31, and stacked observations of more distant galaxies. We find that the simulated disc galaxies have spheroids of the correct shape (oblate with a median axis ratio of ~0.6 at radii of < ~30 kpc, but note there is significant system-to-system scatter in this quantity) and that the kinematics show evidence for two components (due to in situ vs. accreted), as observed. Our findings therefore add considerable weight to the importance of dissipative processes in the formation of stellar haloes and to the notion of a 'dual stellar halo'.


title: Global Lipschitz regularizing effects for linear and nonlinear parabolic equations

abstract: In this paper we prove global bounds on the spatial gradient of viscosity solutions to second order linear and nonlinear parabolic equations in $(0,T) \times \R^N$. Our assumptions include the case that the coefficients be both unbounded and with very mild local regularity (possibly weaker than the Dini continuity), the estimates only depending on the parabolicity constant and on the modulus of continuity of coefficients (but not on their $L^{\infty}$-norm). Our proof provides the analytic counterpart to the probabilistic proof used in Priola and Wang (J. Funct. Anal. 2006) to get this type of gradient estimates in the linear case. We actually extend such estimates to the case of possibly unbounded data and solutions as well as to the case of nonlinear operators including Bellman-Isaacs equations. We investigate both the classical regularizing effect (at time $t>0$) and the possible conservation of Lipschitz regularity from $t=0$, and similarly we prove global H\"older estimates under weaker assumptions on the coefficients. The estimates we prove for unbounded data and solutions seem to be new even in the classical case of linear equations with bounded and H\"older continuous coefficients. Applications to Liouville type theorems are also given in the paper. Finally, we compare in an appendix the analytic and the probabilistic approach discussing the analogy between the doubling variables method of viscosity solutions and the probabilistic coupling method.


title: An invitation to algebraic topological string theory

abstract: The purpose of this note is to provide a short invitation to the universal algebraic approach to topological string theory. In the first section we make an attempt to explain the origin of this approach and how it fits into the bigger picture of full string theory, while in the second half of this note we will introduce the relevant notions in more detail and discuss some of our main results on bulk-deformed open topological string amplitudes.


title: Near Linear-Work Parallel SDD Solvers, Low-Diameter Decomposition, and Low-Stretch Subgraphs

abstract: We present the design and analysis of a near linear-work parallel algorithm for solving symmetric diagonally dominant (SDD) linear systems. On input of a SDD $n$-by-$n$ matrix $A$ with $m$ non-zero entries and a vector $b$, our algorithm computes a vector $\tilde{x}$ such that $\norm[A]{\tilde{x} - A^+b} \leq \vareps \cdot \norm[A]{A^+b}$ in $O(m\log^{O(1)}{n}\log{\frac1\epsilon})$ work and $O(m^{1/3+\theta}\log \frac1\epsilon)$ depth for any fixed $\theta > 0$. The algorithm relies on a parallel algorithm for generating low-stretch spanning trees or spanning subgraphs. To this end, we first develop a parallel decomposition algorithm that in polylogarithmic depth and $\otilde(|E|)$ work, partitions a graph into components with polylogarithmic diameter such that only a small fraction of the original edges are between the components. This can be used to generate low-stretch spanning trees with average stretch $O(n^{\alpha})$ in $O(n^{1+\alpha})$ work and $O(n^{\alpha})$ depth. Alternatively, it can be used to generate spanning subgraphs with polylogarithmic average stretch in $\otilde(|E|)$ work and polylogarithmic depth. We apply this subgraph construction to derive a parallel linear system solver. By using this solver in known applications, our results imply improved parallel randomized algorithms for several problems, including single-source shortest paths, maximum flow, minimum-cost flow, and approximate maximum flow.


title: Effective Medium Theory of Filamentous Triangular Lattice

abstract: We present an effective medium theory that includes bending as well as stretching forces, and we use it to calculate mechanical response of a diluted filamentous triangular lattice. In this lattice, bonds are central-force springs, and there are bending forces between neighboring bonds on the same filament. We investigate the diluted lattice in which each bond is present with a probability $p$. We find a rigidity threshold $p_b$ which has the same value for all positive bending rigidity and a crossover characterizing bending-, stretching-, and bend-stretch coupled elastic regimes controlled by the central-force rigidity percolation point at $p_{\textrm{CF}} \simeq 2/3$ of the lattice when fiber bending rigidity vanishes.


title: New Method for 3D Shape Retrieval

abstract: The recent technological progress in acquisition, modeling and processing of 3D data leads to the proliferation of a large number of 3D objects databases. Consequently, the techniques used for content based 3D retrieval has become necessary. In this paper, we introduce a new method for 3D objects recognition and retrieval by using a set of binary images CLI (Characteristic level images). We propose a 3D indexing and search approach based on the similarity between characteristic level images using Hu moments for it indexing. To measure the similarity between 3D objects we compute the Hausdorff distance between a vectors descriptor. The performance of this new approach is evaluated at set of 3D object of well known database, is NTU (National Taiwan University) database.


title: Electronic beam shifts in monolayer graphene superlattice

abstract: Electronic analogue of generalized Goos-H\"{a}nchen shifts is investigated in the monolayer graphene superlattice with one-dimensional periodic potentials of square barriers. It is found that the lateral shifts for the electron beam transmitted through the monolayer graphene superlattice can be negative as well as positive near the band edges of zero-$\bar{k}$ gap, which are different from those near the band edges of Bragg gap. These negative and positive beam shifts have close relation to the Dirac point. When the condition $q_A d_A= -q_B d_B= m \pi$ ($m=1,2,3...$) is satisfied, the beam shifts can be controlled from negative to positive when the incident energy is above the Dirac point, and vice versa. In addition, the beam shifts can be greatly enhanced by the defect mode inside the zero-$\bar{k}$ gap. These intriguing phenomena can be verified in a relatively simple optical setup, and have potential applications in the graphene-based electron wave devices.


title: Electronic band gap and transport in Fibonacci quasi-periodic graphene superlattice

abstract: We investigate electronic band gap and transport in Fibonacci quasi-periodic graphene superlattice. It is found that such structure can possess a zero-$\bar{k}$ gap which exists in all Fibonacci sequences. Different from Bragg gap, zero-$\bar{k}$ gap associated with Dirac point is less sensitive to the incidence angle and lattice constants. The defect mode appeared inside the zero-$\bar{k}$ gap has a great effect on transmission, conductance and shot noise, which can be applicable to control the electron transport.


title: Propagating Lyapunov Functions to Prove Noise--induced Stabilization

abstract: We investigate an example of noise-induced stabilization in the plane that was also considered in (Gawedzki, Herzog, Wehr 2010) and (Birrell, Herzog, Wehr 2011). We show that despite the deterministic system not being globally stable, the addition of additive noise in the vertical direction leads to a unique invariant probability measure to which the system converges at a uniform, exponential rate. These facts are established primarily through the construction of a Lyapunov function which we generate as the solution to a sequence of Poisson equations. Unlike a number of other works, however, our Lyapunov function is constructed in a systematic way, and we present a meta-algorithm we hope will be applicable to other problems. We conclude by proving positivity properties of the transition density by using Malliavin calculus via some unusually explicit calculations.


title: On fixed points of a generalized multidimensional affine recursion

abstract: Let $G$ be a multiplicative subsemigroup of the general linear group $\Gl(\mathbb{R}^d)$ which consists of matrices with positive entries such that every column and every row contains a strictly positive element. Given a $G$--valued random matrix $A$, we consider the following generalized multidimensional affine equation R\stackrel{\mathcal{D}}{=}\sum_{i=1}^N A_iR_i+B, where $N\ge2$ is a fixed natural number, $A_1,...,A_N$ are independent copies of $A$, $B\in\mathbb{R}^d$ is a random vector with positive entries, and $R_1,...,R_N$ are independent copies of $R\in\mathbb{R}^d$, which have also positive entries. Moreover, all of them are mutually independent and $\stackrel{\mathcal{D}}{=}$ stands for the equality in distribution. We will show with the aid of spectral theory developed by Guivarc'h and Le Page and Kesten's renewal theorem, that under appropriate conditions, there exists $\chi>0$ such that $\P(\{<R, u>>t\})\asymp t^{-\chi},$ as $t\to\8$, for every unit vector $u\in\mathbb{S}^{d-1}$ with positive entries.


title: The DPOSS II distant compact group survey: the EMMI-NTT spectroscopic sample

abstract: This paper presents the results of the redshift survey of 138 candidate compact groups from the DPOSS II catalog (Iovino et al., 2003; de Carvalho et al. 2005), which extends the available redshift range of spectroscopically confirmed compact groups of galaxies to z~0.2. The aims of the survey are to confirm group membership via spectroscopic redshift information, to measure the characteristic parameters of the confirmed groups, namely mass, radius, luminosity, velocity dispersion and crossing time, and to compare them with those of nearby compact groups. Using available information from the literature, we also studied the surrounding group environment and searched for additional, previously unknown, group members, or larger scale structures to whom the group might be associated. Of the 138 observed groups, 96 had three or more concordant galaxies, i.e. a 70% success rate. Of these 96, 62 are isolated on the sky, while the other 34 are close on the sky to a larger scale structure. The remaining objects turned out to be couple of pairs or chance projection of galaxies on the sky. Group environment and the characteristics parameters (mass, crossing time, velocity dispersion) are evaluated and discussed.


title: IVOA Recommendation: Observation Data Model Core Components and its Implementation in the Table Access Protocol Version 1.0

abstract: This document defines the core components of the Observation data model that are necessary to perform data discovery when querying data centers for observations of interest. It exposes use-cases to be carried out, explains the model and provides guidelines for its implementation as a data access service based on the Table Access Protocol (TAP). It aims at providing a simple model easy to understand and to implement by data providers that wish to publish their data into the Virtual Observatory. This interface integrates data modeling and data access aspects in a single service and is named ObsTAP. It will be referenced as such in the IVOA registries. There will be a separate document to cover the full Observation data model. In this document, the Observation Data Model Core Components (ObsCoreDM) defines the core components of queryable metadata required for global discovery of observational data. It is meant to allow a single query to be posed to TAP services at multiple sites to perform global data discovery without having to understand the details of the services present at each site. It defines a minimal set of basic metadata and thus allows for a reasonable cost of implementation by data providers. The combination of the ObsCoreDM with TAP is referred to as an ObsTAP service. As with most of the VO Data Models, ObsCoreDM makes use of STC, Utypes, Units and UCDs. The ObsCoreDM can be serialized as a VOTable. ObsCoreDM can make reference to more complete data models such as ObsProvDM (the Observation Provenance Data Model, to come), Characterisation DM, Spectrum DM or Simple Spectral Line Data Model (SSLDM).


title: Separation of a Slater determinant wave function with a neck structure into spatially localized subsystems

abstract: A method to separate a Slater determinant wave function with a two-center neck structure into spatially localized subsystems is proposed, and its potential applications are presented. An orthonormal set of spatially localized single-particle wave functions is obtained by diagonalizing the coordinate operator for the major axis of a necked system. Using the localized single-particle wave functions, the wave function of each subsystem is defined. Therefore, defined subsystem wave functions are used to obtain density distributions, mass centers, and energies of subsystems. The present method is applied to separations of Margenau--Brink cluster wave functions of $\alpha + \alpha$, $^{16}$O + $^{16}$O, and $\alpha + ^{16}$O into their subsystems, and also to separations of antisymmetrized molecular dynamics wave functions of $^{10}$Be into $\alpha$ + $^6$He subsystems. The method is simple and applicable to the separation of general Slater determinant wave functions that have neck structures into subsystem wave functions.


title: An Interpretation of "Superluminal Neutrino" Compatible with Relativity in the Framework of Standard Model

abstract: According to the measurement of muon-neutrino experiment done by the OPERA collaboration, the speed of high-energy neutrino exceeds that of light in vacuum by 25ppm. Assuming that this result is correct, a possible resolution of the dilemma between it and the validity of relativity is proposed without changing the framework of the Standard Model of elementary particles. The essential idea is based on a possible resolution, proposed previously, of the color confinement problem of quantum chromodynamics.


title: Asymptotic Behavior for a Nonlocal Diffusion Equation in Domains with Holes

abstract: The paper deals with the asymptotic behavior of solutions to a non-local diffusion equation, $u_t=J*u-u:=Lu$, in an exterior domain, $\Omega$, which excludes one or several holes, and with zero Dirichlet data on $\mathbb{R}^N\setminus\Omega$. When the space dimension is three or more this behavior is given by a multiple of the fundamental solution of the heat equation away from the holes. On the other hand, if the solution is scaled according to its decay factor, close to the holes it behaves like a function that is $L$-harmonic, $Lu=0$, in the exterior domain and vanishes in its complement. The height of such a function at infinity is determined through a matching procedure with the multiple of the fundamental solution of the heat equation representing the outer behavior. The inner and the outer behavior can be presented in a unified way through a suitable global approximation.


title: A generalised porosity formalism for isotropic and anisotropic effective opacity and its effects on X-ray line attenuation in clumped O star winds

abstract: We present a generalised formalism for treating the porosity-associated reduction in continuum opacity that occurs when individual clumps in a stochastic medium become optically thick. We consider geometries resulting in either isotropic or anisotropic effective opacity, and, in addition to an idealised model in which all clumps have the same local overdensity and scale, we also treat an ensemble of clumps with optical depths set by Markovian statistics. This formalism is then applied to the specific case of bound-free absorption of X-rays in hot star winds, a process not directly affected by clumping in the optically thin limit. We find that the Markov model gives surprisingly similar results to those found previously for the single clump model, suggesting that porous opacity is not very sensitive to details of the assumed clump distribution function. Further, an anisotropic effective opacity favours escape of X-rays emitted in the tangential direction (the `venetian blind' effect), resulting in a 'bump' of higher flux close to line centre as compared to profiles computed from isotropic porosity models. We demonstrate how this characteristic line shape may be used to diagnose the clump geometry, and we confirm previous results that for optically thick clumping to significantly influence X-ray line profiles, very large porosity lengths, defined as the mean free path between clumps, are required. Moreover, we present the first X-ray line profiles computed directly from line-driven instability simulations using a 3-D patch method, and find that porosity effects from such models also are very small. This further supports the view that porosity has, at most, a marginal effect on X-ray line diagnostics in O stars, and therefore that these diagnostics do indeed provide a good `clumping insensitive' method for deriving O star mass-loss rates.


title: Coronal ejections from convective spherical shell dynamos

abstract: We present a three-dimensional model of rotating convection combined with a simplified corona in spherical coordinates. The motions in the convection zone generate a large-scale magnetic field that gets sporadically ejected into the outer layers above. Our model corona is approximately isothermal, but it includes density stratification due to gravity.