Source author record

Zhen Tian

Zhen Tian appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2026arXiv

HEDP: A Hybrid Energy-Distance Prompt-based Framework for Domain Incremental Learning

Domain Incremental Learning is a critical scenario that requires models to continuously adapt to new data domains without retraining. However, domain shifts often cause severe performance degradation. To address this, we propose Hybrid Energy-Distance Prompt, a domain-incremental framework inspired by Helmholtz free energy. HEDP introduces an energy regularization loss to enhance the separability of domain representations and a hybrid energy-distance weighted mechanism that fuses energy-based and distance-based cues to improve domain selection and generalization. Experiments on multiple benchmarks, including CORe50, show that HEDP achieves superior performance on unseen domains with a 2.57\% accuracy gain, effectively mitigating catastrophic forgetting and enhancing open-world adaptability. Our code is \href{https://github.com/dannis97500/HEDP/}{available here}.

preprint2022arXiv

RecBole 2.0: Towards a More Up-to-Date Recommendation Library

In order to support the study of recent advances in recommender systems, this paper presents an extended recommendation library consisting of eight packages for up-to-date topics and architectures. First of all, from a data perspective, we consider three important topics related to data issues (i.e., sparsity, bias and distribution shift), and develop five packages accordingly: meta-learning, data augmentation, debiasing, fairness and cross-domain recommendation. Furthermore, from a model perspective, we develop two benchmarking packages for Transformer-based and graph neural network (GNN)-based models, respectively. All the packages (consisting of 65 new models) are developed based on a popular recommendation framework RecBole, ensuring that both the implementation and interface are unified. For each package, we provide complete implementations from data loading, experimental setup, evaluation and algorithm implementation. This library provides a valuable resource to facilitate the up-to-date research in recommender systems. The project is released at the link: https://github.com/RUCAIBox/RecBole2.0.

preprint2020arXiv

A plan quality control method of treatment planning for Gamma Knife radiosurgery

With many variables to adjust, conventional manual forward planning for Gamma Knife (GK) radiosurgery is very complicated and cumbersome. The resulting plan quality heavily depends on planners skills, experiences and devoted efforts, and varies significantly among cases, planners, and institutions. Quality control for GK planning is desired to consistently provide high-quality plan to each patient. In this study, we proposed a quality control method for GK planning by building a database of high-quality GK plans. Patient anatomy was described by target volume, target shape complexity, and spatial relationship between target and nearby organs, which determine GK planning difficulty level. Plan quality was evaluated using target coverage, selectivity, intermediate dose spillage, maximum dose to 0.1 cc of brainstem, mean dose of ipsilateral cochlea, and beam-on time. When a new plan is created, a high-quality plan that has the most similar target volume size and shape complexity will be identified from the database. A model has also been built to predict the dose to brainstem and cochlea based on their overlap volume histograms. The identified reference plan and the predicted organ dose will help planners to make quality control decisions accordingly. To validate this method, we have built a database for vestibular schwannoma, which are considered to be challenging for GK planning due to the irregularly-shaped target and its proximity to brainstem and cochlea. Five cases were tested, among which one case was considered to be of high quality and four cases had a lower plan quality than prediction. These four cases were replanned and got substantially improved. Our results have demonstrated the efficacy of our proposed quality control method. This method may also be used as a plan quality prediction method to facilitate the development of automatic treatment planning for GK radiosurgery.

preprint2020arXiv

An institutional study on plan quality and variation of manual forward planning for Gamma Knife radiosurgery for vestibular schwannoma

Due to the complexity and cumbersomeness of Gamma Knife (GK) manual forward planning, the quality of the resulting treatment plans heavily depends on the planners skill, experience and the amount of effort devoted to plan development. Hence, GK plan quality may vary significantly among institutions and planners, and even for a same planner at different cases. This is particularly a concern for challenging cases with complicated geometry, such as vestibular schwannoma cases. The purpose of this retrospective study is to investigate the plan quality and variation in the manually forward planned, clinically acceptable GK treatment plans of 22 previous vestibular schwannoma cases. Considering the impacts of different patient geometry and different trade-offs among the planning objectives in GK planning, it is difficult to objectively assess the plan quality across different cases. To reduce these confounding factors on plan quality assessment, we employed our recently developed multiresolution-level inverse planning algorithm to generate a golden plan for each case, which is expected to be on or close to the pareto surface with a similar trade-off as used in the manual plan. The plan quality of the manual plan is then quantified in terms of its deviation from the golden plan. A scoring criterion between 0-100 was designed to calculate a final score for each manual plan to simplify our analysis. Large quality variation was observed in these 22 cases, with two cases having a score lower than 75, three cases scoring between 80 and 85, two cases between 85 and 90, eight cases between 90 and 95, and seven cases higher than 95. Inter- and intra- planner variability was also observed in our study. This large variation in GK manual planning deserves high attention, and merits further investigation on how to reduce the variation in GK treatment plan quality.

preprint2020arXiv

Deep learning-based Real-time Volumetric Imaging for Lung Stereotactic Body Radiation Therapy: A Proof of Concept Study

Due to the inter- and intra- variation of respiratory motion, it is highly desired to provide real-time volumetric images during the treatment delivery of lung stereotactic body radiation therapy (SBRT) for accurate and active motion management. In this proof-of-concept study, we propose a novel generative adversarial network integrated with perceptual supervision to derive instantaneous volumetric images from a single 2D projection. Our proposed network, named TransNet, consists of three modules, i.e., encoding, transformation and decoding modules. Rather than only using image distance loss between the generated 3D images and the ground truth 3D CT images to supervise the network, perceptual loss in feature space is integrated into loss function to force the TransNet to yield accurate lung boundary. Adversarial supervision is also used to improve the realism of generated 3D images. We conducted a simulation study on 20 patient cases, who had received lung SBRT treatments in our institution and undergone 4D-CT simulation, and evaluated the efficacy and consistency of our method for four different projection angles, i.e., 0, 30, 60 and 90 degree. For each 3D CT image set of a breathing phase, we simulated its 2D projections at these angles.Then for each projection angle, a patient's 3D CT images of 9 phases and the corresponding 2D projection data were used for training, with the remaining phase used for testing. The mean absolute error, normalized MAE, peak signal-to-noise ratio and structural similarity index metric achieved by our method are 99.3 HU, 0.032, 23.4 dB and 0.949, respectively. These results demonstrate the feasibility and efficacy of our 2D-to-3D method for lung cancer patients, which provides a potential solution for in-treatment real-time on-board volumetric imaging for accurate dose delivery to ensure the effectiveness of lung SBRT treatment.

preprint2020arXiv

Interlayer Decoupling in 30° Twisted Bilayer Graphene Quasicrystal

Stacking order has strong influence on the coupling between the two layers of twisted bilayer graphene (BLG), which in turn determines its physical properties. Here, we report the investigation of the interlayer coupling of the epitaxially grown single-crystal 30° twisted BLG on Cu(111) at the atomic scale. The stacking order and morphology of BLG is controlled by a rationally designed two-step growth process, that is, the thermodynamically controlled nucleation and kinetically controlled growth. The crystal structure of the 30°-twisted bilayer graphene (30°-tBLG) is determined to have the quasicrystal like symmetry. The electronic properties and interlayer coupling of the 30°-tBLG is investigated using scanning tunneling microscopy (STM) and spectroscopy (STS). The energy-dependent local density of states (DOS) with in-situ electrostatic doping shows that the electronic states in two graphene layers are decoupled near the Dirac point. A linear dispersion originated from the constituent graphene monolayers is discovered with doubled degeneracy. This study contributes to controlled growth of twist-angle-defined BLG, and provides insights of the electronic properties and interlayer coupling in this intriguing system.

preprint2020arXiv

Potassium Isotope Compositions of Carbonaceous and Ordinary Chondrites: Implications on the Origin of Volatile Depletion in the Early Solar System

Solar system materials are variably depleted in moderately volatile elements (MVEs) relative to the proto-solar composition. To address the origin of this MVE depletion, we conducted a systematic study of high-precision K isotopic composition on 16 carbonaceous chondrites (CCs) of types CM1-2, CO3, CV3, CR2, CK4-5 and CH3 and 28 ordinary chondrites (OCs) covering petrological types 3 to 6 and chemical groups H, L, and LL. We observed significant overall K isotope (delta41K) variations (-1.54 to 0.70 permil). The K isotope compositions of CCs are largely higher than the Bulk Silicate Earth (BSE) value, whereas OCs show typically lower values than BSE. Neither CCs nor OCs show resolvable correlations between K isotopes and chemical groups, petrological types, shock levels, exposure ages, fall or find occurrence, or terrestrial weathering. The lack of a clear trend between K isotopes and K content indicates that the K isotope fractionations were decoupled from the relative elemental K depletions. The range of K isotope variations in the CCs is consistent with a four-component (chondrule, refractory inclusion, matrix and water) mixing model that is able to explain the bulk elemental and isotopic compositions of the main CC groups, but requires a fractionation in K isotopic compositions in chondrules. We propose that the major control of the isotopic compositions of group averages is condensation or vaporization in nebular environments that is preserved in the compositional variation of chondrules. Parent-body processes (aqueous alteration, thermal metamorphism, and metasomatism) can mobilize K and affect the K isotopes in individual samples. In the case of the OCs, the full range of K isotopic variations can only be explained by the combined effects of the size and relative abundances of chondrules, parent-body aqueous and thermal alteration.

preprint2020arXiv

Potassium Isotopic Compositions of Enstatite Meteorites

Enstatite chondrites and aubrites are meteorites that show the closest similarities to the Earth in many isotope systems that undergo mass-independent and mass-dependent isotope fractionations. Due to the analytical challenges to obtain high-precision K isotopic compositions in the past, potential differences in K isotopic compositions between enstatite meteorites and the Earth remained uncertain. We report the first high-precision K isotopic compositions of eight enstatite chondrites and four aubrites and find that there is a significant variation of K isotopic compositions among enstatite meteorites (from -2.34 permil to -0.18 permil). However, K isotopic compositions of nearly all enstatite meteorites scatter around the Bulk Silicate Earth (BSE) value. The average K isotopic composition of the eight enstatite chondrites (-0.47 +/- 0.57 permil) is indistinguishable from the BSE value (-0.48 +/- 0.03 permil), thus further corroborating the isotopic similarity between Earth' building blocks and enstatite meteorite precursors. We found no correlation of K isotopic compositions with the chemical groups, petrological types, shock degrees, and terrestrial weathering conditions; however, the variation of K isotopes among enstatite meteorite can be attributed to the parent body processing. Our sample of the main group aubrite MIL 13004 is exceptional and has an extremely light K isotopic composition (delta 41K= -2.34 +/- 0.12 permil). We attribute this unique K isotopic feature to the presence of abundant djerfisherite inclusions in our sample because this K-bearing sulfide mineral is predicted to be enriched in 39K during equilibrium exchange with silicates.

preprint2019arXiv

A preliminary study on a multi-resolution-level inverse planning algorithm for Gamma Knife radiosurgery

Manual forward planning for GK radiosurgery is complicated and time-consuming, particularly for cases with large or irregularly shaped targets. Inverse planning eases GK planning by solving an optimization problem. However, due to the vast search space, most inverse planning algorithms have to decouple the planning process to isocenter preselection and sector duration optimization. This sequential scheme does not necessarily lead to optimal isocenter locations and hence optimal plans. In this study, we attempt to optimize the isocenter positions, beam shapes and durations simultaneously by proposing a multi-resolution-level (MRL) strategy to handle the large-scale GK optimization problem. In our approach, several rounds of optimizations were performed with a progressively increased spatial resolution for isocenter candidate selection. The isocenters selected from last round and their neighbors on a finer resolution were used as new isocenter candidates for next round of optimization. After plan optimization, shot sequencing was performed to group the optimized sectors to deliverable shots supported by GK treatment units. We have tested our algorithm on 6 GK cases previously treated in our institution (2 meningioma cases, 3 cases with single metastasis and 1 case with 6 metastases). Compared with manual planning, achieving same coverage and similar selectivity, our algorithm improved the gradient index averagely from 3.1 to 2.9 and reduced the maximum dose of brainstem from 8.0Gy to 5.6Gy. The average beam-on time was also reduced by from 103.8 mins to 87.4 mins. Our method was also compared with the inverse planning algorithm provided in Leksell GammaPlan planning system, and outperformed it with better plan quality for all the 6 cases.This preliminary study has demonstrated the effectiveness and feasibility of our MRL inverse planning approach for GK radiosurgery.

preprint2019arXiv

Electron and positron spectra in the three dimensional spatial-dependent propagation model

The spatial-dependent propagation model has been successfully used to explain diverse observational phenomena, including the spectral hardening of cosmic-ray nuclei above $200$ GV, the large-scale dipole anisotropy and the diffusive gamma distribution. In this work, we further apply the spatial-dependent propagation model to both electrons and positrons. To account for the excess of positrons above $10$ GeV, an additional local source is introduced. And we also consider a more realistic spiral distribution of background sources. We find that due to the gradual hardening above $10$ GeV, the hardening of electron spectrum above tens of GeV can be explained in the SDP model and both positron and electron spectra less than TeV energies could be naturally described. The spatial-dependent propagation with spiral-distributed sources could conforms with the total electron spectrum in the whole energy. Meanwhile compared with the conventional model, the spatial-dependent propagation with spiral-distributed sources could produce larger background positron flux, so that the multiplier of background positron flux is $1.42$, which is much smaller than the required value by the conventional model. Thus the shortage of background positron flux could be solved. Furthermore we compute the anisotropy of electron under spatial-dependent propagation model, which is well below the observational limit of Fermi-LAT experiment.

preprint2016arXiv

2D SnS: a phosphorene analogue with strong in-plane electronic anisotropy

We study the anisotropic electronic properties of 2D SnS, an analogue of phosphorene, grown by physical vapor transport. With transmission electron microscope and polarized Raman spectroscopy, we identify the zigzag and armchair directions of the as-grown 2D crystals. 2D SnS field-effect transistors with a cross-Hall-bar structure are fabricated. They show heavily hole-doped (~10$^{19}$ cm$^{-3}$) conductivity with strong in-plane anisotropy. At room temperature the mobility along the zigzag direction exceeds 20 cm$^{2}$V$^{-1}$s$^{-1}$, which can be up to 1.7 times of that in the armchair direction. This strong anisotropy is then explained by the effective-mass ratio along the two directions and agrees well with previous theoretical predictions. Temperature-dependent carrier density is used to find out the acceptor energy level to be ~45 meV above the valence band maximum. This value matches with a calculated defect level of 42 meV for Sn vacancies, indicating that Sn deficiency is the main cause of the p-type conductivity.

preprint2016arXiv

Field-effect transistors of high-mobility few-layer SnSe$_{2}$

We report the transport properties of mechanically exfoliated few-layer SnSe$_{2}$ flakes, whose mobility is found with four probe measurements to be ~ 85 cm$^{2}$V$^{-1}$s$^{-1}$ at 300 K, higher than those of the majority of few-layer transitional metal dichalcogenides (TMDs). The mobility increases strongly with decreased temperature, indicating a phonon limited transport. The conductivity of the semiconducting SnSe$_{2}$ shows a metallic behavior, which is explained by two competing factors involving the different temperature dependence of mobility and carrier density. The Fermi level is found to be 87 meV below the conduction band minima (CBM) at 300 K and 12 meV below the CBM at 78 K, resulting from a heavy n-type doping. Previous studies have found SnSe$_{2}$ field-effect transistors (FETs) to be very difficult to turn off. We find the limiting factor to be the flake thickness compared with the maximum depletion width. With fully depleted devices, we are able to achieve a current on-off ratio of ~10$^{5}$. These results demonstrate the great potential of SnSe2 as a two dimensional (2D) semiconducting material and are helpful for our understanding of other heavily doped 2D materials.

preprint2015arXiv

An Analytic Linear Accelerator Source Model for Monte Carlo Dose Calculations. I. Model Representation and Construction

Monte Carlo (MC) simulation is considered as the most accurate method for radiation dose calculations. Accuracy of a source model for a linear accelerator is critical for the overall dose calculation accuracy. In this paper, we presented an analytical source model that we recently developed for GPU-based MC dose calculations. A key concept called phase-space-ring (PSR) was proposed. It contained a group of particles that are of the same type and close in energy and radial distance to the center of the phase-space plane. The model parameterized probability densities of particle location, direction and energy for each primary photon PSR, scattered photon PSR and electron PSR. For a primary photon PSRs, the particle direction is assumed to be from the beam spot. A finite spot size is modeled with a 2D Gaussian distribution. For a scattered photon PSR, multiple Gaussian components were used to model the particle direction. The direction distribution of an electron PSRs was also modeled as a 2D Gaussian distribution with a large standard deviation. We also developed a method to analyze a phase-space file and derive corresponding model parameters. To test the accuracy of our linac source model, dose distributions of different open fields in a water phantom were calculated using our source model and compared to those directly calculated using the reference phase-space file. The average distance-to-agreement (DTA) was within 1 mm for the depth dose in the build-up region and beam penumbra regions. The root-mean-square (RMS) dose difference was within 1.1% for dose profiles at inner and outer beam regions. The maximal relative difference of output factors was within 0.5%. Good agreements were also found in an IMRT prostate patient case and an IMRT head-and-neck case. These results demonstrated the efficacy of our source model in terms of accurately representing a reference phase-space file.

preprint2015arXiv

An Analytic Linear Accelerator Source Model for Monte Carlo dose calculations. II. Model Utilization in a GPU-based Monte Carlo Package and Automatic Source Commissioning

We recently built an analytical source model for GPU-based MC dose engine. In this paper, we present a sampling strategy to efficiently utilize this source model in GPU-based dose calculation. Our source model was based on a concept of phase-space-ring (PSR). This ring structure makes it effective to account for beam rotational symmetry, but not suitable for dose calculations due to rectangular jaw settings. Hence, we first convert PSR source model to its phase-space let (PSL) representation. Then in dose calculation, different types of sub-sources were separately sampled. Source sampling and particle transport were iterated. So that the particles being sampled and transported simultaneously are of same type and close in energy to alleviate GPU thread divergence. We also present an automatic commissioning approach to adjust the model for a good representation of a clinical linear accelerator . Weighting factors were introduced to adjust relative weights of PSRs, determined by solving a quadratic minimization problem with a non-negativity constraint. We tested the efficiency gain of our model over a previous source model using PSL files. The efficiency was improved by 1.70 ~ 4.41, due to the avoidance of long data reading and transferring. The commissioning problem can be solved in ~20 sec. Its efficacy was tested by comparing the doses computed using the commissioned model and the uncommissioned one, with measurements in different open fields in a water phantom under a clinical Varian Truebeam 6MV beam. For the depth dose curves, the average distance-to-agreement was improved from 0.04~0.28 cm to 0.04~0.12 cm for build-up region and the root-mean-square (RMS) dose difference after build-up region was reduced from 0.32%~0.67% to 0.21%~0.48%. For lateral dose profiles, RMS difference was reduced from 0.31%~2.0% to 0.06%~0.78% at inner beam and from 0.20%~1.25% to 0.10%~0.51% at outer beam.

preprint2013arXiv

GPU-based Monte Carlo radiotherapy dose calculation using phase-space sources

A novel phase-space source implementation has been designed for GPU-based Monte Carlo dose calculation engines. Due to the parallelized nature of GPU hardware, it is essential to simultaneously transport particles of the same type and similar energies but separated spatially to yield a high efficiency. We present three methods for phase-space implementation that have been integrated into the most recent version of the GPU-based Monte Carlo radiotherapy dose calculation package gDPM v3.0. The first method is to sequentially read particles from a patient-dependent phase-space and sort them on-the-fly based on particle type and energy. The second method supplements this with a simple secondary collimator model and fluence map implementation so that patient-independent phase-space sources can be used. Finally, as the third method (called the phase-space-let, or PSL, method) we introduce a novel strategy to pre-process patient-independent phase-spaces and bin particles by type, energy and position. Position bins located outside a rectangular region of interest enclosing the treatment field are ignored, substantially decreasing simulation time. The three methods were validated in absolute dose against BEAMnrc/DOSXYZnrc and compared using gamma-index tests (2%/2mm above the 10% isodose). It was found that the PSL method has the optimal balance between accuracy and efficiency and thus is used as the default method in gDPM v3.0. Using the PSL method, open fields of 4x4, 10x10 and 30x30 cm2 in water resulted in gamma passing rates of 99.96%, 99.92% and 98.66%, respectively. Relative output factors agreed within 1%. An IMRT patient plan using the PSL method resulted in a passing rate of 97%, and was calculated in 50 seconds using a single GPU compared to 8.4 hours (per CPU) for BEAMnrc/DOSXYZnrc.

preprint2013arXiv

The fixed-point iteration method for IMRT optimization with truncated dose deposition coefficient matrix

In the treatment plan optimization for intensity modulated radiation therapy (IMRT), dose-deposition coefficient (DDC) matrix is often pre-computed to parameterize the dose contribution to each voxel in the volume of interest from each beamlet of unit intensity. However, due to the limitation of computer memory and the requirement on computational efficiency, in practice matrix elements of small values are usually truncated, which inevitably compromises the quality of the resulting plan. A fixed-point iteration scheme has been applied in IMRT optimization to solve this problem, which has been reported to be effective and efficient based on the observations of the numerical experiments. In this paper, we aim to point out the mathematics behind this scheme and to answer the following three questions: 1) whether the fixed-point iteration algorithm converges or not? 2) when it converges, whether the fixed point solution is same as the original solution obtained with the complete DDC matrix? 3) if not the same, whether the fixed point solution is more accurate than the naive solution of the truncated problem obtained without the fixed-point iteration? To answer these questions, we first performed mathematical analysis and deductions using a simplified fluence map optimization (FMO) model. Then we conducted numerical experiments on a head-and-neck patient case using both the simplified and the original FMO model. Both our mathematical analysis and numerical experiments demonstrate that with proper DDC matrix truncation, the fixed-point iteration can converge. Even though the converged solution is not the one that we obtain with the complete DDC matrix, the fixed-point iteration scheme could significantly improve the plan accuracy compared with the solution to the truncated problem obtained without the fixed-point iteration.

preprint2012arXiv

Four-dimensional Cone Beam CT Reconstruction and Enhancement using a Temporal Non-Local Means Method

Four-dimensional Cone Beam Computed Tomography (4D-CBCT) has been developed to provide respiratory phase resolved volumetric imaging in image guided radiation therapy (IGRT). Inadequate number of projections in each phase bin results in low quality 4D-CBCT images with obvious streaking artifacts. In this work, we propose two novel 4D-CBCT algorithms: an iterative reconstruction algorithm and an enhancement algorithm, utilizing a temporal nonlocal means (TNLM) method. We define a TNLM energy term for a given set of 4D-CBCT images. Minimization of this term favors those 4D-CBCT images such that any anatomical features at one spatial point at one phase can be found in a nearby spatial point at neighboring phases. 4D-CBCT reconstruction is achieved by minimizing a total energy containing a data fidelity term and the TNLM energy term. As for the image enhancement, 4D-CBCT images generated by the FDK algorithm are enhanced by minimizing the TNLM function while keeping the enhanced images close to the FDK results. A forward-backward splitting algorithm and a Gauss-Jacobi iteration method are employed to solve the problems. The algorithms are implemented on GPU to achieve a high computational efficiency. The reconstruction algorithm and the enhancement algorithm generate visually similar 4D-CBCT images, both better than the FDK results. Quantitative evaluations indicate that, compared with the FDK results, our reconstruction method improves contrast-to-noise-ratio (CNR) by a factor of 2.56~3.13 and our enhancement method increases the CNR by 2.75~3.33 times. The enhancement method also removes over 80% of the streak artifacts from the FDK results. The total computation time is ~460 sec for the reconstruction algorithm and ~610 sec for the enhancement algorithm on an NVIDIA Tesla C1060 GPU card.

preprint2011arXiv

GPU-based Low Dose CT Reconstruction via Edge-preserving Total Variation Regularization

High radiation dose in CT scans increases a lifetime risk of cancer and has become a major clinical concern. Recently, iterative reconstruction algorithms with Total Variation (TV) regularization have been developed to reconstruct CT images from highly undersampled data acquired at low mAs levels in order to reduce the imaging dose. Nonetheless, TV regularization may lead to over-smoothed images and lost edge information. To solve this problem, in this work we develop an iterative CT reconstruction algorithm with edge-preserving TV regularization to reconstruct CT images from highly undersampled data obtained at low mAs levels. The CT image is reconstructed by minimizing an energy consisting of an edge-preserving TV norm and a data fidelity term posed by the x-ray projections. The edge-preserving TV term is proposed to preferentially perform smoothing only on non-edge part of the image in order to avoid over-smoothing, which is realized by introducing a penalty weight to the original total variation norm. Our iterative algorithm is implemented on GPU to improve its speed. We test our reconstruction algorithm on a digital NCAT phantom, a physical chest phantom, and a Catphan phantom. Reconstruction results from a conventional FBP algorithm and a TV regularization method without edge preserving penalty are also presented for comparison purpose. The experimental results illustrate that both TV-based algorithm and our edge-preserving TV algorithm outperform the conventional FBP algorithm in suppressing the streaking artifacts and image noise under the low dose context. Our edge-preserving algorithm is superior to the TV-based algorithm in that it can preserve more information of fine structures and therefore maintain acceptable spatial resolution.

preprint2011arXiv

GPU-based Low-dose 4DCT Reconstruction via Temporal Non-local Means

Four-dimensional computed tomography (4DCT) has been widely used in cancer radiotherapy for accurate target delineation and motion measurement for tumors in thorax and upper abdomen areas. However, 4DCT simulation is associated with much higher imaging dose than conventional CT simulation, which is a major concern in its clinical application. Conventionally, each phase of 4DCT is reconstructed independently using the filtered backprojection (FBP) algorithm. The basic idea of our new algorithm is that, by utilizing the common information among different phases, the input information required to reconstruct image of high quality, and thus the imaging dose, can be reduced. We proposed a temporal non-local means (TNLM) method to explore the inter-phase similarity. All phases of the 4DCT images are reconstructed simultaneously by minimizing a cost function consisting of a data fidelity term and a TNLM regularization term. We utilized a forward-backward splitting algorithm and a Gauss-Jacobi iteration method to efficiently solve the minimization problem. The algorithm was also implemented on graphics processing unit (GPU) to achieve a high computational speed. Our reconstruction algorithm has been tested on a digital NCAT thorax phantom in three low dose scenarios. Our new algorithm generates visually much better CT images containing less image noise and streaking artifacts compared with the standard FBP algorithm. Quantitative analysis shows that much higher contrast-to-noise ratio and signal-to-noise ratio can be obtained using our algorithm. The total reconstruction time for all 10 phases of a slice ranges from 90 to 140 seconds on an NVIDIA Tesla C1060 GPU card.

preprint2010arXiv

Cryogenic temperatures as a path towards high-Q terahertz metamaterials

Optical properties of planar thin film metamaterials were measured at room and liquid nitrogen temperatures using terahertz time-domain spectroscopy. The operation of metamaterials at cryogenic temperatures is anticipated to be a promising path towards low-loss metamaterials since nonradiative losses are strongly suppressed due to higher charge mobility. A 14% increase in the quality factor of the resonances was experimentally observed. It was limited by the high electron scattering rate due to defects in thin films. Supplementary simulations assuming metamaterials made of thick films reveal a temperature controlled behaviour and a 40% increase of the Q-factor at 10 K.

preprint2010arXiv

Superconductor terahertz metamaterial

We characterize the behaviour of split ring resonators made up of high-transition temperature YBCO superconductor using terahertz time domain spectroscopy. The superconductor metamaterial shows sharp change in the transmission spectrum at the fundamental inductive-capacitive resonance and the dipole resonance as the temperature dips below the transition temperature. Our results reveal that the high performance of such a metamaterial is limited by material imperfections and defects such as cracks, voids and secondary phases which play dominant role in partially impeding the flow of current causing dissipation of energy and electrical resistance to appear in the superconductor film.

preprint2010arXiv

Terahertz superconducting plasmonic hole array

We demonstrate thermally tunable superconductor hole array with active control over their resonant transmission induced by surface plasmon polaritons . The array was lithographically fabricated on high temperature YBCO superconductor and characterized by terahertz-time domain spectroscopy. We observe a clear transition from the virtual excitation of the surface plasmon mode to the real surface plasmon mode. The highly tunable superconducting plasmonic hole arrays may have promising applications in the design of low-loss, large dynamic range amplitude modulation, and surface plasmon based terahertz devices.

preprint2009arXiv

Random terahertz metamaterials

Using terahertz time domain spectroscopy we investigate the normal incidence transmission through periodically and randomly arranged planar split ring resonators (SRRs). Introduction of positional disorder in metamaterials has no effect on the quality factor of the fundamental Inductive-Capacitive (LC) resonance. The dipole resonances undergo broadening and shift in their resonance frequencies. The experiment reveals that the randomly distributed SRR structures interact incoherently at LC resonance but couple coherently at the higher frequency dipole resonance.