This study conducted a comprehensive analysis of trace element concentrations in the Upper Indus River Basin (UIRB), a glacier-fed region in the Western Himalayas (WH), aiming to discern their environmental and anthropogenic sources and implications. Despite limited prior data, 69 samples were collected in 2019 from diverse sources within the UIRB, including mainstream, tributaries, and groundwater, to assess trace element concentrations. Enrichment factor (EF) results and comparisons with regional and global averages suggest that rising levels of Zn, Cd, and As may pose safety concerns for drinking water quality. Advanced multivariate statistical techniques such as principal component analysis (PCA), absolute principal component scores (APCS-MLR), Monte Carlo simulation (MCS), etc were applied to estimate the associated human health hazards and also identified key sources of trace elements. The 95th percentile of the MCS results indicates that the estimated total cancer risk for children is significantly greater than (>1000 times) the USEPA\'s acceptable risk threshold of 1.0 × 10-6. The results classified most of the trace elements into two distinct groups: Group A (Li, Rb, Sr, U, Cs, V, Ni, TI, Sb, Mo, Ge), linked to geogenic sources, showed lower concentrations in the lower-middle river reaches, including tributaries and downstream regions. Group B (Pb, Nb, Cr, Zn, Be, Al, Th, Ga, Cu, Co), influenced by both geogenic and anthropogenic activities, exhibited higher concentrations near urban centers and midstream areas, aligning with increased municipal waste and agricultural activities. Furthermore, APCS-MLR source apportionment indicated that trace elements originated from natural geogenic processes, including rock-water interactions and mineral dissolution, as well as anthropogenic activities. These findings underscore the need for targeted measures to mitigate anthropogenic impacts and safeguard water resources for communities along the IRB and WH.

A Monte Carlo study of the compatibilization of A/B polymer blends has been performed using the bond fluctuation model. The considered compatibilizers are copolymer molecules composed of A and B blocks. Different types of copolymer structures have been included, namely, linear diblock and 4-block alternating copolymers, star block copolymers, miktoarm stars, and zipper stars. Zipper stars are composed of two arms of diblock copolymers arranged in alternate order (AB and BA) from the central unit, along with two homogeneous arms of A and B units. The compatibilization performance has been characterized by analyzing the equilibration of repulsion energy, the simulated scattering intensity obtained with opposite refractive indices for A and B, the profiles along a coordinate axis, the radial distribution functions, and the compatibilizer aggregation numbers. According to the results, linear alternate block copolymers, star block copolymers, and zipper stars exhibit significantly better compatibilization, with zipper stars showing slightly but consistently better performance.

BACKGROUND: Dosimetry-based personalized therapy was shown to have clinical benefits e.g. in liver selective internal radiation therapy (SIRT). Yet, there is no consensus about its introduction into clinical practice, mainly as Monte Carlo simulations (gold standard for dosimetry) involve massive computation time. We addressed the problem of computation time and tested a patch-based approach for Monte Carlo simulations for internal dosimetry to improve parallelization. We introduce a physics-inspired cropping layout for patch-based MC dosimetry, and compare it to cropping layouts of the literature as well as dosimetry using organ-S-values, and dose kernels, taking whole-body Monte Carlo simulations as ground truth. This was evaluated in five patients receiving Yttrium-90 liver SIRT.
RESULTS: The patch-based Monte Carlo approach yielded the closest results to the ground truth, making it a valid alternative to the conventional approach. Our physics-inspired cropping layout and mosaicking scheme yielded a voxel-wise error of < 2% compared to whole-body Monte Carlo in soft tissue, while requiring only ≈  10% of the time.
CONCLUSIONS: This work demonstrates the feasibility and accuracy of physics-inspired cropping layouts for patch-based Monte Carlo simulations.

Tacrolimus (TAC) has high pharmacokinetic (PK) variability during the early transplantation period. The relationships between whole-blood and intracellular TAC concentrations and clinical outcomes remain controversial. This study identifies the factors affecting the PK variability of TAC and characterizes the relationships between whole-blood and intracellular TAC concentrations. Data regarding whole-blood TAC concentrations of 1,787 samples from 215 renal transplant recipients (<90 days postoperative) across two centers and intracellular TAC concentrations (648 samples) digitized from previous studies were analyzed using nonlinear mixed-effects modeling. The effects of potential covariates were screened, and the distribution of whole-blood to intracellular TAC concentration ratios (RWB:IC) was estimated. The final model was evaluated using bootstrap, goodness of fit, and prediction-corrected visual predictive checks. The optimal dosing regimens and target ranges for each type of immune cell subsets were determined using Monte Carlo simulations. A two-compartment model adequately described the data, and the estimated mean TAC CL/F was 23.6 L·h-1 (relative standard error: 11.5 %). The hematocrit level, CYP3A5*3 carrier status, co-administration with Wuzhi capsules, and tapering prednisolone dose may contribute to the high variability of TAC PK variability during the early post-transplant period. The estimated RWB:IC of all TAC concentrations in peripheral blood mononuclear cells (PBMCs) was 4940, and inter-center variability of PBMCs was observed. The simulated TAC target range in PBMCs was 20.2-85.9 pg·million cells-1. Inter-center variability in intracellular concentrations should be taken into account in further analyses. TAC dosage adjustments can be guided based on PK/PD variability and simulated intracellular concentrations.

Phase behavior of nanoconfined fluids adsorbed in metal-organic frameworks is of paramount importance for the design of advanced materials for energy and gas storage, separations, electrochemical devices, sensors, and drug delivery, as well as for the pore structure characterization. Phase transformations in adsorbed fluids often involve long-lasting metastable states and hysteresis that has been well-documented in gas adsorption-desorption and nonwetting fluid intrusion-extrusion experiments. However, theoretical prediction of the observed nanophase behavior remains a challenging problem. The mesoscopic canonical, or mesocanonical, ensemble (MCE) is devised to study the nanophase behavior under conditions of controlled fluctuations to stabilize metastable and labile states. Here, we implement and apply the MCE Monte Carlo (MCEMC) simulation scheme to predict the origins of reversible and hysteric adsorption phase transitions in a series of practical MOF materials, including IRMOF-1, ZIF-412, UiO-66, Cu-BTC, IRMOF-74-V, VII, and IX. The MCEMC method, called the gauge cell method, allows to produce Van der Waals type isotherms with distinctive swings around the phase transition regions. The constructed isotherms determine the positions of phase equilibrium and spinodals, as well as the nucleation barriers separating metastable states. We demonstrate the unique capabilities of the MCEMC method in quantitative predictions of experimental observations compared with the conventional grand canonical and canonical ensemble simulations. The MCEMC method is implemented in the open-source RASPA and LAMMPS software packages and recommended for studies of adsorption behavior and pore structure characterization of MOFs and other nanoporous materials.

OBJECTIVE: This study aimed to assess the accuracy of patient-specific absorbed dose calculations for tumours and organs at risk in radiopharmaceutical therapy planning, utilizing hybrid planar-SPECT/CT imaging.
METHODS: Three Monte Carlo (MC) simulated digital patient phantoms were created, with time-activity data for mIBG labelled to I-123 (LEHR and ME collimators) and I-131 (HE collimator). The study assessed the accuracy of the mean absorbed doses for I-131-mIBG therapy treatment planning. Multiple planar whole-body (WB) images were simulated (between 1 to 72 h post-injection (p.i)). The geometric-mean image of the anterior and posterior WB images was calculated, with scatter and attenuation corrections applied. Time-activity curves were created for regions of interest over the liver and two tumours (diameters: 3.0 cm and 5.0 cm) in the WB images. A corresponding SPECT study was simulated at 24 h p.i and reconstructed using the OS-EM algorithm, incorporating scatter, attenuation, collimator-detector response, septal scatter and penetration corrections. MC voxel-based absorbed dose rate calculations used two image sets, (i) the activity distribution represented by the SPECT images and (ii) the activity distribution from the SPECT images distributed uniformly within the volume of interest. Mean absorbed doses were calculated considering photon and charged particle emissions, and beta emissions only. True absorbed doses were calculated by MC voxel-based dosimetry of the known activity distributions for reference.
RESULTS: Considering photon and charged particle emissions, mean absorbed dose accuracies across all three radionuclide-collimator combinations of 3.8 ± 5.5% and 0.1 ± 0.9% (liver), 5.2 ± 10.0% and 4.3 ± 1.7% (3.0 cm tumour) and 15.0 ± 5.8% and 2.6 ± 0.6% (5.0 cm tumour) were obtained for image set (i) and (ii) respectively. Considering charged particle emissions, accuracies of 2.7 ± 4.1% and 5.7 ± 0.7% (liver), 3.2 ± 10.2% and 9.1 ± 1.7% (3.0 cm tumour) and 13.6 ± 5.7% and 7.0 ± 0.6% (5.0 cm tumour) were obtained for image set (i) and (ii) respectively.
CONCLUSIONS: The hybrid WB planar-SPECT/CT method proved accurate for I-131-mIBG dosimetry, suggesting its potential for personalized treatment planning.

Singlet-to-triplet intersystem crossing (ISC) in organic molecules is intimately connected with their geometries: by modifying the molecular shape, symmetry selection rules pertaining to spin-orbit coupling can be partially relieved, leading to extra matrix elements for increased ISC. As an analog to this molecular design concept, the study finds that the lattice symmetry of supramolecular polymers also defines their triplet formation efficiencies. A supramolecular polymer self-assembled from weakly interacting molecules is considered. Its 2D oblique unit cell effectively renders it as a coplanar array of 1D molecular columns weakly bound to each other. Using momentum-resolved photoluminescence imaging in combination with Monte Carlo simulations, the study found that photogenerated charge carriers in the supramolecular polymer predominantly recombine as spin-uncorrelated carrier pairs through inter-column charge transfer states. This lattice-defined recombination pathway leads to a substantial triplet formation efficiency (≈60%) in the supramolecular polymer. These findings suggest that lattice symmetry of micro-/macroscopic structures relying on intermolecular interactions can be strategized for controlled triplet formation.

Autoregressive models in time series are useful in various areas. In this article, we propose a skew-t autoregressive model. We estimate its parameters using the expectation-maximization (EM) method and develop the influence methodology based on local perturbations for its validation. We obtain the normal curvatures for four perturbation strategies to identify influential observations, and then to assess their performance through Monte Carlo simulations. An example of financial data analysis is presented to study daily log-returns for Brent crude futures and investigate possible impact by the COVID-19 pandemic.

Researchers often estimate the association between the hazard of a time-to-event outcome and the characteristics of individuals and the clusters in which individuals are nested. Lin and Wei\'s robust variance estimator is often used with a Cox regression model fit to clustered data. Recently, alternative variance estimators have been proposed: the Fay-Graubard estimator, the Kauermann-Carroll estimator, and the Mancl-DeRouen estimator. Using Monte Carlo simulations, we found that, when fitting a marginal Cox regression model with both individual-level and cluster-level covariates: (i) in the presence of weak to moderate within-cluster homogeneity of outcomes, the Lin-Wei variance estimator can result in estimates of the SE with moderate bias when the number of clusters is fewer than 20-30, while in the presence of strong within-cluster homogeneity, it can result in biased estimation even when the number of clusters is as large as 100; (ii) when the number of clusters was less than approximately 20, the Fay-Graubard variance estimator tended to result in estimates of SE with the lowest bias; (iii) when the number of clusters exceeded approximately 20, the Mancl-DeRouen estimator tended to result in estimated standard errors with the lowest bias; (iv) the Mancl-DeRouen estimator used with a t-distribution tended to result in 95% confidence that had the best performance of the estimators; (v) when the magnitude of within-cluster homogeneity in outcomes was strong or very strong, all methods resulted in confidence intervals with lower than advertised coverage rates even when the number of clusters was very large.

Objective.There is an increasing interest in calculating and measuring linear energy transfer (LET) spectra in particle therapy in order to assess their impact in biological terms. As such, the accuracy of the particle fluence energy spectra becomes paramount. This study focuses on quantifying energy depositions of distinct proton, helium, carbon, and oxygen ion beams using a silicon pixel detector developed at CERN to determine LET spectra in silicon.Approach.While detection systems have been investigated in this pursuit, the scarcity of detectors capable of providing per-ion data with high spatial and temporal resolution remains an issue. This gap is where silicon pixel detector technology steps in, enabling online tracking of single-ion energy deposition. The used detector consisted of a 300µm thick silicon sensor operated in partial depletion.Main results.During post-processing, artifacts in the acquired signals were identified and methods for their corrections were developed. Subsequently, a correlation between measured and Monte Carlo-based simulated energy deposition distributions was performed, relying on a two-step recalibration approach based on linear and saturating exponential models. Despite the observed saturation effects, deviations were confined below 7% across the entire investigated range of track-averaged LET values in silicon from 0.77 keVµm-1to 93.16 keVµm-1.Significance.Simulated and measured mean energy depositions were found to be aligned within 7%, after applying artifact corrections. This extends the range of accessible LET spectra in silicon to clinically relevant values and validates the accuracy and reliability of the measurements. These findings pave the way towards LET-based dosimetry through an approach to translate these measurements to LET spectra in water. This will be addressed in a future study, extending functionality of treatment planning systems into clinical routine, with the potential of providing ion-beam therapy of utmost precision to cancer patients.