Next-to-Minimal SOFTSUSY

We describe an extension to the SOFTSUSY program that provides for the calculation of the sparticle spectrum in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), where a chiral superfield that is a singlet of the Standard Model gauge group is added to the Minimal Supersymmetric Standard Model (MSSM) fields. Often, a Z3 symmetry is imposed upon the model. SOFTSUSY can calculate the spectrum in this case as well as the case where general Z3 violating (denoted as ) terms are added to the soft supersymmetry breaking terms and the superpotential. The user provides a theoretical boundary condition for the couplings and mass terms of the singlet. Radiative electroweak symmetry breaking data along with electroweak and CKM matrix data are used as weak-scale boundary conditions. The renormalisation group equations are solved numerically between the weak scale and a high energy scale using a nested iterative algorithm. This paper serves as a manual to the NMSSM mode of the program, detailing the approximations and conventions used.

Image-based biomechanical assessment of vertebral body and intervertebral disc in the human lumbar spine

Life expectancy continuously increases but our society faces age-related conditions. Among musculoskeletal diseases, osteoporosis associated with risk of vertebral fracture and degenerative intervertebral disc (IVD) are painful pathologies responsible for tremendous healthcare costs. Hence, reliable diagnostic tools are necessary to plan a treatment or follow up its efficacy. Yet, radiographic and MRI techniques, respectively clinical standards for evaluation of bone strength and IVD degeneration, are unspecific and not objective. Increasingly used in biomedical engineering, CT-based finite element (FE) models constitute the state-of-art for vertebral strength prediction. However, as non-invasive biomechanical evaluation and personalised FE models of the IVD are not available, rigid boundary conditions (BCs) are applied on the FE models to avoid uncertainties of disc degeneration that might bias the predictions. Moreover, considering the impact of low back pain, the biomechanical status of the IVD is needed as a criterion for early disc degeneration. Thus, the first FE study focuses on two rigid BCs applied on the vertebral bodies during compression test of cadaver vertebral bodies, vertebral sections and PMMA embedding. The second FE study highlights the large influence of the intervertebral disc’s compliance on the vertebral strength, damage distribution and its initiation. The third study introduces a new protocol for normalisation of the IVD stiffness in compression, torsion and bending using MRI-based data to account for its morphology. In the last study, a new criterion (Otsu threshold) for disc degeneration based on quantitative MRI data (axial T2 map) is proposed. The results show that vertebral strength and damage distribution computed with rigid BCs are identical. Yet, large discrepancies in strength and damage localisation were observed when the vertebral bodies were loaded via IVDs. The normalisation protocol attenuated the effect of geometry on the IVD stiffnesses without complete suppression. Finally, the Otsu threshold computed in the posterior part of annulus fibrosus was related to the disc biomechanics and meet objectivity and simplicity required for a clinical application. In conclusion, the stiffness normalisation protocol necessary for consistent IVD comparisons and the relation found between degeneration, mechanical response of the IVD and Otsu threshold lead the way for non-invasive evaluation biomechanical status of the IVD. As the FE prediction of vertebral strength is largely influenced by the IVD conditions, this data could also improve the future FE models of osteoporotic vertebra.

Relating in situ gas measurements to the surface outgassing properties of cometary nuclei

The sensitivity of the gas flow field to changes in different initial conditions has been studied for the case of a highly simplified cometary nucleus model. The nucleus model simulated a homogeneously outgassing sphere with a more active ring around an axis of symmetry. The varied initial conditions were the number density of the homogeneous region, the surface temperature, and the composition of the flow (varying amounts of H2O and CO2) from the active ring. The sensitivity analysis was performed using the Polynomial Chaos Expansion (PCE) method. Direct Simulation Monte Carlo (DSMC) was used for the flow, thereby allowing strong deviations from local thermal equilibrium. The PCE approach can be used to produce a sensitivity analysis with only four runs per modified input parameter and allows one to study and quantify non-linear responses of measurable parameters to linear changes in the input over a wide range. Hence the PCE allows one to obtain a functional relationship between the flow field properties at every point in the inner coma and the input conditions. It is for example shown that the velocity and the temperature of the background gas are not simply linear functions of the initial number density at the source. As probably expected, the main influence on the resulting flow field parameter is the corresponding initial parameter (i.e. the initial number density determines the background number density, the temperature of the surface determines the flow field temperature, etc.). However, the velocity of the flow field is also influenced by the surface temperature while the number density is not sensitive to the surface temperature at all in our model set-up. Another example is the change in the composition of the flow over the active area. Such changes can be seen in the velocity but again not in the number density. Although this study uses only a simple test case, we suggest that the approach, when applied to a real case in 3D, should assist in identifying the sensitivity of gas parameters measured in situ by, for example, the Rosetta spacecraft to the surface boundary conditions and vice versa.

Brand Authenticity: An Integrative Framework and Measurement Scale

This article develops an integrative framework of the concept of perceived brand authenticity (PBA) and sheds light on PBA’s (1) measurement, (2) drivers, (3) consequences, as well as (4) an underlying process of its effects and (5) boundary conditions. A multi-phase scale development process resulted in a 15-item PBA scale to measure its four dimensions of credibility, integrity, symbolism, and continuity. PBA is influenced by indexical, existential, and iconic cues, whereby the latter’s influence is moderated by consumers’ level of marketing skepticism. Results also suggest that PBA drives brand choice likelihood through self-congruence for consumers high in self-authenticity.

A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

North Atlantic Eddy-Driven Jet in interglacial and glacial winter climates

The atmospheric westerly flow in the North Atlantic (NA) sector is dominated by atmospheric waves or eddies generating via momentum flux convergence, the so-called eddy-driven jet. The position of this jet is variable and shows for the present-day winter climate three preferred latitudinal states: a northern, central, and southernposition in the NA. Here, the authors analyze the behavior of the eddy-driven jet under different glacial and interglacial boundary conditions using atmosphere–land-only simulations with the CCSM4 climate model. As state-of-the-art climate models tend to underestimate the trimodality of the jet latitude, the authors apply a bias correction and successfully extract the trimodal behavior of the jet within CCSM4. The analysis shows that during interglacial times (i.e., the early Holocene and the Eemian) the preferred jet positions are rather stable and the observed multimodality is the typical interglacial character of the jet. During glacial times, the jet is strongly enhanced, its position is shifted southward, and the trimodal behavior vanishes. This is mainly due to the presence of the Laurentide ice sheet (LIS). The LIS enhances stationary waves downstream, thereby accelerating and displacing the NA eddy-driven jet by anomalous stationary momentum flux convergence. Additionally, changes in the transient eddy activity caused by topography changes as well as other glacial boundary conditions lead to an acceleration of the westerly winds over the southern NA at the expenseof more northernareas. Consequently, bothstationaryand transient eddiesfoster the southward shift of the NA eddy-driven jet during glacial winter times.

Fate of accidental symmetries of the relativistic hydrogen atom in a spherical cavity

The non-relativistic hydrogen atom enjoys an accidental SO(4) symmetry, that enlarges the rotational SO(3) symmetry, by extending the angular momentum algebra with the Runge–Lenz vector. In the relativistic hydrogen atom the accidental symmetry is partially lifted. Due to the Johnson–Lippmann operator, which commutes with the Dirac Hamiltonian, some degeneracy remains. When the non-relativistic hydrogen atom is put in a spherical cavity of radius R with perfectly reflecting Robin boundary conditions, characterized by a self-adjoint extension parameter γ, in general the accidental SO(4) symmetry is lifted. However, for R=(l+1)(l+2)a (where a is the Bohr radius and l is the orbital angular momentum) some degeneracy remains when γ=∞ or γ = 2/R. In the relativistic case, we consider the most general spherically and parity invariant boundary condition, which is characterized by a self-adjoint extension parameter. In this case, the remnant accidental symmetry is always lifted in a finite volume. We also investigate the accidental symmetry in the context of the Pauli equation, which sheds light on the proper non-relativistic treatment including spin. In that case, again some degeneracy remains for specific values of R and γ.

Quantifying differences in land use emission estimates implied by definition discrepancies

The quantification of CO2 emissions from anthropogenic land use and land use change (eLUC) is essential to understand the drivers of the atmospheric CO2 increase and to inform climate change mitigation policy. Reported values in synthesis reports are commonly derived from different approaches (observation-driven bookkeeping and process-modelling) but recent work has emphasized that inconsistencies between methods may imply substantial differences in eLUC estimates. However, a consistent quantification is lacking and no concise modelling protocol for the separation of primary and secondary components of eLUC has been established. Here, we review differences of eLUC quantification methods and apply an Earth System Model (ESM) of Intermediate Complexity to quantify them. We find that the magnitude of effects due to merely conceptual differences between ESM and offline vegetation model-based quantifications is ~ 20 % for today. Under a future business-as-usual scenario, differences tend to increase further due to slowing land conversion rates and an increasing impact of altered environmental conditions on land-atmosphere fluxes. We establish how coupled Earth System Models may be applied to separate secondary component fluxes of eLUC arising from the replacement of potential C sinks/sources and the land use feedback and show that secondary fluxes derived from offline vegetation models are conceptually and quantitatively not identical to either, nor their sum. Therefore, we argue that synthesis studies should resort to the "least common denominator" of different methods, following the bookkeeping approach where only primary land use emissions are quantified under the assumption of constant environmental boundary conditions.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p ≤ 0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Non-self-adjoint graphs

On finite metric graphs we consider Laplace operators, subject to various classes of non-self-adjoint boundary conditions imposed at graph vertices. We investigate spectral properties, existence of a Riesz basis of projectors and similarity transforms to self-adjoint Laplacians. Among other things, we describe a simple way to relate the similarity transforms between Laplacians on certain graphs with elementary similarity transforms between matrices defining the boundary conditions.

Pressure from the vacuum of confined spinor matter

Charged spinor matter field is quantized in a spatial region bounded by two parallel neutral plates. The most general set of boundary conditions ensuring the confinement of matter within the plates is considered. We study a response of the vacuum of the confined matter to the background uniform magnetic field which is directed orthogonally to the plates. It is proven that, in the case of a sufficiently strong magnetic field, the vacuum pressure onto the plates is positive and independent of the boundary condition, as well as of the distance between the plates.

Influence of quantized massive matter fields on the Casimir effect

Charged massive matter fields of spin-0 and spin- 1/2 are quantized in the presence of an external uniform magnetic field in a spatial region bounded by two parallel plates. The most general set of boundary conditions at the plates, that is required by mathematical consistency and the self-adjointness of the Hamiltonian operator, is employed. The vacuum fluctuations of the matter field in the case of the magnetic field orthogonal to the plates are analyzed, and it is shown that the pressure from the vacuum onto the plates is positive and independent of the boundary condition, as well as of the distance between the plates. Possibilities of the detection of this new-type Casimir effect are discussed. Read More: http://www.worldscientific.com/doi/10.1142/S0217732315500996

Casimir effect with quantized charged spinor matter in background magnetic field

We study the influence of a background uniform magnetic field and boundary conditions on the vacuum of a quantized charged spinor matter field confined between two parallel neutral plates; the magnetic field is directed orthogonally to the plates. The admissible set of boundary conditions at the plates is determined by the requirement that the Dirac Hamiltonian operator be self-adjoint. It is shown that, in the case of a sufficiently strong magnetic field and a sufficiently large separation of the plates, the generalized Casimir force is repulsive, being independent of the choice of a boundary condition, as well as of the distance between the plates. The detection of this effect seems to be feasible in the foreseeable future.

Hemodynamics analyses in treated and untreated carotid arteries of the same patient: A preliminary study based on three patient cases

Carotid atherosclerotic disease is highly related to cerebrovascular events. Carotid endarterectomy is the common operation method to treat this disease. In this study, hemodynamics analyses are performed on the carotid arteries in three patients, whose right carotid artery had been treated by carotid endarterectomy and the left carotid artery remained untreated. Flow and loading conditions are compared between these treated and untreated carotid arteries and evaluation of the operative results is discussed. Patient-specific models are reconstructed from MDCT data. Intraoperative ultrasound flow measurements are performed on the treated carotid arteries and the obtained data are used as the boundary conditions of the models and the validations of the computational results. Finite volume method is employed to solve the transport equations and the flow and loading conditions of the models are reported. The results indicate that: (i) in two of the three patients, the internal-to-external flow rate ratio in the untreated carotid artery is larger than that in the treated one, and the average overall flow split ratio by summing up the data of both the left and right carotid arteries is about 2.15; (ii) in the carotid bulb, high wall shear stress occurs at the bifurcation near the external carotid artery in all of the cases without hard plaques; (iii) the operated arteries present low time-averaged wall shear stress at the carotid bulb, especially for the treated arteries with patch technique, indicating the possibility of the recurrence of stenosis; (iv) high temporal gradient of wall shear stress (>35 Pa/s) is shown in the narrowing regions along the vessels; and (v) in the carotid arteries without serious stenosis, the maximum velocity magnitude during mid-diastole is 32~37% of that at systolic peak, however, in the carotid artery with 50% stenosis by hard plaques, this value is nearly doubled (64%). The computational work quantifies flow and loading distributions in the treated and untreated carotid arteries of the same patient, contributing to evaluation of the operative results and indicating the recurrent sites of potential atheromatous plaques.

hp-DGFEM for second-order mixed elliptic problems polyhedra

We prove exponential rates of convergence of hp-version discontinuous Galerkin (dG) interior penalty finite element methods for second-order elliptic problems with mixed Dirichlet-Neumann boundary conditions in axiparallel polyhedra. The dG discretizations are based on axiparallel, σ-geometric anisotropic meshes of mapped hexahedra and anisotropic polynomial degree distributions of μ-bounded variation. We consider piecewise analytic solutions which belong to a larger analytic class than those for the pure Dirichlet problem considered in [11, 12]. For such solutions, we establish the exponential convergence of a nonconforming dG interpolant given by local L 2 -projections on elements away from corners and edges, and by suitable local low-order quasi-interpolants on elements at corners and edges. Due to the appearance of non-homogeneous, weighted norms in the analytic regularity class, new arguments are introduced to bound the dG consistency errors in elements abutting on Neumann edges. The non-homogeneous norms also entail some crucial modifications of the stability and quasi-optimality proofs, as well as of the analysis for the anisotropic interpolation operators. The exponential convergence bounds for the dG interpolant constructed in this paper generalize the results of [11, 12] for the pure Dirichlet case.