Optimal Lewenstein-Sanpera Decomposition for some Biparatite Systems

It is shown that for a given bipartite density matrix and by choosing a suitable separable set (instead of product set) on the separable-entangled boundary, optimal Lewenstein-Sanpera (L-S) decomposition can be obtained via optimization for a generic entangled density matrix. Based on this, We obtain optimal L-S decomposition for some bipartite systems such as $2\otimes 2$ and $2\otimes 3$ Bell decomposable states, generic two qubit state in Wootters basis, iso-concurrence decomposable states, states obtained from BD states via one parameter and three parameters local operations and classical communications (LOCC), $d\otimes d$ Werner and isotropic states, and a one parameter $3\otimes 3$ state. We also obtain the optimal decomposition for multi partite isotropic state. It is shown that in all $2\otimes 2$ systems considered here the average concurrence of the decomposition is equal to the concurrence. We also show that for some $2\otimes 3$ Bell decomposable states the average concurrence of the decomposition is equal to the lower bound of the concurrence of state presented recently in [Buchleitner et al, quant-ph/0302144], so an exact expression for concurrence of these states is obtained. It is also shown that for $d\otimes d$ isotropic state where decomposition leads to a separable and an entangled pure state, the average I-concurrence of the decomposition is equal to the I-concurrence of the state. Keywords: Quantum entanglement, Optimal Lewenstein-Sanpera decomposition, Concurrence, Bell decomposable states, LOCC} PACS Index: 03.65.UdComment: 31 pages, Late

Improved Energy Resolution Measurements of Electron Precipitation Observed During an IPDP‐Type EMIC Event

High energy resolution DEMETER satellite observations from the Instrument for the Detection of Particle (IDP) are analyzed during an electromagnetic ion cyclotron (EMIC)-induced electron precipitation event. Analysis of an Interval Pulsation with Diminishing Periods (IPDP)-type EMIC wave event, using combined satellite observations to correct for incident proton contamination, detected an energy precipitation spectrum ranging from ∼150 keV to ∼1.5 MeV. While inconsistent with many theoretical predictions of >1 MeV EMIC-induced electron precipitation, the finding is consistent with an increasing number of experimentally observed events detected using lower resolution integral channel measurements on the POES, FIREBIRD, and ELFIN satellites. Revised and improved DEMETER differential energy fluxes, after correction for incident proton contamination shows that they agree to within 40% in peak flux magnitude, and 85 keV (within 40%) for the energy at which the peak occurred as calculated from POES integral channel electron precipitation measurements. This work shows that a subset of EMIC waves found close to the plasmapause, that is, IPDP-type rising tone events, can produce electron precipitation with peak energies substantially below 1 MeV. The rising tone features of IPDP EMIC waves, along with the association with the high cold plasma density regime, and the rapidly varying electron density gradients of the plasmapause may be an important factor in the generation of such low energy precipitation, co-incident with a high energy tail. Our work highlights the importance of undertaking proton contamination correction when using the high-resolution DEMETER particle measurements to investigate EMIC-driven electron precipitation

hp-adaptive Galerkin Time Stepping Methods for Nonlinear Initial Value Problems

This work is concerned with the derivation of an a posteriori error estimator for Galerkin approximations to nonlinear initial value problems with an emphasis on finite-time existence in the context of blow-up. The structure of the derived estimator leads naturally to the development of both h and hp versions of an adaptive algorithm designed to approximate the blow-up time. The adaptive algorithms are then applied in a series of numerical experiments, and the rate of convergence to the blow-up time is investigated

Modeling Field Line Curvature Scattering Loss of 1–10 MeV Protons During Geomagnetic Storms

The proton radiation belt contains high fluxes of adiabatically trapped protons varying in energy from ∼one to hundreds of megaelectron volts (MeV). At large radial distances, magnetospheric field lines become stretched on the nightside of Earth and exhibit a small radius of curvature RC near the equator. This leads protons to undergo field line curvature (FLC) scattering, whereby changes to the first adiabatic invariant accumulate as field strength becomes nonuniform across a gyroorbit. The outer boundary of the proton belt at a given energy corresponds to the range of magnetic L shell over which this transition to nonadiabatic motion takes place, and is sensitive to the occurrence of geomagnetic storms. In this work, we first find expressions for nightside equatorial RC and field strength Be as functions of Dst and L* to fit the TS04 field model. We then apply the Tu et al. (2014, https://doi.org/10.1002/2014ja019864) condition for nonadiabatic onset to solve the outer boundary L*, and refine our expression for RC to achieve agreement with Van Allen Probes observations of 1–50 MeV proton flux over the 2014–2018 era. Finally, we implement this nonadiabatic onset condition into the British Antarctic Survey proton belt model (BAS-PRO) to solve the temporal evolution of proton fluxes at L ≤ 4. Compared with observations, BAS-PRO reproduces storm losses due to FLC scattering, but there is a discrepancy in mid-2017 that suggests a ∼5 MeV proton source not accounted for. Our work sheds light on outer zone proton belt variability at 1–10 MeV and demonstrates a useful tool for real-time forecasting

Transcriptomic Profiling of DNA Damage Response in Patient-Derived Glioblastoma Cells before and after Radiation and Temozolomide Treatment

Glioblastoma is a highly aggressive, invasive and treatment-resistant tumour. The DNA damage response (DDR) provides tumour cells with enhanced ability to activate cell cycle arrest and repair treatment-induced DNA damage. We studied the expression of DDR, its relationship with standard treatment response and patient survival, and its activation after treatment. The transcriptomic profile of DDR pathways was characterised within a cohort of isocitrate dehydrogenase (IDH) wild-type glioblastoma from The Cancer Genome Atlas (TCGA) and 12 patient-derived glioblastoma cell lines. The relationship between DDR expression and patient survival and cell line response to temozolomide (TMZ) or radiation therapy (RT) was assessed. Finally, the expression of 84 DDR genes was examined in glioblastoma cells treated with TMZ and/or RT. Although distinct DDR cluster groups were apparent in the TCGA cohort and cell lines, no significant differences in OS and treatment response were observed. At the gene level, the high expression of ATP23, RAD51C and RPA3 independently associated with poor prognosis in glioblastoma patients. Finally, we observed a substantial upregulation of DDR genes after treatment with TMZ and/or RT, particularly in RT-treated glioblastoma cells, peaking within 24 h after treatment. Our results confirm the potential influence of DDR genes in patient outcome. The observation of DDR genes in response to TMZ and RT gives insight into the global response of DDR pathways after adjuvant treatment in glioblastoma, which may have utility in determining DDR targets for inhibition

Recirculating Flows Involving Short Fiber Suspensions: Numerical Difficulties and Efficient Advanced Micro-Macro Solvers

Numerical modelling of non-Newtonian flows usually involves the coupling between equations of motion characterized by an elliptic character, and the fluid constitutive equation, which defines an advection problem linked to the fluid history. There are different numerical techniques to treat the hyperbolic advection equations. In non-recirculating flows, Eulerian discretizations can give a convergent solution within a short computing time. However, the existence of steady recirculating flow areas induces additional difficulties. Actually, in these flows neither boundary conditions nor initial conditions are known. In this paper we compares different advanced strategies (some of them recently proposed and extended here for addressing complex flows) when they are applied to the solution of the kinetic theory description of a short fiber suspension fluid flows