@inproceedings {731,
	title = {Application of the thermodynamic state index: Mapping the Unified Mechanics Theory onto the critical phenomena of the 2D Ising model},
	booktitle = {Proceedings of the 44th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2026},
	month = {17-20 Jun 2026},
	pages = {SPP-2026-1D-03},
	address = {Los Ba{\~n}os},
	abstract = {Current predictive models for material failure rely heavily on empirical curve-fitting with derivations that do not start from first-principles{\textquoteright} physics. This study proposes a physics-informed framework for reliability by validating the Unified Mechanics Theory (UMT) against the exact analytical solution of the 2D Ising Model. The Thermodynamic State Index (TSI), a normalised measure of entropic disorder, was mapped onto the Onsager solution for the Ising lattice. The result demonstrates that the phase transition, representing the system{\textquoteright}s "structural" collapse, occurs at a specific critical state index Φc ≈ 0.264 for the square lattice. At this threshold, the magnetic susceptibility diverges, confirming that the TSI accurately tracks the onset of system-wide instability. This work bridges the gap between statistical mechanics and continuum damage mechanics, offering a rigorous, non-empirical metric for quantifying the proximity of a system to catastrophic failure or critical phenomena.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2026-1D-03},
	author = {Kandhalvi M. Asaali and Cristine D. Villagonzalo}
}
@inproceedings {729,
	title = {Implementation of the measurement-based variational quantum eigensolver on the 1D anisotropic XXZ spin-chain},
	booktitle = {Proceedings of the 44th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2026},
	month = {17-20 Jun 2026},
	pages = {SPP-2026-2H-04},
	address = {Los Ba{\~n}os},
	abstract = {The measurement-based variational quantum eigensolver (MBVQE) was applied to the 1D anisotropic XXZ Heisenberg model. Using Graphix, a Python-based software library capable of simulating MBQC patterns, and a classical optimizer, the anisotropy parameter was varied from J<sub>z</sub> = -1.5 to J<sub>z</sub> = 1.5 in half-integer increments for N = 2, 3, 4, 5 qubits. Optimization techniques were also applied to the system to ensure proper flow structure. The results showed that the constructed MBVQE has strong approximations for cases N = 2 and N = 3 whilst the cases N = 4 and N = 5 had weaker correlations between the approximation and the exact ground state energies. The study also observed that the MBVQE is not necessarily more costly than its circuit-based counterpart if proper flow structure and optimization techniques are applied.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2026-2H-04},
	author = {Matt Laurence S Mendoza and Cristine D. Villagonzalo}
}
@inproceedings {732,
	title = {Method comparison: Time-independent and time-dependent coupling on qubitdyne protocol with higher photon states},
	booktitle = {Proceedings of the 44th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2026},
	month = {17-20 Jun 2026},
	pages = {SPP-2026-PB-21},
	address = {Los Ba{\~n}os},
	abstract = {The qubitdyne protocol allows for quadrature detection of bosonic modes in stationary systems via repeated weak interactions with an auxiliary qubit. This approach circumvents the loss of photons in the cavity. While the earlier implementation of qubitdyne protocol utilizes constant coupling, Strandberg, et al. [<em>PRL</em> <strong>133</strong>, 063601 (2024)], suggests time-dependent coupling can significantly enhance the measurement process. In this study, we quantitatively analyze the efficiency of the qubitdyne protocol with constant and linearly increasing time-dependent couplings. The simulations were carried out with Fock, cat, and coherent states with increasing values of photons in the cavity (α = 4, 6, 8). The accuracy of the reconstruction was quantified by calculating the Kolmogorov-Smirnov (KS) statistic with respect to the ideal homodyne distributions. The results reveal that the constant coupling results in large deviations and high KS values up to ≈0.9 for the reconstructed quadrature distributions attributed to inefficiency in depleting the cavity. However, the time-dependent coupling results in low KS values (< 0.1) for all the reconstructed distributions. This reduces measurements by up to six times.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2026-PB-21},
	author = {Dharlene Kristina M. Rosete and Cristine D. Villagonzalo}
}
@inproceedings {733,
	title = {Quantum simulation for the 1D heat equation with a central Gaussian barrier potential via Schr{\"o}dingerization},
	booktitle = {Proceedings of the 44th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2026},
	month = {17-20 Jun 2026},
	pages = {SPP-2026-3C-05},
	address = {Los Ba{\~n}os},
	abstract = {We present a quantum simulation of the 1D heat equation with a central Gaussian barrier using the Schr{\"o}dingerization method. The non-unitary nature of the heat equation is subverted through a warped phase transformation. This transforms the heat equation into a unitary operator that can be evolved in a quantum circuit. Quantum simulations were performed using Qiskit with nx = 6 position registers and np = 2, 4, 6 momentum registers. The accuracy of the quantum simulation was evaluated by getting the l{\texttwosuperior} error relative to the classical matrix exponentiation method. Results show that increasing the momentum registers np and reducing the timestep Δt both reduce the simulation error, correctly reproducing the M-shaped solution induced by the Gaussian barrier.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2026-3C-05},
	author = {Red Christian F. Hernandez and Cristine D. Villagonzalo}
}
@inproceedings {712,
	title = {Dynamics of an open quantum Otto engine with a working substance in a harmonic oscillator potential well},
	booktitle = {Proceedings of the 43rd Samahang Pisika ng Pilipinas Physics Conference},
	year = {2025},
	month = {25{\textendash}28 Jun 2025},
	pages = {SPP-2025-1D-06},
	address = {Quezon City},
	abstract = {This study investigates the performance of an open system quantum Otto heat engine, whose working medium consists of particles confined in a harmonic oscillator potential. The expansion and compression of this system is controlled by the frequency of the oscillator. Thermal energy is transferred from the heat baths connected to the isochoric branches. The dynamical trajectories were plotted for each thermodynamic step and the work done by the engine was calculated. From the ratio of the useful work and input heat, we obtain the efficiency when the working substance consists of either bosons or fermions. Here bosons are found to be more efficient as they allow multiple occupancy of states.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2025-1D-06},
	author = {Gene Michael P. Herman and Cristine D. Villagonzalo}
}
@inproceedings {713,
	title = {External electromagnetic effects on the spin transport of the open one-dimensional single-band Fermi-Hubbard lattice},
	booktitle = {Proceedings of the 43rd Samahang Pisika ng Pilipinas Physics Conference},
	year = {2025},
	month = {25{\textendash}28 Jun 2025},
	pages = {SPP-2025-2E-03},
	address = {Quezon City},
	abstract = {Experimentally, conductivities help probe phases of matter among strongly correlated systems. In this study, the effect of external electromagnetic field on the spin conductivity of an open one-dimensional single-band Fermi-Hubbard lattice is calculated with the assumption of linear response. It is shown that for repulsive two-particle interactions, the spin conductivity can be optimized at finite low-temperature when the lattice is subject to a strong magnetic field and a more positive electric field. On the other hand, external electromagnetic fields drive the spin conductivity to an asymptotic behavior when the two-particle interactions are attractive.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2025-2E-03},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@inproceedings {711,
	title = {One dimensional P{\"o}schl-Teller potential scattering as implemented in a quantum computing simulator},
	booktitle = {Proceedings of the 43rd Samahang Pisika ng Pilipinas Physics Conference},
	year = {2025},
	month = {25{\textendash}28 Jun 2025},
	pages = {SPP-2025-1D-04},
	address = {Quezon City},
	abstract = {The scattering of a quantum particle in one dimension from a P{\"o}schl-Teller potential is demonstrated using a quantum simulator, Qiskit, consisting of 10 qubits. Solving the Schr{\"o}dinger equation with this potential is implemented as the action on the qubits by Fourier transforms and diagonal operators represented as controlled-phase shift gates. Hence, a quantum algorithm is constructed that is capable of modulating the reflection and transmission of a quantum particle by varying the shape of the P{\"o}schl-Teller potential using a single parameter λ. The results obtained from Qiskit are verified using symbolic architecture simulations.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2025-1D-04},
	author = {Walter S. Go and Dharlene Kristina M. Rosete and Red Christian F. Hernandez and Cristine D. Villagonzalo}
}
@article {694,
	title = {Enhancing thermoelectric efficiency in twisted bilayer PdSe$_{2}$},
	journal = {ACS Appl. Energy Mater.},
	volume = {7},
	year = {2024},
	pages = { 7798},
	abstract = {PdSe$_{2}$ is a promising two-dimensional (2D) material that has garnered significant attention for thermoelectric energy conversion applications due to its exceptional electronic and thermal properties. This study investigates the effect of twisting on the electron and phonon transport, as well as the thermoelectric properties of monolayer{\textendash}twisted bilayer{\textendash}monolayer PdSe$_{2}$ devices using full self-consistent ab initio nonequilibrium Green{\textquoteright}s function methods combined with density functional theory. Our findings reveal that 90{\textdegree} twist angle significantly enhances the figure of merit (ZT), achieving values of 3.95 at 300 K and 6.25 at 600 K, an 81.19\% increase at room temperature compared to nontwisted configurations. This enhancement is primarily due to the tunability of the band gap, weakened interlayer coupling, and boundary scattering at the interface. Additionally, our analysis suggests that adjusting the interlayer contact length could further optimize ZT. Achieving ZT values exceeding 3 in twisted bilayer PdSe$_{2}$ emphasizes the potential of exploiting twist modulations in 2D materials, representing significant progress toward developing highly efficient nanodevices capable of converting waste heat into electrical energy.},
	doi = {10.1021/acsaem.4c01336},
	author = {Jariyanee Prasongkit and Sirichok Jungthawan and Cristine D. Villagonzalo}
}
@inproceedings {675,
	title = {Implementation of a quantum Otto thermal machine using an XX spin chain model},
	booktitle = {Proceedings of the 42nd Samahang Pisika ng Pilipinas Physics Conference},
	year = {2024},
	month = {3{\textendash}6 Jul 2024},
	pages = {SPP-2024-PC-12},
	address = {Batangas City},
	abstract = {Quantum thermal machine models are essential to understanding the thermodynamics of nanoscale systems. A quantum Otto thermal machine with a spin XX chain working substance was considered. The different working regimes was demonstrated for a chain of 10 spins by calculating the work and heat absorbed for one cycle. It was found that at high bath temperatures, the quantum Otto engine mainly operates as either accelerator or a heat engine, while at low bath temperatures it operates mainly as a heater. Additionally, sinusoidal fluctuations in the heat and work plots are more apparent at low bath temperatures. This is presumably attributed to the low dimension of the system, which dissappears in the thermodynamic limit. These findings are precursor to optimal device design using Otto engines based on the spin XX chain model.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2024-PC-12},
	author = {Geronimo Allan Jerome Acosta and Cristine D. Villagonzalo}
}
@article {700,
	title = {Optimal polarization and susceptibility profiles of a half-filled open magnetoelectric Fermi-Hubbard optical dimer},
	journal = {Phys. Lett. A},
	volume = {525},
	year = {2024},
	pages = {129898},
	abstract = {Magnetoelectric materials, such as thin films and magnetoelectric antennae, exhibit tunable magnetic polarization using external electric field and electric polarization using external magnetic field at room temperature by introducing strain, for instance. A more straightforward method to optimize magnetoelectric effect is by adjusting the applied electromagnetic fields in systems modeled by the Fermi-Hubbard dimer. In this study, direct and converse magnetoelectric effects are theoretically observed in an open Fermi-Hubbard dimer by tuning the external electromagnetic fields. Results show that an external magnetic field can enhance the electric polarization and susceptibility of the dimer particles without the necessity of mechanical strain. Conversely, gradients in an external electric field modulate the magnetic polarization and susceptibility. This striking behavior in the Fermi-Hubbard dimer can serve as a physical model for developing quantum control of complex materials such as in biosensors, biomedical processes, and energy harvesting.},
	doi = {10.1016/j.physleta.2024.129898},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@inproceedings {676,
	title = {Quantum circuit implementation of the first-order perturbation in transverse-field Ising systems},
	booktitle = {Proceedings of the 42nd Samahang Pisika ng Pilipinas Physics Conference},
	year = {2024},
	month = {3{\textendash}6 Jul 2024},
	pages = {SPP-2024-PB-23},
	address = {Batangas City},
	abstract = {The transverse-field Ising (TI) model is the standard prototype from which the properties of magnetic systems are often scrutinized. In addition, inhomogeneities in these systems can be treated as perturbations to the TI model. Incorporating perturbation escalates the difficulty of simulating such magnetic systems, specially if it is accompanied by an increase in system size or when higher dimensions or more complex lattice structures are considered. In this work, a quantum computational approached is demonstrated to carry out perturbation theory in a TI system consisting of four spins with nearest neighbor interactions along the x-axis with a transverse field along the z-axis. Considered separately as perturbations were the next-nearest neighbor interactions along the x-axis and the nearest neighbor interactions along the y-axis. Through a series of simple quantum gates implemented in Qiskit, IBM{\textquoteright}s quantum computing software platform, the first order correction to the ground state of the TI system was obtained in both cases. This work on quantum circuits provides a foundation upon which more complex magnetic systems can be simulated.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2024-PB-23},
	author = {Mark Gabriel Danganan and Cristine D. Villagonzalo}
}
@conference {699,
	title = {Thermodynamics performance optimization of a Lipkin-Meshkov-Glick heat engine model with anisotropy},
	year = {2024},
	month = {10{\textendash}16 Nov 2024},
	publisher = {Asian Network School and Workshop on Complex Condensed Matter Systems 2024},
	address = {Chulalongkorn University, Bangkok and Suranaree University of Technology, Nakhon Ratchasima},
	abstract = {Heat engine implementation in the atomic or nanoscale involves thermodynamic and quantum properties of the working substance. Interest in the optimization of these quantum effects is driven by the rise of electronic devices in the nanoscale and the search for devices that can power quantum technologies in the future. It has been shown that quantum heat engines outperform their classical counterparts, and several physically realizable engines have been constructed using quantized spin models. We investigate the performance of the quantum heat engine following a Carnot cycle with four strokes {\textendash} two reversible isothermal and isomagnetic steps, connected to constant temperature heat baths. The working medium of this quantum device is the Lipkin-Meshkov-Glick (LMG) model consists of two or three spins. In this work, anisotropy is introduced to the working medium. The energy eigenstates are calculated by solving Schr{\"o}dinger{\textquoteright}s equation and the thermodynamic performance of each process is obtained using the canonical ensemble. The effect of anisotropy in the operation of the engine is seen in the energy eigenvalues as work production is maximized for specific anisotropy and applied magnetic field strengths. Optimization of the applied magnetic field for a given anisotropy and the change in the adjacent energy level spacing determines the efficiency of a quantum Carnot heat engine.},
	url = {https://indico.cern.ch/event/1450400/},
	author = {Gene Michael P. Herman and Cristine D. Villagonzalo}
}
@conference {644,
	title = {Caloric and spin transport profile of Fermi-Hubbard systems},
	year = {2023},
	month = {6{\textendash}10 Nov 2023},
	pages = {P.16},
	publisher = {Asian Network School and Workshop on Complex Condensed Matter Systems 2023},
	address = {IOP VAST, Vietnam},
	abstract = {An open one-dimensional Fermi-Hubbard optical lattice exposed to external electromagnetic fields is investigated using the grand canonical ensemble formalism. As the system was reduced to a dimer, the eigenstates were determined using the exact diagonalization method while tuning the two-particle interactions of the lattice, a varying chemical potential, and applied external electromagnetic fields. Entropy and heat capacity gradients exist due to the interaction of the system with the external magnetic and electric fields, also known as magnetocaloric and electrocaloric effects, respectively. At particular temperatures, the isothermal entropic change was calculated to determine the necessary conditions for direct or inverse caloric effects to occur. The effects of more lattice sites and the elastic two-particle interaction can be incorporated through analytical methods. It was shown that at low temperatures, spins tend to scatter under the influence of repulsive two-particle interactions, while backflow occurs in the attractive regime. Furthermore, suitable temperature ranges for spin conductivity to occur can be found by adjusting the electromagnetic fields.},
	url = {https://iop.vast.vn/asian/2023/program.php?t=1\&n=16},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@inproceedings {624,
	title = {Effects of interparticle interaction on the spin transport in a half-filled open elastic Fermi-Hubbard optical lattice},
	booktitle = {Proceedings of the 41st Samahang Pisika ng Pilipinas Physics Conference},
	year = {2023},
	month = {19{\textendash}21 July 2023},
	pages = {SPP-2023-1F-05},
	address = {Siargao},
	abstract = {Spin is a particle property whose study of transport receives attention due to lossless transfer of information. In this study, spin transport is dependent on interparticle interactions in a one-dimensional half-filled open elastic Fermi-Hubbard optical lattice. Near absolute zero temperature temperature, spins flow in the repulsive regime and are not transported in the attractive regime. In addition, as the temperature approaches zero, the abrupt increase of the conductivity for the repulsive regime and the decrease to zero of the conductivity for the attractive regime take place at relatively higher temperatures.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2023-1F-05},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@inproceedings {626,
	title = {Election opinion dynamics through the 2D Ising model},
	booktitle = {Proceedings of the 41st Samahang Pisika ng Pilipinas Physics Conference},
	year = {2023},
	month = {19{\textendash}21 July 2023},
	pages = {SPP-2023-PB-45},
	address = {Siargao},
	abstract = {The Ising model and its variations have seen regular applications in describing the behavior of systems with many interacting bodies of varying complexity. We explored a zero-field Ising model that includes next nearest neighbor interactions with strength J$_{2}$ to investigate the opinion dynamics of a voting population in a two-party election. We constructed lattices filled with voters preferring to vote for either x candidate or x{\textquoteright} candidate in random configurations. We then allowed these sites to switch preferences using the Metropolis algorithm while varying the influence strength J$_{2}$, of the "clique" of each voter. Our work demonstrates that, when exposing the voting population to biased information/disinformation, increasing J$_{2}$ leads to a decrease in the critical bias of the voters where the overall "alignment" of their voting preferences drops to zero. We also investigated how this introduction of next nearest neighbor interactions is affected by changes in voting population size. Here, we found out that the rate of change at which the critical bias decreases is reduced as the population becomes larger due to the sheer amount of voters inhibiting long-range order among them.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2023-PB-45},
	author = {Jose Antonio R. Gamboa and Cristine D. Villagonzalo}
}
@conference {645,
	title = {Enhanced efficiency of a quantum Otto engine based on an XX spin chain near its dynamical phase transition},
	year = {2023},
	month = {6{\textendash}10 Nov 2023},
	pages = {P.26},
	publisher = {Asian Network School and Workshop on Complex Condensed Matter Systems 2023},
	address = {IOP VAST, Vietnam},
	abstract = {Implementation of heat engines in the nanoscale must consider quantum effects. In this study, a quantum Otto engine with an XX spin chain as a working substance was investigated near its dynamical phase transition for its efficiency enhancement. We used two isothermal and two isomagnetic branches for the thermodynamic cycle of the quantum Otto engine. The working substance was operated at the ground state and at a thermal state. The isothermal branches were implemented at zero temperature and inverse temperature βJ >= 1, while the isomagnetic branches at constant  |λ/J| < 2. The quantum heat engine was found to achieve a 100\% peak efficiency near the ground state curve in the chemical potential diagram and operates at constant 100\% efficiency for extended regions when the system undergoes dynamical phase transition.},
	url = {https://iop.vast.vn/asian/2023/program.php?t=1\&n=26},
	author = {Geronimo Allan Jerome Acosta and Cristine D. Villagonzalo}
}
@inproceedings {625,
	title = {Entropy and specific heat features of a four spin system at its energy level crossings},
	booktitle = {Proceedings of the 41st Samahang Pisika ng Pilipinas Physics Conference},
	year = {2023},
	month = {19{\textendash}21 July 2023},
	pages = {SPP-2023-PA-33},
	address = {Siargao},
	abstract = {A previous study showed that the efficiency of a quantum heat engine utilizing the Lipkin-Meshkov-Glick (LMG) spin system can reach the Carnot limit at low temperatures when the external field reaches the energy level cross points of the ground states. These points occur when different energy levels become degenerate. To delve more into the quantum behavior at the crossing points, this work considered the thermodynamic characteristics of an LMG model consisting of four spins. Using the eigenenergies for the N = 4 spin system with anisotropy, the entropy, and specific heat were numerically obtained in the canonical ensemble and studied as a function of its energy level crossings. It is found that lower values of energy level crossings correspond to more stable entropies and higher specific heat values of the system. These features make the LMG spin system suitable as a working substance in heat engine applications.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2023-PA-33},
	author = {Andrei Espencer M Ordonio and Cristine D. Villagonzalo}
}
@article {649,
	title = {Magnetocaloric and electrocaloric heat engine and refrigeration cycles in the open Fermi-Hubbard optical dimer with attractive interactions},
	journal = {Phys. Lett. A},
	volume = {481},
	year = {2023},
	pages = {129011},
	abstract = {Quantum heat engines and refrigerators can be realized in a Fermi-Hubbard optical dimer through immersion in a mean-field Fermi gas bath using the magnetocaloric and electrocaloric effects. It is demonstrated that in the presence of onsite and heat bath-induced attractive interactions, singlet-to-triplet transitions are achieved with a finite external magnetic field. Similarly, triplet-to-singlet transitions occur in the presence of an applied finite electric field. Results show that to lift the degeneracies, considerable electric potential magnitudes are needed. These magnitudes are greater than the magnetic potential strengths where the degeneracies occur. Furthermore, in the attractive regime, combinations of direct and inverse caloric effects can be combined to operate quantum heat engines and refrigerators, like what happens in the repulsive regime.},
	doi = {10.1016/j.physleta.2023.129011},
	url = {https://www.sciencedirect.com/science/article/abs/pii/S0375960123003912},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@inproceedings {623,
	title = {Quantum heat engine efficiency of a small spin system with anisotropy},
	booktitle = {Proceedings of the 41st Samahang Pisika ng Pilipinas Physics Conference},
	year = {2023},
	month = {19{\textendash}21 July 2023},
	pages = {SPP-2023-3H-05},
	address = {Siargao},
	abstract = {The rise of electronic devices lower than the nanoscale shifts the focus to quantum heat engines as they have improved features than their classical counterparts. In this work, the two-spin and three-spin systems using the Lipkin-Meshkov-Glick (LMG) model as the working substance of the heat engine are considered. The efficiency of the quantum heat engine comprised of two reversible isothermal processes and two isomagnetic processes with varying spin anisotropy is investigated. The configuration with the highest anisotropy was found to yield the highest efficiency for the two-spin system, and the most isotropic case was the most efficient for the three-spin system. This is attributed to the effect of anisotropy on the energy levels of the LMG spin system.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2023-3H-05},
	author = {Gene Michael P. Herman and Cristine D. Villagonzalo}
}
@article {577,
	title = {Caloric effects in an open Fermi{\textendash}Hubbard optical dimer due to onsite and heat bath-induced two-particle interactions},
	journal = {J. Magn. Magn. Mater.},
	volume = {564},
	year = {2022},
	pages = {170094},
	abstract = {The immersion of the Fermi{\textendash}Hubbard optical dimer in a mean field Fermi gas, acting as a heat bath, modifies the behavior of the statistical mechanical properties of the particles trapped in it. This is done by varying the two-particle interactions induced by the heat bath and the onsite two-particle interaction <em>U</em>. It is shown that in the presence of the heat bath-induced interaction, the partition function <em>Z</em> drops to a minimum finite value in the repulsive regime of <em>U</em> at the normalized temperature region <em>kT</em> ≲ <em>t/</em>2, where <em>t</em> is the tunneling parameter. In this temperature region, the convergence points of the respective isolines of the entropy <em>S</em> and the heat capacity <em>Cᵥ</em> move towards each other. Such shift in their features finetune <em>S</em> and <em>Cᵥ</em> relative to their behavior in the absence of two-particle interactions with the bath.},
	doi = {10.1016/j.jmmm.2022.170094},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@inproceedings {583,
	title = {Chemical potential-induced effects on caloric properties of a Fermi-Hubbard dimer},
	booktitle = {Proceedings of the 40th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2022},
	month = {19{\textendash}21 Oct 2022},
	pages = {SPP-2022-1G-04},
	address = {Legazpi City},
	abstract = {Chemical potential is one of the tunable properties of an optical lattice system when its temperature is varied. In this study, the effects of the chemical potential on the entropy and the heat capacity of a Fermi-Hubbard dimer immersed in a mean-field Fermi gas bath are calculated. Repulsive two-particle interactions are found to require lesser chemical potentials to observe a rapid increase in the entropy at low temperatures. Such phenomena signify quantum level crossings. Furthermore, half-filling of the optical dimer is unnecessary for such phenomena to occur.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2022-1G-04},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@inproceedings {581,
	title = {Crossing and anti-crossing of the Lipkin-Meshkov-Glick model for N = 4 spin system},
	booktitle = {Proceedings of the 40th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2022},
	month = {19{\textendash}21 Oct 2022},
	pages = {SPP-2022-2G-06},
	address = {Legazpi City},
	abstract = {Energy level crossings are degeneracies of two or more energy levels which occur when two levels coincide due to the tuning of an external parameter. In this work, we observe the "motion" of the energy levels in terms of the applied external magnetic field and anisotropy. We demonstrate how the eigenenergies of the anisotropic Lipkin-Meshkov-Glick (LMG) model for the N = 4 spin system shift as a function of the aforementioned parameters. The resulting plots show multiple level crossings as the system becomes isotropic where the LMG model has degenerate ground states along the critical lines determined by the anisotropy parameter. This work has also verified that for the isotropic case, energy level crossings occur exactly at the critical points.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2022-2G-06},
	author = {Andrei Espencer M Ordonio and Klarence Anne M. Feri and Cristine D. Villagonzalo}
}
@inproceedings {585,
	title = {Exploring quantum phenomena of various fermionic quantum models through their thermodynamic properties},
	booktitle = {Proceedings of the 40th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2022},
	month = {19{\textendash}21 Oct 2022},
	pages = {SPP-2022-INV-1A-01},
	address = {Legazpi City},
	abstract = {Thermodynamic properties are usually associated with the macroscopic or bulk properties of a system. In the nanoscale of condensed fermionic systems, their thermodynamics properties such as the thermoelectric coefficient, thermal conductivity or heat capacity, among others, can be used to probe quantum phenomena. Different theoretical models for fermions that have a distinct feature, such as (i) disorder, (ii) two-particle interaction, or (iii) spin-orbit coupling, are explored at low temperature to show how a temperature gradient affects their thermodynamics properties. In turn, these offer insight into the quantum nature of systems described by these fermionic models.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2022-INV-1A-01},
	author = {Cristine Villagonzalo}
}
@inproceedings {582,
	title = {Quantum Carnot heat engine with an effective external electric field perturbation},
	booktitle = {Proceedings of the 40th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2022},
	month = {19{\textendash}21 Oct 2022},
	pages = {SPP-2022-1D-03},
	address = {Legazpi City},
	abstract = {The classical approach to heat engines is not enough to study devices in the atomic scale as it utilizes spin-systems as working substances whose behavior is dictated by quantum mechanics. We constructed a quantum heat engine that follows the Carnot cycle with two reversible isothermal and isomagnetic processes in contact with heat baths of two different temperatures, based on the Lipkin-Meshkov-Glick (LMG) model. The engine runs between two magnetic fields with one kept constant and the other magnetic field varied. A perturbation in the form of an effective external electric field in the x-direction, δJx, was introduced to the system where δ << 1. Up to the second-order energy correction was determined using the time-independent perturbation theory and the unperturbed LMG eigenstates as basis. The entropy, internal energy, and the work done in each cycle were determined and used to obtain the efficiency of the system with spins N = 2. For a fixed set of values for the high-temperature heat bath and the low-temperature heat bath, this work demonstrates that the system{\textquoteright}s efficiency is dependent on δ and the varying magnetic field. The system{\textquoteright}s maximum efficiency was found to approach the value of the Carnot limit.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2022-1D-03},
	author = {Gene Michael P. Herman and Cristine D. Villagonzalo}
}
@article {574,
	title = {Refrigeration using magnetocaloric and electrocaloric effects in a Fermi-Hubbard optical dimer exposed to a heat bath},
	journal = {Physica A},
	volume = {600},
	year = {2022},
	pages = {127540},
	abstract = {We consider a Fermi{\textendash}Hubbard dimer exposed to external electromagnetic fields and a Fermi gas bath under the assumptions of a repulsive two-particle interaction and the half-filling condition. Using the exact diagonalization method, it was found that the entropy and heat capacity rapidly rise within distinct localized regions of the magnetic and electric fields at low temperatures. The behavior of the magnetic (electric) Gr{\"u}neisen ratio denotes that the direction of the external field determines the magnetization or demagnetization (polarization or depolarization) due to the magnetocaloric (electrocaloric) effect among the lattice particles at a given temperature. At the same time, the isothermal entropic change as a function of the temperature indicates that refrigeration is optimized when the electric field along the line joining the two sites of the optical lattice is assisted by a transverse magnetic field.},
	doi = {10.1016/j.physa.2022.127540},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@conference {593,
	title = {Spin transport in a Fermi-Hubbard optical lattice interacting with a Fermi gas bath},
	year = {2022},
	month = {18{\textendash}24 Aug 2022},
	pages = {P23-SF4-07},
	publisher = {29th International Conference on Low Temperature Physics},
	address = {Presented online},
	url = {https://www.lt29.jp/},
	author = {Vladimir P. Villegas and Cristine D. Villagonzalo}
}
@inproceedings {villegas2021,
	title = {Canonical ensemble profile of a Fermi-Hubbard model treated as a double qubit chain},
	booktitle = {Proceedings of the 39th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2021},
	month = {20{\textendash}22 Oct 2021},
	pages = {SPP-2021-2B-03},
	address = {Presented online},
	abstract = {Through the Jordan-Wigner transformation, a Fermi-Hubbard (FH) optical lattice can be mapped onto a double qubit chain system (QS). In this study, the partition function and the Helmholtz free energy were calculated for a FH dimer and a double QS double well. The differences of the values of the two properties as functions of temperature were taken and are shown as saturating to a constant as the temperature increases. This implies that electrons can be emulated as qubits.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2021-2B-03},
	author = {Vladimir P. Villegas and Cristine Villagonzalo}
}
@article {566,
	title = {Optimal fixed-point quantum search in an interacting Ising spin system},
	journal = {Quantum Inf. Process.},
	volume = {20},
	year = {2021},
	pages = {90},
	abstract = {In Grover{\textquoteright}s search algorithm, a priori knowledge of the number of target states is needed to effectively find a solution. This is due to the inherent oscillatory nature of unitary gates in the algorithm. A fixed-point quantum search is introduced in [Phys. Rev. Lett. 113, 210501 (2014)] to mitigate this oscillation even without knowing the size of the target space. This is done by modifying the phase inversion and oracle of the original Grover{\textquoteright}s algorithm so that the phase can assume intermediate values in [-π,π]. Naturally, qubits in a quantum computer are interacting among themselves and this might introduce an error in the phase assignment. In this work, we used the Ising spin chain to simulate the interactions among the qubits and demonstrate how to implement the fixed-point quantum search. By inspecting the state vectors after some finite number of iterations, we found that the error is primarily due to the phase difference of the ideal target state with that of the final state. Further investigation shows that high fidelity between the ideal and real evolution can be achieved by using a low Rabi frequency.},
	doi = {10.1007/s11128-021-03023-1},
	author = {Neris I. Sombillo and Ronald S. Banzon and Cristine Villagonzalo}
}
@inproceedings {feri2021,
	title = {Specific heat behavior in finite-sized Lipkin-Meshkov-Glick spin system},
	booktitle = {Proceedings of the 39th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2021},
	month = {20{\textendash}22 Oct 2021},
	pages = {SPP-2021-1F-05},
	address = {Presented online},
	abstract = {A quantum heat engine whose working substance is described by the Lipkin-Meshkov-Glick (LMG) model has been shown to improve its efficiency when the LMG system is near its quantum phase transition. It is of interest to know how the specific heat of such an engine will behave under similar conditions. Therefore, in this work, we construct the specific heat curves of the LMG spin system consisting of spin-1/2 particles in a plane of sizes N = 2, 4, 6, 8. Effects of the field strength, temperature, and system size on the specific heat are analyzed. Salient features of the specific heat curve as a function of field strength are observed in relation to the ground state degeneracy that occurs in key field strength values.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2021-1F-05},
	author = {Klarence Anne M. Feri and Cristine Villagonzalo}
}
@inproceedings {552,
	title = {Characterization of phonon density of states of a graphene junction beyond nearest neighbor interactions},
	booktitle = {Proceedings of the 38th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2020},
	month = {19{\textendash}23 Oct 2020},
	pages = {SPP-2020-4G-05},
	address = {Presented online},
	abstract = {Graphene has exceptional thermal properties that has presented a wide range of applications. Ballistic thermal transport in graphene due to phonons was made possible by placing a zigzag graphene junction in the middle of a heat sink and a heat source- which induced a temperature gradient that caused heat to flow through the junction. Following the harmonic oscillator model for the graphene junction, interlayer and intralayer interactions were derived beyond nearest neighbors. This paper aims to characterize the zigzag graphene junction{\textquoteright}s phonon density of states and examine the effects of adding next nearest neighbors and varying the bond stretching force constant, ks, on the next nearest neighbor interactions. We observed a wider distribution of energies when adding next nearest neighbor interactions at equal ks. Varying ks showed similar results as this indicates the strength of interaction between the next nearest neighboring atoms.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2020-4G-05},
	author = {Krishna Lyn Delima and Cristine Villagonzalo}
}
@inproceedings {551,
	title = {Drude model electron motion in 2D space with applied electric fields},
	booktitle = {Proceedings of the 38th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2020},
	month = {19{\textendash}23 Oct 2020},
	pages = {SPP-2020-4G-04},
	address = {Presented online},
	abstract = {The key to the research on emergent materials is the investigation of their response to modulated fields. Here in this work, the differential equation for electron motion is solved and plotted using Fortran and gnuplot to analyze the charge-carrier dynamical behavior in two-dimensional (2D) materials through the Drude model. Momentum decay of an electron in motion inside a material is compared with three other cases of applied electric fields. Properties such as speed and current density of the electron are also determined to show how this old model still functions in today{\textquoteright}s science.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2020-4G-04},
	author = {Tim Jeff M. Rodriguez and Cristine Villagonzalo}
}
@inproceedings {549,
	title = {Free energy in a spin transport among particles in a Fermi-Hubbard optical lattice},
	booktitle = {Proceedings of the 38th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2020},
	month = {19{\textendash}23 Oct 2020},
	pages = {SPP-2020-5C-03},
	address = {Presented online},
	abstract = {The Fermi-Hubbard model investigates the dynamics of fermions confined in a one-dimensional (1D) optical lattice. In turn, the lattice is confined to a heat bath consisting of a fermionic quantum gas. Here, we use the path integral method to investigate the effects of the coupling at zero temperature and determine the five processes among bath particles which facilitate spin transport among optical lattice particles. Furthermore, we compute for the free energy per lattice site during such process.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2020-5C-03},
	author = {Vladimir P. Villegas and Cristine Villagonzalo}
}
@inproceedings {550,
	title = {Temperature dependence of three-dimensional thermoelectric properties of a free electron gas-like material},
	booktitle = {Proceedings of the 38th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2020},
	month = {19{\textendash}23 Oct 2020},
	pages = {SPP-2020-4G-01},
	address = {Presented online},
	abstract = {<p>Regeneration of electricity from waste heat can be done through the use of thermoelectric media. In this work, how the bulk thermoelectric properties vary with temperature were simulated using the Fermi energy for Bi<sub>2</sub>Te<sub>3</sub>{\textemdash}a promising material for thermoelectric generators. The numerical calculations were based on the equations derived by J. Muldera in 2006 for a free electron gas in three dimensions. The temperature dependence of the thermoelectric properties, such as the thermoelectric power, thermal conductivity and the figure of merit were determined by keeping the number density of charge carriers constant at their Fermi level value. The results show that such a material with a Bi<sub>2</sub>Te<sub>3</sub> Fermi level has thermoelectric properties suited for room temperature applications.</p>},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2020-4G-01},
	author = {Jairus Carlo F. Calpatura and Cristine Villagonzalo}
}
@inproceedings {acostaSPP2019,
	title = {Current amplification and rectification in a chain of four spins individually coupled to different heat baths},
	booktitle = {Proceedings of the 37th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2019},
	month = {29 May 2019},
	pages = {SPP-2019-PB-33},
	address = {Tagbilaran City, Philippines},
	abstract = {In this work, we considered a linear chain consisting of four spins individually coupled to different heat baths. We calculated numerically the heat currents by solving the transition rate equations as previously discussed in Ref. [1]. The temperature of the heat bath of the first spin is held constant. The temperatures of the second and third heat baths are set to be the average of the heat baths at the edges while the temperature of the fourth heat bath at the edge is allowed to vary from zero to twice the temperature of the first heat bath. We demonstrate that for a chain of spins conjugate to this heat bath temperature configuration, current amplification and rectification are exhibited.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2019-PB-33},
	author = {Geronimo Allan Jerome Acosta and Cristine Villagonzalo}
}
@inproceedings {mactalSPP2019,
	title = {Effects of spatial disorder in modelling the anomalous diffusion of tracer particles on a graphene lattice},
	booktitle = {Proceedings of the 37th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2019},
	month = {29 May 2019},
	pages = {SPP-2019-PB-18},
	address = {Tagbilaran City, Philippines},
	abstract = {Modern microscopic techniques such as high resolution transmission electron microscopy have been developed by experimentalists in order to track the anomalous diffusion of tracer particles. By using the site percolation model, this paper investigated the effects of spatial disorder on a graphene lattice and how parameters such as the percolation threshold contributes to a significant deviation from the linear diffusion equation. In order to obtain the numerical percolation threshold, a cluster-identifying algorithm was developed. This value was found to be 0.7139 which has a theoretical deviation of 2.542\%. After establishing this, simulations of the particle diffusing on the percolating system with p=0.7139 were done and its mean squared displacement (MSD) at each time step t was extracted. After analysis of the log-log plot of the MSD vs. number of steps taken, it was found that the particle exhibited subdiffusive behavior for the different cases tested because it violated Einstein{\textquoteright}s third condition for linear diffusion.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2019-PB-18},
	author = {Gerard Antonio L Mactal and Cristine Villagonzalo}
}
@inproceedings {delimaSPP2019,
	title = {Interlayer and intralayer interactions in a graphene junction beyond nearest neighbors},
	booktitle = {Proceedings of the 37th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2019},
	month = {29 May 2019},
	pages = {SPP-2019-PB-34},
	address = {Tagbilaran City, Philippines},
	abstract = {The potential applications of the ballistic property of the thermal transport due to phonons in a zigzag graphene junction has prompted scientists to study the material{\textquoteright}s thermal properties. A zigzag graphene junction placed between a heat source and a heat sink has an induced temperature gradient wherein ballistic thermal transport becomes possible. Thermal transport studies for this system are usually characterized at the molecular level by nearest neighbor interactions between the carbon atoms within the junction. This paper aims to derive interlayer and intralayer interactions beyond nearest neighbors using a harmonic interaction potential and highlight their effects. We successfully determined these interactions by taking the Hamiltonian interaction matrix and illustrated a systematic way to add more interactions between graphene interlayers.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2019-PB-34},
	author = {Krishna Lyn Delima and Cristine Villagonzalo}
}
@inproceedings {naragSPP2019,
	title = {Investigation on the effect of tensile and compressive strain on the electronic structure of silicene using the density functional theory},
	booktitle = {Proceedings of the 37th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2019},
	month = {29 May 2019},
	pages = {SPP-2019-PA-16},
	address = {Tagbilaran City, Philippines},
	abstract = {We investigate the effect of tensile and compressive strain on the electronic structure of silicene using the density functional theory as implemented in the Quantum Espresso software. It was found that a tensile strain of about 7\% induces a semimetal-metal transition. Such a transition also occurs for 5\% applied compressive strain. Furthermore, an approximately linear relationship between an increasing strain and conduction (valence) band-Fermi energy difference was obtained for both types of strain. Engineering the band gap can widen the potential applications for mechanical modulation of the electronic properties of silicene.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2019-PA-16},
	author = {Dean Von Johari C. Narag and Cristine Villagonzalo}
}
@inproceedings {villanoSPP2019,
	title = {Monte Carlo approach to one-dimensional cell proliferation within a tissue-cell population},
	booktitle = {Proceedings of the 37th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2019},
	month = {29 May 2019},
	pages = {SPP-2019-PB-35},
	address = {Tagbilaran City, Philippines},
	abstract = {The interplay of biological, chemical, and physical factors of cancer cells give rise to its complex behavior. This complexity is one of the main reasons why the cure for cancer remains elusive. Finding a way to predict and understand this stochastic behavior would enable us to develop new and effective treatments that would hopefully lead to a cure for cancer. With this in mind, we constructed an algorithm for a stochastic one-dimensional cell proliferation through agent-based modeling and Monte Carlo methods. The proliferation of cells in a one-dimensional axis is demonstrated for different growth rates and cell-cell attraction strengths.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2019-PB-35},
	author = {Celso T. Villano Jr and Cristine Villagonzalo}
}
@inproceedings {coSPP2019,
	title = {Simulation of continuous-time random walks with anharmonicity},
	booktitle = {Proceedings of the 37th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2019},
	month = {29 May 2019},
	pages = {SPP-2019-PB-48},
	address = {Tagbilaran City, Philippines},
	abstract = {Anomalous diffusion processes are best described and analyzed by a continuous-time random walk (CTRW) formalism. The latter is a generalization of a transport model using random walk processes that are not limited to Markovian systems. Using a robust numerical generation of continuous trajectories based on Langevin equations, sample paths from a CTRW system encountering a nonlinear net force were generated and analyzed. The resultant force comprised of a linear repelling term and an anharmonic term, that is, f(x) = -k$_{1}$ x + k$_{2}$ x{\textthreesuperior}. The effect of anharmonicity through varying the k$_{2}$ values on the trajectories was manifested on the waiting time. This is the time at which the particle is relatively stationary. At a fixed stability index α, the results show that generally, waiting times would increase as k$_{2}$ values increased and that there is an optimum value of k$_{2}$ where the waiting time is longest. These results are compared to the case when only the linear force is present.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2019-PB-48},
	author = {Co, Richmond and Cristine Villagonzalo}
}
@conference {493,
	title = {Effect of error tolerance probability to the fidelity of a 4-qubit quantum search simulation},
	year = {2018},
	month = {5{\textendash}9 Nov 2018},
	pages = {P.16},
	publisher = {2018 ICTP Asian Network School and Workshop on Complex Condensed Matter Systems},
	address = {Suranaree University of Technology, Thailand},
	abstract = {We simulate the optimized fixed-point quantum search (OFPQS) algorithm with first- and second-nearest neighbor interactions in a 4-qubit Ising spin system using a nonrefocusing technique. In order to perform a single qubit gate, we utilize the coupling interactions to control the neighboring spins. The phase-marking scheme of the OFPQS is implemented using a single rotation about the transverse axis instead of the usual rotation about the field axis. Furthermore, the state of the ancilla is initialized in a superposition of the basis states to avoid the generation of unwanted states. This procedure gives us a probability of finding the target state within the prescribed bound of the OFPQS. Based on this scheme, we investigate the effect of error tolerance probability of the OFPQS to the fidelity of the nonrefocusing implementation.},
	url = {https://indico.cern.ch/event/745791/contributions/3202508/},
	author = {Neris Sombillo and Ronald Banzon and Cristine Villagonzalo}
}
@inproceedings {459,
	title = {Effect of varying the atomic bond stiffness in the heat baths to the thermal conductance of a graphene junction},
	booktitle = {Proceedings of the 36th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2018},
	month = {6{\textendash}9 June 2018},
	pages = {SPP-2018-3B-04},
	address = {Puerto Princesa City, Philippines},
	abstract = {We simulate ballistic thermal transport due to phonons in a graphene junction by letting it interact with heat baths located at its edges and thereby exposing it to a temperature gradient. The bond stretching force constant for the left (L, heat source), central (C, graphene junction) and right (R, heat sink) regions should ideally be the same since they are all made out of graphene with interactions described by the same interatomic potential. But in an experiment, the graphene sheets acting as heat baths are mounted on silicon oxide while the graphene junction in the central region is suspended with a trench underneath it to reduce its coupling to the silicon substrate. In such a case, the force constants may vary from one region to the next. In this paper, we investigate the effect of varying the atomic bond stiffness in the heat baths to the calculated thermal conductance of a graphene junction using the nonequilibrium Green{\textquoteright}s function method. We find that varying the force constants do not have a significant effect as the junction becomes larger and approaches the bulk behavior. At small junctions, the finite-size effect becomes relevant.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2018-3B-04},
	author = {Nica Jane B. Ferrer and Cristine Villagonzalo}
}
@inproceedings {466,
	title = {An effective field analysis of magnetic orbit-to-orbit interactions in a two-dimensional electron gas},
	booktitle = {Proceedings of the 36th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2018},
	month = {6{\textendash}9 June 2018},
	pages = {SPP-2018-PB-14},
	address = {Puerto Princesa City, Philippines},
	abstract = {An effective field Schr{\"o}dinger equation accounting for the effect of the magnetic orbit-to-orbit interaction between electrons in a two-dimensional electron gas (2DEG) is considered. This is formulated out from the effective single-particle equation using the Hartree approximation. As a specific example, a 2DEG in a Hall bar with rectangular geometry is considered, where such equation is reduced to a one-dimensional Schr{\"o}dinger equation. The effect of constant and inverse position-dependent magnetic orbit-to-orbit interaction potentials to the energy levels of the 2DEG is calculated perturbatively.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2018-PB-14},
	author = {Daniel F. Marquez and Cristine Villagonzalo}
}
@inproceedings {468,
	title = {Harmonic disturbance on the topological phases of a Chern insulator},
	booktitle = {Proceedings of the 36th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2018},
	month = {6{\textendash}9 June 2018},
	pages = {SPP-2018-PB-24},
	address = {Puerto Princesa City, Philippines},
	abstract = {We demonstrate the effect of a harmonic oscillator perturbing potential on the topological phases of a Chern insulator formed from a square lattice with two quantum states per site. For a region with a non-trivial topological phase characterized by a non-zero Chern number, the perturbation introduced an energy gap closing and divided the non-trivial region into three, which possessed a topologically trivial region and two non-trivial regions. The perturbation induced singularities on the Berry curvature which led to the division of the original topologically non-trivial region.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2018-PB-24},
	author = {Dean Von Johari C. Narag and Cristine Villagonzalo}
}
@inproceedings {469,
	title = {Impact of the absorption and reflection of signals on the radiation-induced bystander effect},
	booktitle = {Proceedings of the 36th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2018},
	month = {6{\textendash}9 June 2018},
	pages = {SPP-2018-PC-22},
	address = {Puerto Princesa City, Philippines},
	abstract = {The radiation-induced bystander effect (RIBE) implies radiation risks to cells or tissues not directly exposed to radiation. Effects were found to be strong in the low-dose range which humans are exposed to in the environment such as in medical activities and unpredicted accidents. One possible mechanism of RIBE is through signal transmission between irradiated cells and non-irradiated cells. This work modeled signal transmission through a 2D random walk (RW) of signals in a cell-culture medium containing 103 non-irradiated cells. The effect of boundary conditions (reflection vs. absorption) were observed. Results for the reflection case was better suited for experimental data than the absorption case. It showed independence of RIBE with distance from irradiated source.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2018-PC-22},
	author = {Zheina J. Ottong and Cristine Villagonzalo}
}
@inproceedings {463,
	title = {Investigating the diffusion of an atom on the graphene lattice},
	booktitle = {Proceedings of the 36th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2018},
	month = {6{\textendash}9 June 2018},
	pages = {SPP-2018-PB-04},
	address = {Puerto Princesa City, Philippines},
	abstract = {Atomic-scale surface diffusion in graphene is of recent interest to determine the viability of single-atom catalysts in functionalized graphene materials. To model such behavior, we study the random walk behavior of an atom in graphene using mean square displacement (MSD). Without the presence of a driving force, the atom performing a random walk on the graphene lattice would expect to follow linear Brownian motion. We confirm this linear behavior by measuring its mean square displacement (MSD). As the atom was given a driving force, the relationship between MSD and time steps begin to follow the properties of anomalous diffusion. The highest values of a was obtained when the atom was restricted to the first six nearest neighbors while the highest diffusion coeffcient Da was obtained when the atom is free to move to all lattice sites. The power-law relationship between MSD and time step can be controlled and serve as a precursor in understanding how the metal atom diffuses on graphene-oxide support.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2018-PB-04},
	author = {Gerard Antonio L Mactal and Cristine Villagonzalo}
}
@conference {482,
	title = {Magnetic orbit to orbit interaction in a two-dimensional electron gas},
	year = {2018},
	month = {18{\textendash}19 June 2018},
	pages = {P130},
	publisher = {ASTHRDP Graduate Scholars{\textquoteright} Conference},
	address = {PICC, Pasay City},
	abstract = {The orbit to orbit magnetic interaction between electrons in a two-dimensional electron gas (2DEG) is considered. The orbit to orbit interaction is realized by constructing a Lorentz force law of a system of electrons accounting for both Coulomb electrostatic and Biot-Savart magnetostatic interactions, and external fields. From this, a generalized potential energy and Lagrangian is constructed. Using the calculated mechanical momentum, the effective Hamiltonian is expressed from the energy functon where a sufficient order of inverse speed of light is only considered. Using the Hartree approximation, the solution if the many-body Schr{\"o}dinger equation with this Hamiltonian is approximated with an effective single-particle Schr{\"o}dinger equation. From this the effective mean field Schr{\"o}dinger equation is calculated. We look into the effects of the magnetic orbit-to-orbit interactiong through several forms of effective field potential energies.},
	author = {Daniel F. Marquez and Cristine Villagonzalo}
}
@inproceedings {472,
	title = {The next-next-nearest neighbor interaction of two Fibonacci anyons},
	booktitle = {Proceedings of the 36th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2018},
	month = {6{\textendash}9 June 2018},
	pages = {SPP-2018-3B-03},
	address = {Puerto Princesa City, Philippines},
	abstract = {The interaction between the next-next-nearest neighboring Fibonacci anyons was observed. The union requires the greatest proximity between the particles, hence the need for two braiding processes. The outcome of the fusion is projected onto the trivial channel, to signify a preference to the ferromagnetic state (i.e., having a negative coupling strength constant). The next-next-nearest neighbor interaction is shown to be more favorable than the next-nearest neighbor interaction.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2018-3B-03},
	author = {Celso T. Villano Jr and Cristine Villagonzalo}
}
@conference {494,
	title = {Numerical method for solving the semiclassical Boltzmann equation on multilayer heterogeneous structures},
	year = {2018},
	month = {5{\textendash}9 Nov 2018},
	pages = {P.27},
	publisher = {2018 ICTP Asian Network School and Workshop on Complex Condensed Matter Systems},
	address = {Suranaree University of Technology, Thailand},
	abstract = {We adapt the method developed by D.R. Penn and M. D. Stiles to implement a numerical solution to the semi-classical Boltzmann Equation on multilayer heterogeneous systems. This involves discretization of the linearized spin-independent Boltzmann Equation where the solution for each layer, considered as a bulk material, is obtained. The solutions for the different layers are then joined using boundary conditions at the interfaces to get the result for the whole multilayer structure. This method does not exploit the relaxation time approximation thus produces wider applicability especially to materials that cannot be considered to have a free-electron structure. The resulting distribution function is then used to calculate for the current density for a multi-layered material under a constant external electric field.},
	author = {April Cortez and Cristine Villagonzalo}
}
@inproceedings {470,
	title = {Radioactive diffusion without complete signal annihilation},
	booktitle = {Proceedings of the 36th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2018},
	month = {6{\textendash}9 June 2018},
	pages = {SPP-2018-PC-19},
	address = {Puerto Princesa City, Philippines},
	abstract = {The radiation-induced bystander effect (RIBE) is a non-contact phenomenon in which cells undergo genetic mutations despite the lack of a direct radiation traversal. In traditional models of the effect, the radiation-carrying signals travel through a medium via random walks and are annihilated upon causing successful mutations. In this paper, we analyze the effect of removing immediate signal annihilation upon successful contact. We find that the removal of immediate signal annihilation is significant only when the initial dosage, the mean number of re-emitted signals, or the probability of mutations are high. Future researchers may take the result into consideration when building models of RIBE.},
	url = {https://paperview.spp-online.org/proceedings/article/view/SPP-2018-PC-19},
	author = {Eduard Renzo Santos and Cristine Villagonzalo}
}
@inproceedings {sombillo,
	title = {Amplitude amplification of the marked states in quantum search simulation},
	booktitle = {Proceedings of the 35th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2017},
	month = {7{\textendash}10 June 2017},
	pages = {SPP-2017-PB-17},
	address = {Cebu City, Philippines},
	abstract = {Yoder et al. [1] proposed a quantum search algorithm that prevents the system to move away from the marked states when the number of target states M is unknown. This method uses a recursive approach that makes the target state act as a fixed-point. We simulate the algorithm in an Ising spin chain with first- and second-nearest neighbor interaction. Our method involves a selective phase-shift rotation per iteration to obtain the target state. We find the probability of success for the fixed-point implementation to fluctuate only about a limited range of values within an indicated tolerance as compared to the large oscillations of the same for Grover{\textquoteright}s non-fixed-point algorithm.},
	url = {http://paperview.spp-online.org/proceedings/article/view/142},
	author = {Neris Sombillo and Ronald Banzon and Cristine Villagonzalo}
}
@inproceedings {acosta,
	title = {Derivation of heat transport equations for a chain of four interacting spins individually coupled to different heat baths},
	booktitle = {Proceedings of the 35th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2017},
	month = {7{\textendash}10 June 2017},
	pages = {SPP-2017-PB-02},
	address = {Cebu City, Philippines},
	abstract = {In this work we considered a linear chain consisting of four (4) spins individually coupled to different heat baths. For this configuration, we used the quantum master equation to derive the heat currents in terms of transition rates. We obtain a set of linear equations that can be solved explicitly to get the transition rates.},
	url = {http://paperview.spp-online.org/proceedings/article/view/42},
	author = {Geronimo Allan Jerome Acosta and Cristine Villagonzalo}
}
@conference {420,
	title = {Fixed-point quantum search as implemented in an Ising spin system},
	year = {2017},
	month = {26{\textendash}29 June 2017},
	pages = {50},
	publisher = {6th Adiabatic Quantum Computing Conference},
	address = {Tokyo, Japan},
	abstract = {The fixed-point quantum search provides a solution to the oscillation problem in a quantum search algorithm through the variation of oracle and inversion phases at each iteration. We simulate the fixed-point algorithm in an Ising spin system with first- and second-nearest neighbor interactions. Since the oracle and the inversion operator acts on O(2), then the oracle qubit is required to be in a superposition of states at the start of the simulation. This provides a way to implement the rotation required in the variation of phases at each iteration. Results show that the probability of success calculated from the algorithm and obtained from the simulation gives a similar behavior which effectively damps the oscillation. We also investigate the effect of tolerance on the fidelity as we increase the number of iterations.},
	url = {http://www.smapip.is.tohoku.ac.jp/~aqc2017/program.html},
	author = {Neris I Sombillo and Ronald S. Banzon and Cristine Villagonzalo}
}
@inproceedings {ferrer,
	title = {Implementing the nonequilibrium Green{\textquoteright}s function method in calculating the phonon transmission in a graphene lattice},
	booktitle = {Proceedings of the 35th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2017},
	month = {7{\textendash}10 June 2017},
	pages = {SPP-2017-1C-06},
	address = {Cebu City, Philippines},
	abstract = {The phonon transmission in a graphene lattice was determined using the nonequilibrium Green{\textquoteright}s function (NEGF) method. The transport system consists of a graphene lattice that was connected to semi-infinite thermal leads. These leads induce a temperature gradient across the central lattice. To achieve ballistic transport, the leads were also chosen to be graphene lattices. Through the NEGF method, the effect of the leads are easily projected to the central lattice via the self-energies Σrα. The decimation method was used to obtain the retarded surface Green{\textquoteright}s functions for the semi-infinite thermal leads grα(ω). After which, the phonon transmission function χ(ω) was calculated using the Caroli formula for the central region. The phonon transmission as a function of ω was obtained for the first and second renormalizations of the Green{\textquoteright}s functions in the decimation method.},
	url = {http://paperview.spp-online.org/proceedings/article/view/145},
	author = {Nica Jane B. Ferrer and Cristine Villagonzalo}
}
@article {377,
	title = {Measurement-enhanced quantum search},
	journal = {AIP Conf. Proc. (8th Jagna International Workshop)},
	volume = {1871},
	year = {2017},
	pages = {020002},
	abstract = {The oscillatory nature of the unitary Grover operator makes quantum searching difficult if one has no prior knowledge of the number of target states in the database. In this work, we coupled the database to an external qubit to transfer the target state in a space where the quantum search operation is not applied. This is followed by a measurement of the external and the ancilla qubit which result to a mixed state with a higher probability of finding the target state. The algorithm is simulated in a spin chain having first and second order neighbor interactions. The result shows that the oscillation is suppressed and the probability of success generally increased with iteration number.},
	doi = {10.1063/1.4996512},
	url = {http://aip.scitation.org/doi/10.1063/1.4996512},
	author = {Neris I Sombillo and Ronald Banzon and Cristine Villagonzalo}
}
@inproceedings {santos,
	title = {The radio-biological bystander effect as a non-linear diffusive phenomenon},
	booktitle = {Proceedings of the 35th Samahang Pisika ng Pilipinas Physics Conference},
	year = {2017},
	month = {7{\textendash}10 June 2017},
	pages = {SPP-2017-3B-03},
	address = {Cebu City, Philippines},
	abstract = {The radiation-induced bystander effect (RIBE) is a phenomenon in which cells undergo genetic damage without having been directly irradiated by an external radiation source. In current literature, it is often modeled via Monte Carlo simulations with the standard Brownian motion as the underlying mode of diffusion. In this study, we replace the standard Brownian motion with the fractional Brownian motion (fBm) and compare it with previous experimental data. We note that the ensemble probability of irradiation is slightly higher in the subdiffusive regime as compared to the standard Brownian motion, which tends to underestimate it. Furthermore, the uniformity of irradiation is preserved in the subdiffusive regime. For these reasons, we conclude that the radio-biological bystander effect is more appropriately examined in the subdiffusive regime (H < 1/2).},
	url = {http://paperview.spp-online.org/proceedings/article/view/74},
	author = {Eduard Renzo Santos and Cristine Villagonzalo}
}
@conference {436,
	title = {The role of orbit to orbit magnetic interaction in a two-dimensional electron gas in the quantum Hall regime},
	year = {2017},
	month = {20{\textendash}24 Nov 2017},
	pages = {P.5},
	publisher = {2017 ICTP Asian Network School and Workshop on Complex Condensed Matter Systems},
	address = {Institute of Physics, Vietnam Academy of Science and Technology, Hanoi, Vietnam},
	abstract = {The orbit to orbit magnetic interaction between electrons in a quantum Hall system is considered. The standard method in the presence of scalar Coulomb interaction potential and external electromagnetic fields is to express the mechanical momentum as p + eA in the effective Hamiltonian. The orbit to orbit interaction is realized by constructing an effective system Hamiltonian possessing a magnetic interaction vector potential in the mechanical momentum. The resulting many-body Schr{\"o}dinger equation with this Hamiltonian is solved to calculate the electron density and potential energy profiles. The aim of this work is to see the effect of the orbit to orbit interaction in calculating the compressibility of the electron gas in the quantum Hall regime.},
	url = {https://www.iop.vast.ac.vn/theor/conferences/asian/2017/program.php},
	author = {Daniel F Marquez and Cristine Villagonzalo}
}
@conference {375,
	title = {Spin-orbit coupled effects on the heat capacity of a two-dimensional electron gas},
	year = {2017},
	month = {4{\textendash}7 Jan 2017},
	publisher = {8th Jagna International Workshop},
	address = {Jagna, Bohol},
	url = {https://www.msuiit.edu.ph/announcements/detail.php?id=1063},
	author = {Cristine Villagonzalo}
}
@inproceedings {316,
	title = {Clustering of irradiated cells due to the bystander effect},
	booktitle = {Proceedings of the 34th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2016},
	pages = {SPP-2016-1C-01},
	address = {University of the Philippines Visayas, Iloilo City},
	abstract = {For quite some time, it has been observed that cells become irradiated despite the fact that no charged particle has ever directly traversed them. This is known as the radiation-induced bystander effect (RIBE). Because RIBE can be characterized as a Markov branching process, some degree of clustering can be expected to occur. By removing the influence of cell positioning on the formation of cancer clusters, we determine that alone, RIBE has a weak, yet non-negligible effect on cluster formation.},
	author = {Eduard Renzo Santos and Cristine Villagonzalo}
}
@conference {340,
	title = {Dependence of Topological Anderson Insulator Phase on the Statistical Distribution of Disorder},
	year = {2016},
	month = {25{\textendash}27 Jan 2016},
	pages = {P1-SC-003},
	publisher = {Annual Meeting of the Physical Society of the Republic of China},
	address = {National Sun Yat-sen University, Kaohsiung},
	abstract = {We consider the dependence to statistical distribution of disorder for attaining disorder-induced topological Anderson insulator (TAI) phase characterized by quantized conductance plateau established in HgTe quantum well simulations. The renormalization of parameters that explains the transition from metal to TAI phase under the effective medium theory proposed by Groth et al. are proportional to the strength of disorder and not explicitly on its statistics. We then suppose that as long as the strength of disorder can be defined, under the appropriate conditions, the system will exhibit the TAI phase. We compare, for instance, the disorder drawn from Gaussian distribution to the case of disorder drawn from a uniform distribution where TAI phase is initially demonstrated. We found that the first case exhibits TAI phase but is not as robust as the uniform distribution case. Inspection of the Gaussian type of disorder is pursued due to vast observation of such distribution in nature. We consider as well the condition when the uniform distribution is discrete, as opposed to being continuous. Lastly, we discuss a possible generalization under the effective medium theory.},
	url = {http://psroc2016.conf.tw},
	author = {Gene Moses M Itable and Cristine Villagonzalo}
}
@inproceedings {336,
	title = {Edge states of the topological Anderson insulator},
	booktitle = {Proceedings of the 34th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2016},
	month = {18{\textendash}21 Aug 2016},
	pages = {SPP-2016-3A-04},
	address = {University of the Philippines Visayas, Iloilo City},
	abstract = {The disorder-induced topological Anderson insulator (TAI) has a similar characteristic to a topological insulator, i.e. insulator in the bulk but conductor in the boundary. In the framework of an effective medium theory, we show that at the transition to TAI from a weaker disorder regime, the conditions for the edge states to exist are appropriately changed yet the physics remains unaltered. This then led us to explore the behavior of edge states at the domain wall of a TAI and a non-TAI through numerical simulations.},
	author = {Gene Moses M Itable and Cristine Villagonzalo}
}
@inproceedings {337,
	title = {k{\textperiodcentered}p Perturbation Model of the Three-terminal Ferromagnet-Insulator-Semiconductor Device with Spin-Orbit Interaction},
	booktitle = {Proceedings of the 34th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2016},
	month = {18{\textendash}21 Aug 2016},
	pages = {SPP-2016-PA-17},
	address = {University of the Philippines Visayas, Iloilo City},
	abstract = {Although the creation of spin polarization in semiconductors in a three-terminal ferromagnet-insulator-normal conductor (3T-FIN) proved successful experimentally, the behavior of spin carriers are not yet well understood. Here we develop a model based on the k{\textperiodcentered}p perturbation theory to provide an insight to the energy behavior of spin carriers in the semiconductor. A potential barrier and with a half-harmonic oscillator model was utilized. It was established that for the model potential, the spin-orbit interaction does not contribute to the additional energy in the 3T-FIN device for both polarized and unpolarized input.},
	author = {April Cortez and Cristine Villagonzalo}
}
@conference {339,
	title = {A Multi-functional Spintronic Device Driven by Spin-orbit Interactions},
	year = {2016},
	month = {25{\textendash}27 Jan 2016},
	pages = {O3-NP1 01},
	publisher = {Annual Meeting of the Physical Society of the Republic of China},
	address = {National Sun Yat-sen University, Kaohsiung},
	abstract = {We propose a simple two-terminal device that operates simultaneously as (i) a spin inverter for an incident polarized electron, (ii) a spin polarizer for an unpolarized source and (iii) a tunable spin current attenuator at certain range of electron energies. The device consists of a U-shaped nanowire with a uniform Dresselhaus spin-orbit coupling and a spatially-modulated Rashba spin-orbit coupling along the curved section. The latter is implemented through a tunable gate voltage at every alternating segment along the wire{\textquoteright}s arc. This multi-purpose device does not need magnetic components or an external magnetic field to carry out the different functions as enumerated above. We discuss how the interplay of the wire{\textquoteright}s geometry and the spin-orbit interaction present in the device drives its spin transport.},
	url = {http://psroc2016.conf.tw},
	author = {Carlos Baldo III and Cristine Villagonzalo}
}
@inproceedings {312,
	title = {Near-resonant approximation in quantum search simulation},
	booktitle = {Proceedings of the 34th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2016},
	month = {18{\textendash}21 Aug 2016},
	pages = {SPP-2016-2C-01},
	address = {University of the Philippines Visayas, Iloilo City},
	abstract = {<div>The simulation of a quantum search algorithm in an Ising spin chain is equivalent to solving a system of 2^<em>L</em> coupled first order differential equations with <em>L</em> being the number of qubits. The near-resonant approximation can be used to simplify the evolution of the states in the database. In this work, we compare the simulation of a quantum search algorithm using the exact evolution according to Schr{\"o}dinger equation and using the near-resonant approximation. Our results show that the error of the near-resonant approximation is comparable to that of the exact evolution in terms of the fidelity measure. The advantage gained in having this small error at the near-resonant approximation is its faster execution time relative to the exact case.</div>},
	author = {Neris D Ilano and Ronald Banzon and Cristine Villagonzalo}
}
@article {baldo2016a,
	title = {Spin inverter and polarizer curved nanowire driven by Rashba and Dresselhaus spin-orbit interactions},
	journal = {Physica E},
	volume = {83},
	year = {2016},
	note = {cited By 0; Article in Press},
	pages = {498-504},
	abstract = {<div style="margin: 0px;">We propose in theory a curved nanowire structure that can both serve as a spin inverter and a spin polarizer driven by a periodic Rashba spin{\textendash}orbit coupling (SOC) and a uniform Dresselhaus SOC. The curved section of the U-shaped quasi-one dimensional nanowire with an arc of radius R and circumferential length πR is divided into segments of equal length initially having only its inherent homogeneous Dresselhaus SOC. Then a Rashba-type SOC is applied at every alternating segment. By tuning the Rashba SOC strength and the incident electron energy, this device can flip the spin at the output of an incoming spin-polarized electron. On the other hand, this same device acts as a spin filter for an unpolarized input for which an outgoing electron with a non-zero polarization can be achieved without the application of an external magnetic field. Moreover, the potential modulation caused by the periodic Rashba SOC enables this device to function as an attenuator for a certain range of incident electron energies that can make the probability current density drop to 10^(-4) of its otherwise magnitude in other regimes.</div>},
	doi = {10.1016/j.physe.2015.11.024},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@inproceedings {320,
	title = {A critically damped quantum search may not necessarily be optimal},
	booktitle = {Proceedings of the AQIS Conference 2015},
	year = {2015},
	month = {24{\textendash}30 Aug 2015},
	pages = {147-148},
	address = {Korea Institute for Advanced Study, Seoul},
	abstract = {We explore the rate at which the target state is obtained in a damped quantum search of an Ising system. As the number of target states is reduced, the damped quantum search gives a higher probability of success than that of the classical search. For any system in general, a new damping parameter was obtained from an optimization of the damped quantum search. Application of the optimum damping parameter per iteration on the damped quantum search shows a significant improvement of the probability of success over the critically damped quantum search.},
	url = {http://aqis-conf.org/2015/program/},
	author = {Neris D Ilano and Ronald Banzon and Cristine Villagonzalo}
}
@conference {350,
	title = {Effect of varying the position of the ancilla qubit and the coupling constants in the simulation of Grover{\textquoteright}s quantum search algorithm},
	year = {2015},
	month = {7{\textendash}11 Jan 2015},
	pages = {P-62},
	publisher = {9th International Conference on Computational Physics},
	address = {National University of Singapore},
	abstract = {We simulate Grover{\textquoteright}s algorithm in an Ising spin chain quantum computer with spatially-dependent coupling constant. The interaction to the nearest and second-nearest neighbor of the spins is considered. The behavior of the fidelity as a function of the number of Grover iteration during the algorithm is determined. For a fixed location of the ancilla qubit, the overall fidelity decreases as we increase the size of the database. Also, as we vary the coupling constant per neighbor, we obtain a high fidelity for a given number of iteration. The simulation of the algorithm using a classical computer system restricts the size of the database that gives a high fidelity result.},
	url = {http://www.physics.nus.edu.sg/iccp9/},
	author = {Neris D Ilano and Cristine Villagonzalo and Ronald Banzon}
}
@inproceedings {cortez-spp-2015,
	title = {Electric field and spin precession effects on spin injection to semiconductors in three-terminal ferromagnet-insulator-nonmagnetic junctions},
	booktitle = {Proceedings of the 33rd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2015},
	month = {3{\textendash}6 Jun 2015},
	pages = {SPP-2015-PA-44},
	address = {University of Northern Philippines, Vigan City},
	abstract = {<div>The application of an electric field and the ensuing spin precession due to internal magnetic fields in the semiconductor of a three-terminal (3T) device is investigated using the standard diffusion theory. Because the externally applied current is confined in one direction in such a device, the effective electric field is assumed to be significant. The dominant precession direction that dictates the spin dynamics is perpendicular to both the injected spins and the electric field. At large effective electric fields, the polarization in the semiconductor is completely drawn from the polarization of the ferromagnet and becomes well-distributed in the semiconductor.</div>},
	author = {April Cortez and Cristine Villagonzalo}
}
@inproceedings {baldo-spp-2015,
	title = {Electron spin polarization along a one-dimensional wire with periodic Rashba and uniform Dresselhaus spin-orbit couplings},
	booktitle = {Proceedings of the 33rd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2015},
	month = {3{\textendash}6 Jun 2015},
	pages = {SPP-2015-3B-05},
	address = {University of Northern Philippines, Vigan City},
	abstract = {<div>We investigate the transmission of an electron spin on a one-dimensional wire with a periodic Rashba spin-orbit couplings (SOC) and a constant contribution from a Dresselhaus SOC. The straight wire is divided into initially uniform segments and the Rashba SOC is turned on every alternating segment. Eigenvalues and eigenfunctions for both without and with Rashba SOC segments are analytically derived. From these obtained values, we numerically calculate the net polarization at the output. We find in our result that regardless of the initial spin orientation, the output polarization turns out to be spin-up polarized. The modulation effect of the segment length of the wire via a periodic Rashba SOC is manifested by the appearance of an energy gap whose width increases with larger Rashba coupling strength. On the other hand, the position of the gap can be tuned by setting the strength of the Dresselhaus SOC. We find that for electrons with Fermi energies falling within this gap, the output spin current is unpolarized.</div>},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@inproceedings {itable-spp-2015,
	title = {Non-equilibrium collective transport in an Ising chain with Kawasaki dynamics},
	booktitle = {Proceedings of the 33rd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2015},
	month = {3{\textendash}6 Jun 2015},
	pages = {SPP-2015-5C-03},
	address = {University of Northern Philippines, Vigan City},
	abstract = {<div>An Ising chain with Kawasaki spin-exchange dynamics in contact with a heat bath with temporal periodic variation between two temperatures is considered. This environment prohibits the system from attaining thermal equilibrium. We demonstrate here via Monte Carlo simulations that such a case exhibits a Brownian motor type collective transport. This confirms our hypothesis that a cycle of cluster and dispersion of domains with fine tuned parameters, in effect, results into spin transport.</div>},
	author = {Gene Moses M Itable and Cristine Villagonzalo}
}
@conference {349,
	title = {Optimal curvilinear transport of a spin-polarized electron with Rashba and Dresselhaus spin-orbit coupling in a uniform magnetic field},
	year = {2015},
	month = {7{\textendash}11 Jan 2015},
	pages = {C3-13},
	publisher = {9th International Conference on Computational Physics},
	address = {National University of Singapore},
	abstract = {The spin-polarized transport of an electron along a curved-one dimensional wire with Rashba and Dresselhaus spin-orbit couplings (SOCs) and a uniform perpendicular magnetic field is considered in this study. Here we analytically derived the appropriate eigenenergies and eigenvalues of the system and numerically determine the output transmission coefficients. Spin polarization, probability current density and conductance are then obtained using transfer matrix approach to examine the dependence of spin transport on the radius, R, of curved wire, on the external magnetic field, B, and on the SOCs strengths.<br>

From our results, we find the condition, B ≈ R^(1/3), at which there would be no spin switching regardless of the SOC strengths. A value of R beyond this condition will no longer have a spin reversal. Moreover, for a material with sufficiently strong SOCs, regardless of its type, the effect of the external magnetic field is reduced. The effective SOC strength follows the standard form equation of a circle in which the Rashba SOC and the Dresselhaus SOC are the component coordinates. From here, we optimize the probability current density by properly tuning the Rashba coupling strength for a given value of B, R and Dresselhaus SOC.},
	url = {http://www.physics.nus.edu.sg/iccp9/},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@inproceedings {ilano-spp-2015,
	title = {Simulation of Grover{\textquoteright}s algorithm with position-variation of the ancilla qubit},
	booktitle = {Proceedings of the 33rd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2015},
	month = {3{\textendash}6 Jun 2015},
	pages = {SPP-2015-2B-06},
	address = {University of Northern Philippines, Vigan City},
	abstract = {<div>We simulate Grover{\textquoteright}s quantum search algorithm in an Ising spin chain quantum computer with nearest and second-nearest neighbor interactions. We determine the optimum Rabi frequency that satisfies the 2πk method by selecting a high fidelity value in the frequency scan for each database size. Also, we vary the position of the ancilla qubit and investigate the effect on the fidelity for different sizes of the database. The behavior of the fidelity as a function of the number of Grover iterations during the algorithm is considered. For a fixed location of the ancilla qubit, the overall fidelity decreases as we increase the size of the database.</div>},
	author = {Neris D Ilano and Ronald S. Banzon and Cristine Villagonzalo}
}
@conference {villagonzalo-qpcmp-2015,
	title = {Spectral and thermodynamic properties of two-dimensional electron systems in the quantum Hall regime},
	year = {2015},
	month = {22{\textendash}23 Apr 2015},
	publisher = {Workshop on Theories in Quantum Phenomena and Condensed Matter Physics},
	address = {University of the Philippines Los Ba{\~n}os, Laguna},
	url = {https://imsp.uplb.edu.ph/index.php/component/content/article/17-imsp-news/research/370-workshop-on-theories-in-quantum-phenomena-and-condensed-matter-physics},
	author = {Cristine Villagonzalo}
}
@conference {355,
	title = {Fundamental studies and simulation of nanostructured materials},
	year = {2014},
	month = {13{\textendash}15 Feb 2014},
	publisher = {6th Asian Computational Materials Design Workshop},
	address = {De La Salle University, Manila},
	url = {http://cmdasiaphilippines.wixsite.com/cmd-asia/6th-cmd-programme-day-1},
	author = {Cristine Villagonzalo}
}
@inproceedings {baldo-spp-2014,
	title = {Magnetoconductance of a spin polarized electron along a curved one-dimensional wire},
	booktitle = {Proceedings of the 32nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2014},
	month = {17{\textendash}20 Oct 2014},
	pages = {SPP2014-1A-6},
	address = {University of the Philippines Diliman, Quezon City},
	abstract = {<div>In this paper, we revisit our transport studies on the motion of a spin-polarized electron along a curved quantum wire under the presence of a uniform perpendicular magnetic field <strong>B</strong> and both the Rashba and Dresselhaus spin-orbit couplings (SOCs). Then from the previously derived output transmission coefficients, we numerically determine the magnetoconductance of the electron motion along the curved wire. Our results agree with our earlier reports that for large effective SOCs, the effects of <strong>B</strong> on spin switching is reduced. Moreover, we found that for a given energy\&nbsp;<em>E</em>, the desired high MC ratio can be achieved by tuning <strong>B</strong>. Finally, with large <em>R</em> the MC reduces significantly due to adiabatic transport.</div>},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@conference {354,
	title = {Simulation of Grover{\textquoteright}s algorithm with spatially-dependent coupling constant},
	year = {2014},
	month = {16{\textendash}25 Sep 2014},
	pages = {19},
	publisher = {Okinawa School in Physics: Coherent Quantum Dynamics},
	address = {Okinawa Institute of Science and Technology, Japan},
	abstract = {We simulate Grover{\textquoteright}s algorithm in an Ising spin chain quantum computer with spatially-dependent coupling constant. The simulation shows that the tendency of two distinct states to transition from the same initial state is reduced. The fidelity with the ideal result and the probability of finding the target state are also investigated.},
	url = {https://groups.oist.jp/cqd/participants-0},
	author = {Neris D Ilano and Cristine Villagonzalo and Ronald Banzon}
}
@article {Baldo201493,
	title = {Spin-orbit coupled transport in a curved quantum wire},
	journal = {Physica E},
	volume = {63},
	year = {2014},
	note = {cited By 0},
	pages = {93-98},
	abstract = {<div>We study the interplay of both Rashba and Dresselhaus spin{\textendash}orbit couplings (SOCs) and a uniform perpendicular magnetic field <strong>B</strong> on the transport of a spin-polarized electron along a curved quantum wire. Eigenenergies and eigenfunctions of the system were analytically solved in the presence of both SOCs for a confinement radius R. From the transmission coefficients, the spin transport properties such as spin polarization, probability current density and spin conductance were computed numerically to determine their dependence on the SOCs, <strong>B</strong> and R. We find the condition for <strong>B</strong> that if it is beyond R<sup>1/3</sup>, no spin reversal will occur. Our results show that for a sufficiently large SOC strength, regardless of its inversion asymmetry origin, the effect of the external magnetic field is reduced. Finding the optimal effective SOC strength is essential in achieving suitable magneto-transport properties for possible spintronic device applications.</div>},
	doi = {10.1016/j.physe.2014.05.013},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@inproceedings {ilano-spp-2014,
	title = {Working precision in a simulation of Grover{\textquoteright}s algorithm in an Ising spin system},
	booktitle = {Proceedings of the 32nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2014},
	month = {17{\textendash}20 Oct 2014},
	pages = {SPP2014-5A-2},
	address = {University of the Philippines Diliman, Quezon City},
	abstract = {<div>We consider the simulation of Grover{\textquoteright}s algorithm in an Ising nuclear spin chain computer with first- and second-nearest neighbor interaction. We calculate the fidelity, a measure of the accuracy of the state, as a function of π/2- and π-pulses and investigate the precision of the success probability of the target state. The fidelity is constrained by the set working precision of the computer system, hence we simulate the algorithm using different working precisions. A reduction of precision is also observed as you increase the number of iterations that restricts the size of the database to effectively simulate the algorithm. The amount of decrease is dependent on the decomposition of pulses.</div>},
	author = {Neris D Ilano and Cristine Villagonzalo and Ronald S. Banzon}
}
@inproceedings {suico-spp-2013,
	title = {Characteristics of non-spanning clusters in a percolating two-dimensional square lattice},
	booktitle = {Proceedings of the 31st Samahang Pisika ng Pilipinas Physics Congress},
	year = {2013},
	month = {23{\textendash}25 Oct 2013},
	pages = {SPP2013-PB-10},
	address = {University of San Carlos, Cebu City},
	abstract = {<div>A study of interest in site percolation is to determine the characteristics of other clusters formed that do not span a given space. In this work we simulate on a square lattice of size <em>L</em> x <em>L</em>\&nbsp;and show the behavior of the mean cluster size and its associated parameters near the percolation threshold. We find that at the percolation threshold the mean cluster size is proportional to <em>L</em>^1.656.</div>},
	author = {Joachim Suico and Ronald S. Banzon and Cristine Villagonzalo}
}
@inproceedings {capili-spp-2013,
	title = {Determination of the percolation critical exponent as measure of dimension},
	booktitle = {Proceedings of the 31st Samahang Pisika ng Pilipinas Physics Congress},
	year = {2013},
	month = {23{\textendash}25 Oct 2013},
	pages = {SPP2013-4C-6},
	address = {University of San Carlos, Cebu City},
	abstract = {<div>The critical exponent\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>\&nbsp;related to the percolation probability describes the dimensionality of a specific lattice topology. The\&nbsp;accurate determination of\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>\&nbsp;will aid in distinguishing different universality classes. This work describes the fine tuning of the\&nbsp;Sigma method developed previously by the authors in obtaining\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>. The approach uses the log-log plot of the percolation\&nbsp;probability versus the difference of the occupation probability and the percolation threshold. Considering a symmetric distribution of points about the midpoint of this plot yields the\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>\&nbsp;with (1) the lowest standard deviation and (2) the\&nbsp;appropriate asymptotic behavior for both one and two-dimensional lattices.</div>},
	author = {Micielle N Capili and Ronald S. Banzon and Cristine Villagonzalo}
}
@inproceedings {ilano-spp-2013,
	title = {The effect of the Hadamard gate decomposition in the simulation of Grover{\textquoteright}s algorithm},
	booktitle = {Proceedings of the 31st Samahang Pisika ng Pilipinas Physics Congress},
	year = {2013},
	month = {23{\textendash}25 Oct 2013},
	pages = {SPP2013-5B-2},
	address = {University of San Carlos, Cebu City},
	abstract = {<div>Simulation of Grover{\textquoteright}s algorithm in an Ising nuclear spin chain computer with first- and second-nearest neighbor interaction is considered. The most frequently used quantum gate in the algorithm is the Hadamard gate. We investigate the effect of its decomposition into several rotation matrices and numerically calculate the fidelity and probability of success of finding the target state. Using an appropriate decomposition, we obtain a high fidelity for different numbers of iterations. The probability of finding the target state also increased and is much closer to the ideal result.</div>},
	author = {Neris D Ilano and Cristine Villagonzalo and Ronald S. Banzon}
}
@conference {356,
	title = {Signatures of Landau level crossings in a two-dimensional electron gas},
	year = {2013},
	month = {1{\textendash}5 Jul 2013},
	pages = {MoP110},
	publisher = {20th International Conference on Electronic Properties of Two-Dimensional Systems},
	address = {Wroc{\l}aw University of Technology, Poland},
	url = {http://www.ep2dsmss.org/},
	author = {Rayda Gammag and Cristine Villagonzalo}
}
@inproceedings {baldo-spp-2013,
	title = {Spin switching on a quantum wire with spin orbit coupling in a uniform magnetic field},
	booktitle = {Proceedings of the 31st Samahang Pisika ng Pilipinas Physics Congress},
	year = {2013},
	month = {23{\textendash}25 Oct 2013},
	pages = {SPP2013-5B-4},
	address = {University of San Carlos, Cebu City},
	abstract = {<div>In this paper, the interplay of a uniform perpendicular magnetic field and the radius of curvature on the spin transport along a curved quantum wire with Rashba and Dresselhaus spin-orbit interactions is investigated. From the analytic calculations of the system{\textquoteright}s eigenenergies and eigenfunctions, transmission coefficients were numerically determined using transfer-matrix approach and were used to calculate the probability current densities and output polarization. Our results show that the spin orientation of the initially spin-polarized electron can be switched for smaller values of radius of curvature of the curved quantum wire. However, the range of spin switching is largely dependent on the threading magnetic field and on the strength of the spin-orbit interactions.</div>},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@article {Gammag201316,
	title = {Two-dimensional electron gas tilt-induced Landau level crossings},
	journal = {Solid State Commun.},
	volume = {156},
	year = {2013},
	note = {cited By 1},
	pages = {16-20},
	abstract = {<div>This work elucidates the novel behavior found in a two-dimensional electron gas (2DEG) under a tilted magnetic field in which the field{\textquoteright}s angle becomes the dominant factor in tuning the spin-splitting rather than the strength of the spin-orbit interaction. The 2DEG eigenvalues are derived with Rashba and Zeeman interactions for various tilt angles and they show crossing-free levels except at very high tilt. Moreover, concomitant with the crossings is the appearance of beats in the 2DEG density of states. The crossings from different levels occur consecutively at around 87{\textdegree}. Similar new observations in Shubnikov{\textendash}de Haas experimental measurements by Hatke et al. (2012) [6] attributed such phenomena to an in-plane-magnetic-field-induced increase in the effective mass. We show here that this behavior is inherent to a 2DEG where spin{\textendash}orbit interaction and the in-plane magnetic field contribution are taken into account.</div>},
	doi = {10.1016/j.ssc.2012.12.006},
	author = {Rayda Gammag and Cristine Villagonzalo}
}
@inproceedings {242,
	title = {The effect of RF distribution on the implementation of quantum CN gate in an Ising spin chain},
	booktitle = {Proceedings of the 30th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2012},
	month = {22{\textendash}24 Oct 2012},
	pages = {SPP-2012-056},
	address = {De La Salle Health Sciences Institute, Dasmari{\~n}as, Cavite},
	abstract = {<div>We consider the dynamics of an Ising spin quantum computer with a radio frequency (RF) pulse distribution. The effect of the frequency spread decreases the resonant transition probability and increases thenear-resonant transition probability. Also, we apply the 2πk method to suppress the near-resonant transition. However, it has an unwanted effect on the probability of the resonant transition. It turns out that the desirable result can be obtained if we use the lowest possible value of k.</div>},
	author = {Neris D Ilano and Cristine Villagonzalo and Ronald S. Banzon}
}
@article {Gammag2012,
	title = {The interplay of Rashba spin-orbit interaction and Landau level broadening on a two-dimensional electron gas under a tilted magnetic field},
	journal = {J. Phys. Conf. Ser.},
	volume = {400},
	number = {PART 4},
	year = {2012},
	note = {cited By 0},
	pages = {042009},
	abstract = {<div>A two-dimensional electron gas in a tilted magnetic field with Rashba spin-orbit interaction (RSOI) was studied. The RSOI is accredited to the asymmetry of the heterostructure where the two-dimensional electron gas is found. The effects of the disorder-attributed Landau level broadening and the RSOI on the spin splitting were identified by simulating the density of states which was assumed to take a Gaussian shape. Increased Landau level broadening obscures the spin splitting and increases the overlap between spin states resulting to stout Gaussian peaks. On the other hand, stronger RSOI amplifies the splitting and lessens the overlap between spin states of the Landau levels. The splitting, however, results to stouter peaks. The similarity in the RSOI and Landau level broadening effects can be explained by recognizing that the asymmetry of the heterostructure is in itself a form of structural disorder.</div>},
	doi = {10.1088/1742-6596/400/4/042009},
	author = {Rayda Gammag and Cristine Villagonzalo}
}
@article {BaldoIII2012,
	title = {Low-lying spin excitations due to next-nearest neighbour interactions in a ferromagnetic lattice with a body-centred cubic (bcc) crystal structure},
	journal = {J. Phys. Conf. Ser.},
	volume = {400},
	number = {PART 3},
	year = {2012},
	note = {cited By 0},
	pages = {032007},
	abstract = {<div>Spin excitations due to the next-nearest-neighbour interactions are investigated for a ferromagnetic lattice with a body-centred cubic (bcc) crystal structure. Using the three-dimensional Heisenberg spin model in a uniform external magnetic field, both the nearest-neighbour (NN) and next-nearest-neighbour (NNN) exchange interactions are considered. The spin wave energies at low temperatures are calculated for which the system is close to its ground state. From the dispersion plots, the additional exchange interactions by the NNN cause an increase in the spin excitation energies. Furthermore, there is a critical NNN/NN exchange ratio for which the effect of the NNN interactions becomes evident as characterized by the presence of two maximum energy peaks in the energy spectra.</div>},
	doi = {10.1088/1742-6596/400/3/032007},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@inproceedings {245,
	title = {Percolation critical exponent of a constrained lattice},
	booktitle = {Proceedings of the 30th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2012},
	month = {22{\textendash}24 Oct 2012},
	pages = {SPP-2012-004},
	address = {De La Salle Health Sciences Institute, Dasmari{\~n}as, Cavite},
	abstract = {<div>We implement the Average Occupation Probability method to treat site percolation on a square lattice with banned sites. It was found that the introduction of banned sites increases the value of the percolation probability while reducing the value of the critical exponent\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>. The deviation of\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>\&nbsp;in site percolation of a constrained lattice from that of a regular system suggests that these\&nbsp;lattice\&nbsp;topologies belong to two different universality classes.</div>},
	author = {Micielle N Capili and Cristine Villagonzalo and Ronald S. Banzon}
}
@article {Gammag2012,
	title = {Persistent spin splitting of a two-dimensional electron gas in tilted magnetic fields},
	journal = {Eur. Phys. J. B},
	volume = {85},
	number = {1},
	year = {2012},
	note = {cited By 1},
	pages = {22},
	abstract = {<div>By varying the orientation of the applied magnetic field with respect to the normal of a two-dimensional electron gas, the chemical potential and the specific heat reveal persistent spin splitting in all field ranges. The corresponding shape of the thermodynamic quantities distinguishes whether the Rashba spin-orbit interaction (RSOI), the Zeeman term or both dominate the splitting. The interplay of the tilting of the magnetic field and RSOI resulted to an amplified splitting even in weak fields. The effects of changing the RSOI strength and the Landau level broadening are also investigated.</div>},
	doi = {10.1140/epjb/e2011-20615-x},
	author = {Rayda Gammag and Cristine Villagonzalo}
}
@article {Gammag2012757,
	title = {Quenching of the DOS beats in a two-dimensional electron gas in tilted magnetic fields},
	journal = {Solid State Commun.},
	volume = {152},
	number = {9},
	year = {2012},
	note = {cited By 2},
	pages = {757-761},
	abstract = {<div>A two-dimensional electron gas exposed to a tilted magnetic field is considered with the Rashba spin{\textendash}orbit interaction and the Zeeman effect. An exact solution for the eigenvalues was obtained assuming that two opposite spin states of adjacent Landau levels have equal probability. No crossings between adjacent eigenenergies were observed, for the tilt angles studied here (θ<=80{\textdegree}), unlike in the perpendicular-magnetic-field case. The absence of crossings leads to quenched beating structures in the oscillations of the density of states (DOS). Persistent spin-splittings were observed at the weak magnetic field region. The splittings, however, can be effectively screened by an increased Landau level broadening. The results shed light on how spins can be controlled through the Rashba interaction strength, the disorder-related broadening and the magnetic field tilt angle.</div>},
	doi = {10.1016/j.ssc.2012.01.047},
	author = {Rayda Gammag and Cristine Villagonzalo}
}
@inproceedings {244,
	title = {Rashba and Dresselhaus spin-orbit interactions in a curvilinear motion of a spin-polarized electron},
	booktitle = {Proceedings of the 30th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2012},
	month = {22{\textendash}24 Oct 2012},
	pages = {SPP-2012-054},
	address = {De La Salle Health Sciences Institute, Dasmari{\~n}as, Cavite},
	abstract = {<div>To see the effects of both Rashba and Dresselhaus spin-orbit interactions on a spin-polarized electron along a curved path, we obtain the system{\textquoteright}s exact energy spectrum taking into account the Zeeman effect. This approximates the motion of an electron in a one-dimensional wire. Our results show that the lifted degeneracy between spins of opposite polarization does not only depend on both Rashba and Dresselhaus coupling constants, but also on the radius of curvature.</div>},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@inproceedings {248,
	title = {Determination of the critical exponent β in site percolation using average occupation probability on a square lattice},
	booktitle = {Proceedings of the 29th Samahang Pisika ng Pilipinas National Physics Congress},
	year = {2011},
	month = {24{\textendash}26 Oct 2011},
	pages = {SPP2011-5A-3},
	address = {University of the Philippines Diliman, Quezon City},
	abstract = {<div>Using the Average Occupation Probability in site percolation, the critical exponent\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>\&nbsp;for a lattice (L {\texttimes} L) is determined. The results for different lattice\&nbsp;sizes show that as the lattice becomes larger, the calculated\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>\&nbsp;agrees better with that obtained from 2D scaling theory. There appears a common\&nbsp;range of p-pc which can be used for any L to compute for\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">β</span>. We estimate that a minimum lattice size of 200 {\texttimes} 200 is necessary for the\&nbsp;procedure to yield an acceptable value of the critical exponent.</div>},
	url = {https://paperview.spp-online.org/proceedings/article/view/3139},
	author = {Micielle N Capili and Cristine Villagonzalo and Ronald S. Banzon}
}
@article {Villagonzalo201143,
	title = {The intrinsic features of the specific heat at half-filled Landau levels of two-dimensional electron systems},
	journal = {J. Low Temp. Phys.},
	volume = {163},
	number = {1-2},
	year = {2011},
	note = {cited By 2},
	pages = {43-52},
	abstract = {<div>The specific heat capacity of a two-dimensional electron gas is derived for two types of the density of states, namely, the Dirac delta function spectrum and that based on a Gaussian function. For the first time, a closed form expression of the specific heat for each case is obtained at half-filling. When the chemical potential is temperature-independent, the temperature is calculated at which the specific heat is a maximum. Here the effects of the broadening of the Landau levels are distinguished from those of the different filling factors. In general, the results derived herein hold for any thermodynamic system having similar resonant states.</div>},
	doi = {10.1007/s10909-010-0259-3},
	author = {Cristine Villagonzalo and Rayda Gammag}
}
@inproceedings {ilano-spp-2011,
	title = {Optimizing the damped quantum search algorithm},
	booktitle = {Proceedings of the 29th Samahang Pisika ng Pilipinas National Physics Congress},
	year = {2011},
	month = {24{\textendash}26 Oct 2011},
	pages = {SPP2011-4A-4},
	address = {University of the Philippines Diliman, Quezon City},
	abstract = {<div>An optimization of the damped quantum search is introduced by calculating the damping parameter that gives the minimum probability of failure per iteration. The resulting damping parameter is compared with the critical damping. The result shows that the optimized damping is generally not equal to the critical damping if the degeneracy of the target state is not comparable to the size of the database. The probability of failure for the optimized damping and critical damping does not differ significantly for the first two iterations. This difference is relatively large for the third iteration. This suggests that for greater number of iterations, the critically damped search is not optimized.</div>},
	url = {https://paperview.spp-online.org/proceedings/article/view/3135},
	author = {Neris D Ilano and Ronald S. Banzon and Cristine Villagonzalo}
}
@inproceedings {249,
	title = {Persistence in the zero temperature dynamics of the 1D Ising model with long range interactions},
	booktitle = {Proceedings of the 29th Samahang Pisika ng Pilipinas National Physics Congress},
	year = {2011},
	month = {24{\textendash}26 Oct 2011},
	pages = {SPP2011-PB-33},
	address = {University of the Philippines Diliman, Quezon City},
	abstract = {<div>We study the zero-temperature dynamics of the 1D Ising Model with long range interactions. In particular, we focus on the persistence behavior as the strength of the long range interaction is varied. Persistence pertains to those spins that did not flip after some time. We observe that whatever the strength of the long range interaction, the system exhibits blocking where some spins do not flip at all. The system size is also varied to see its effect on these properties.</div>},
	url = {https://paperview.spp-online.org/proceedings/article/view/3240},
	author = {Hilton Lazo and Ronald Banzon and Cristine Villagonzalo}
}
@conference {358,
	title = {Long range disorder effects on the heat capacity of two-dimensional electron systems},
	year = {2010},
	month = {12{\textendash}15 Dec 2010},
	pages = {66-67},
	publisher = {9th Asia-Pacific Workshop on Materials Physics},
	address = {Hanoi University of Science, Vietnam},
	abstract = {Broadening of Landau levels of two dimensional electron systems at low temperature and at high perpendicular magnetic fields causes perceptible variations in their interlevel and intralevel excitations as seen in their heat capacity. At a characteristic temperature, the interlevel excitations become negligible and the effect of a small broadening on the energy spectrum vanishes. In this work, a long-range disorder is considered. This can be modeled for the two-dimensional electron gas via a Lorentzian density of states or a Gaussian density of states whose broadening is much larger than kT. The temperature dependence of the heat capacity is derived at half-filled Landau levels or when the chemical potential is a constant. The heat capacity behavior of a system with a long-range disorder is compared to the case of impurity-free quantum Hall systems.},
	url = {http://www.iop.vast.ac.vn/theor/conferences/apwmp/files/Detailed\%20Program.pdf},
	author = {Cristine Villagonzalo and Rayda P Gammag}
}
@inproceedings {baldospp2010,
	title = {Spin waves in a simple ferromagnetic cubic lattice using the Heisenberg model with nearest and next-nearest neighbor interactions},
	booktitle = {Proceedings of the 28th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2010},
	month = {25{\textendash}27 Oct 2010},
	pages = {SPP-2010-006},
	address = {Meralco Management and Leadership Development Center, Antipolo City},
	abstract = {<div>We study the combined effect of the nearest and the next-nearest neighbor interactions on the spin wave dispersion of an isotropic cubic ferromagnet in a uniform magnetic field using the Heisenberg model. Near the groundstate, the components of the spin angular momentum perpendicular to the applied field are negligible as compared to the parallel spin component. Hence, the dispersion relation and the ground state energy estimate of the system are obtained by considering the spin fluctuations along the plane normal to the magnetic field. The results herein are compared to the case where only the effect of the nearest neighbors are present.</div>},
	author = {Carlos F Baldo III and Cristine Villagonzalo}
}
@inproceedings {257,
	title = {Uncertainty relation of the momentum components of a two-dimensional electron system under a tilted magnetic field},
	booktitle = {Proceedings of the 28th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2010},
	month = {25{\textendash}27 Oct 2010},
	pages = {SPP-2010-005},
	address = {Meralco Management and Leadership Development Center, Antipolo City},
	abstract = {<div>A two-dimensional electron gas along the <em>xy</em>-plane is considered. When an applied perpendicular magnetic field, directed along the <em>z</em>-axis, is turned on, the electrons move in cyclotron orbits and the components of the momentum along the <em>x</em>- and <em>y</em>-axes do not commute. This work shows that if the applied magnetic field is tilted at an angle\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">θ</span>\&nbsp;with respect to the <em>z</em>-axis, the uncertainty in simultaneous measurements of the momentum components becomes dependent on\&nbsp;<span style="color: $\#$545454; font-family: arial, sans-serif; font-size: small; line-height: 18.2px;">θ</span>.</div>},
	author = {Rayda P Gammag and Cristine Villagonzalo}
}
@article {gammag2009a,
	title = {Crossover of the specific heat from the quantum Hall to the classical regime},
	journal = {AIP Conf. Proc.},
	volume = {1119},
	year = {2009},
	pages = {213},
	abstract = {<p>The specific heat, <em>C<sub>V</sub></em>, of a two-dimensional electron system in the quantum Hall effect (QHE) regime was calculated numerically. A peak was observed at temperature <em>T<sub>p</sub></em> uniform to all filling factors. This <em>T<sub>p</sub></em> signified the shift of dominance from the extended states to the localized states. After the peak decays <em>C<sub>V</sub></em> rose again, this time reflecting the destruction of the QHE and the emergence of the classical counterpart. The classical Hall was characterized by a <em>T</em>-invariant chemical potential.</p>},
	doi = {http://dx.doi.org/10.1063/1.3137811},
	author = {Rayda Gammag and Cristine Villagonzalo}
}
@inproceedings {ilano-spp-2009,
	title = {Damped quantum search in an Ising spin system},
	booktitle = {Proceedings of the 27th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2009},
	month = {28{\textendash}30 Oct 2009},
	pages = {SPP-2009-006},
	address = {Development Academy of the Philippines Convention Center, Tagaytay City},
	abstract = {<div>We seek on 8 and a 12 spins Ising systems at least one of their eigenstates of a certain eigenvalue by applying the damped quantum search. In particular, without the knowledge of the number of degeneracy, it requires 150 queries to find at least one of the characteristic states having an eigenvalue of -11ε in a database of 4096 items with a probability of 99.1\%. The average number of queries are compared with that of the undamped quantum search. It is found that for a small number of target states, the damped quantum search still preserves the quantum results, while for a large number of target states the classical limit is approached. Moreover, the damped quantum search is shown to increase monotonically as the number of iteration is increased.</div>},
	author = {Neris D Ilano and Ronald S. Banzon and Cristine Villagonzalo}
}
@article {jacosalem2009a,
	title = {Numerical investigations of symmetric embedded quantum well systems using the shooting method and transfer matrix technique},
	journal = {AIP Conf. Proc.},
	volume = {1119},
	year = {2009},
	pages = {219},
	abstract = {<div>Embedded quantum well systems are constructed via a superposition of modified P{\"o}schl-Teller (MPT) potentials. To illustrate the advantages and/or disadvantages of each system, a comparison is made on their energy spectra and scattering properties. The bound states of each system are obtained using the numerical shooting method. The accuracy of this method depends on the step length δz, which can be completely controlled, allowing an improved comparison with analytical solutions to be achieved. A systematic approach of finding the eigenenergies En (n = 0, 1, 2, 3, {\textellipsis}) is also developed in this work. The approach utilizes the established bisection method{\textemdash}a root-finding algorithm that repeatedly divides an interval in half then selects the subinterval in which a root exists. Here selective bisection is implemented for computational efficiency. This entire procedure in finding the solutions of the Schr{\"o}dinger equation is tested for a single MPT well where the solution is analytically known. The result shows a perfect agreement between the numerical and the analytical method. The transfer matrix approach is used to obtain the transmission probability of a particle passing through the embedded quantum well systems. This technique also shows a perfect agreement with the analytical solution.</div>},
	doi = {http://dx.doi.org/10.1063/1.3137825},
	author = {Naomi Jane Jacosalem and Cristine Villagonzalo}
}
@inproceedings {johnson-spp-2009,
	title = {One dimensional migration of DNA in the presence of point particle traps},
	booktitle = {Proceedings of the 27th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2009},
	month = {28{\textendash}30 Oct 2009},
	pages = {SPP-2009-071},
	address = {Development Academy of the Philippines Convention Center, Tagaytay City},
	abstract = {<div>This paper studies the one-dimensional motion of a DNA molecule through a microchannel with a series of point particle electric traps. The three particle molecule model uses three beads connected via harmonic springs in the presence of a Weeks Chandler-Andersen repulsive potential and a Brownian force. Using the stochastic fourth order Runge-Kutta method, the trajectory of the center of mass is observed as the field strengths were varied. System parameters are identified that can impede the flow of a DNA molecule.</div>},
	author = {Sarah Johnson and Cristine Villagonzalo}
}
@inproceedings {gammag-spp-2009,
	title = {Sensitivity of the specific heat to the electronic density of states},
	booktitle = {Proceedings of the 27th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2009},
	month = {28{\textendash}30 Oct 2009},
	pages = {SPP-2009-062},
	address = {Development Academy of the Philippines Convention Center, Tagaytay City},
	abstract = {<div>The effect of an anomalous shape of the electronic density of states near the Fermi energy on the specific heat for bulk solids is investigated in this work. The density of states was simulated using a superposition of the usual square-root-energy-dependent function, and an additional Gaussian function centered in the region about the Fermi energy. The results show an enhanced specific heat when the anomalous density of states{\textquoteright} peak is located after but near the Fermi energy. If this peak is narrow, the specific heat displays an extra maximum with temperature. A large specific heat is suitable for applications wherein it is preferred to keep a system at constant temperature despite a large thermal energy exchange with its environment.</div>},
	author = {Rayda P Gammag and Cristine Villagonzalo}
}
@inproceedings {lazo-spp-2009,
	title = {Site percolation of a fixed lattice size with growing regions of disallowed occupancy},
	booktitle = {Proceedings of the 27th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2009},
	month = {28{\textendash}30 Oct 2009},
	pages = {SPP-2009-036},
	address = {Development Academy of the Philippines Convention Center, Tagaytay City},
	abstract = {<div>We study site percolation in a rectangular lattice with increasing size of two clumped groups of banned sites. The percolation threshold was found to increase as the size of the banned sites is increased. The results were compared to the previously studied lattice with a system with a single group of clumped banned sites with the same number of available sites for occupation. It was observed that the percolation threshold for the dual case is greater than that of the single case. However, the thresholds follow the same trend as the number of available sites is increased.</div>},
	author = {Hilton Lazo and Cristine Villagonzalo and Ronald S. Banzon}
}
@inproceedings {panganiban-spp-2009,
	title = {Thresholds of percolating systems with two symmetric forbidden zones having a constant separation},
	booktitle = {Proceedings of the 27th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2009},
	month = {28{\textendash}30 Oct 2009},
	pages = {SPP-2009-034},
	address = {Development Academy of the Philippines Convention Center, Tagaytay City},
	abstract = {<div>In this paper, we investigate the percolation thresholds for a two-dimensional system with two square banned sites, its topological equivalence with another system with a single square banned site, and the scaling property between their thresholds. An eight-corner check algorithm was implemented to evaluate the percolation of this dual banned sites system. We found that percolation threshold for the dual case follows the same trend as the single case, where the threshold increases as the banned site dimension (<em>X</em>) is increased. Our results suggest that the thresholds for the single and dual case are related by some finite scaling laws.</div>},
	author = {Earl Angelo Panganiban and Cristine Villagonzalo and Ronald S. Banzon}
}
@article {villagonzalo2009a,
	title = {Women in physics in the Philippines: Quantized yet taking steps toward a mature science culture},
	journal = {AIP Conf. Proc.},
	volume = {1119},
	year = {2009},
	pages = {155-156},
	abstract = {<div>Scientific culture in the Philippines is young and physics is no exception. There are only four physics PhD-granting universities with research laboratories. More than 10 universities offer a bachelor{\textquoteright}s degree or master{\textquoteright}s degree in Physics. Like the world trend, these physics institutions are male dominated. However, four of the leading universities already have female PhD faculty members in physics occupying positions of an assistant professor or better. On a positive note, female physicists are no longer limited to work in the national capital region but have carved out their careers in other parts of the country. Also, female physicists have spread into other non-physics-degree-granting universities or found work in the industrial sector. The number of female graduates in physics in the undergraduate and graduate level have slowly but steadily increased since 2002. With the observed increase in number, a working group for women in physics in the Philippines was created this year. In order to provide recommendations to regulators and policy makers, the group{\textquoteright}s first step is to monitor the number of female students and physicists, their study and work environments, and the scholarships and opportunities for development that are available to them.</div>},
	doi = {http://dx.doi.org/10.1063/1.3137751},
	author = {Cristine Villagonzalo and Jinky Bornales and Jelly Grace Betoya-Nonesa}
}
@inproceedings {gammag-spp-2008,
	title = {The 2DES chemical potential: Revealing the link from the classical to the quantum Hall regime},
	booktitle = {Proceedings of the 26th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2008},
	month = {22{\textendash}24 Oct 2008},
	pages = {SPP-2008-026},
	address = {University of the Phlippines Baguio, Baguio City},
	abstract = {<div>A two-dimensional electron system under a magnetic field was studied numerically particularly, the temperature dependence of its chemical potential. For all fillings ν the chemical potential μ saturated to a common value at high temperature T - the classical Hall region. Approaching low T, the quantum Hall regime, μ followed varying paths depending on ν. For the complete fillings, μ showed a knee temperature Tν pointing to the juncture between the two regimes. This crossover, characterized by Tν being constant, indicates the balance between the localization and delocalization effects of the magnetic field and temperature, respectively.</div>},
	author = {Rayda P Gammag and Cristine Villagonzalo}
}
@article {Gammag2008487,
	title = {The interplay of Landau level broadening and temperature on two-dimensional electron systems},
	journal = {Solid State Commun.},
	volume = {146},
	number = {11-12},
	year = {2008},
	note = {cited By 5},
	pages = {487-490},
	abstract = {<div>This work investigates the influence of low temperature and broadened Landau levels on the thermodynamic properties of two-dimensional electron systems. The interplay between these two physical parameters on the magnetic field dependence of the chemical potential, the specific heat and the magnetization is calculated. In the absence of a complete theory that explains the Landau level broadening, experimental and theoretical studies in literature perform different model calculations of this parameter. Here it is presented that different broadening parameters of Gaussian-shaped Landau levels cause width variations in their contributions to interlevel and intralevel excitations. Below a characteristic temperature, the interlevel excitations become negligible. Likewise, at this temperature range, the effect of the Landau level broadening vanishes.</div>},
	doi = {10.1016/j.ssc.2008.03.042},
	author = {Rayda P Gammag and Cristine Villagonzalo}
}
@inproceedings {jacosalem-spp-2008,
	title = {Realignment of energy spectra of modified P{\"o}schl-Teller potentials in a composite quantum well-in-a-well system},
	booktitle = {Proceedings of the 26th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2008},
	month = {22{\textendash}24 Oct 2008},
	pages = {SPP-2008-021},
	address = {University of the Phlippines Baguio, Baguio City},
	abstract = {<div>Quantum well-in-a-well (QWW) systems are constructed via a superposition of modified Poschl-Teller potential wells. The systems are constructed such that the bound states of the component side wells correspond exactly to the excited states of the constituent center well. The energy spectra of these systems are obtained using the shooting method. The contribution of each component well on the total number of bound states of the constructed composite system is presented. Results show that the middle component well is the dominant contributor to the total number of bound states of the QWW system. All contributions are accounted for, however, they are shifted relative to the initial energy level.</div>},
	author = {Naomi Jane Jacosalem and Cristine Villagonzalo}
}
@article {762,
	title = {Realignment of energy spectra of modified P{\"o}schl-Teller potentials in a composite quantum well-in-a-well system},
	journal = {Pisika},
	volume = {2(1)},
	year = {2008},
	pages = {44},
	abstract = {Quantum well-in-a-well (QWW) systems are constructed via a superposition of modified P{\"o}schl-Teller potential wells. The systems are constructed such that the bound states of the component side wells correspond exactly to the excited states of the constituent center well. The energy spectra of these systems are obtained using the shooting method. The contribution of each component well on the total number of bound states of the constructed composite system is presented. Results show that the middle component well is the dominant contributor to the total number of bound states of the QWW system. All contributions are accounted for, however, they are shifted relative to the initial energy level.},
	url = {https://paperview.spp-online.org/pisika-qc/article/view/v2n1p44},
	author = {Naomi Jane Jacosalem and Cristine Villagonzalo}
}
@inproceedings {johnson-spp-2008,
	title = {Simulation of DNA mobility in one-dimension using the bead-spring model},
	booktitle = {Proceedings of the 26th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2008},
	month = {22{\textendash}24 Oct 2008},
	pages = {SPP-2008-138},
	address = {University of the Phlippines Baguio, Baguio City},
	abstract = {<div>This paper studies the one-dimensional motion of a DNA molecule through a microchannel with a single electric trap. The three particle molecule model uses three beads connected via harmonic springs in the presence of a Weeks-Chandler-Andersen repulsive potential and a Brownian force. Using the stochastic fourth order Runge-Kutta method, the spring constant is sought that allows motion through the microchannel without overlap between the three beads.</div>},
	author = {Sarah Johnson and Cristine Villagonzalo}
}
@inproceedings {amarra-spp-2008,
	title = {Simulation of Grover{\textquoteright}s algorithm on a one-dimensional two-level Ising spin system},
	booktitle = {Proceedings of the 26th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2008},
	month = {22{\textendash}24 Oct 2008},
	pages = {SPP-2008-004},
	address = {University of the Phlippines Baguio, Baguio City},
	abstract = {<div>We present an application of Grover{\textquoteright}s search algorithm for quantum computers to the one-dimensional Ising model. In particular, a linear, open ferromagnetic system of four spins is considered, and the algorithm is used to find the six states having energy -ε. By using the Hamiltonian operator, the action of a phase-flipping oracle is defined. An 84\% probability of success is obtained after one iteration.</div>},
	author = {Anthony Amarra and Cristine Villagonzalo and Ronald S. Banzon}
}
@inproceedings {hernal-spp-2008,
	title = {Spanning a space in a quasi-one dimensional lattice vs. two-dimensional lattice with square banned sites},
	booktitle = {Proceedings of the 26th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2008},
	month = {22{\textendash}24 Oct 2008},
	pages = {SPP-2008-163},
	address = {University of the Phlippines Baguio, Baguio City},
	abstract = {<div>Percolation thresholds in quasi- one and two dimensional systems were compared. In the former, one-way spanning is employed with a spanning cluster defined as an array that extends horizontally from one end of the lattice to the other. A four-corner checking procedure is introduced in this study to measure the percolation threshold in the corresponding two-dimensional systems with square banned sites (2DSBS). The results indicate that the threshold values from the two methods converge suggesting that the quasi-1D system and the 2DSBS are equivalent topologically.</div>},
	author = {Mary Joyce Hernal and Cristine Villagonzalo and Ronald S. Banzon}
}
@inproceedings {castro-spp-2007,
	title = {Analysis of symmetric solutions to the Ginzburg-Landau equations for superconducting films},
	booktitle = {Proceedings of the 25th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2007},
	month = {24{\textendash}26 Oct 2007},
	pages = {SPP-2007-177},
	address = {University of the Phlippines Los Ba{\~n}os, Laguna},
	abstract = {<div>Superconductivity in films in a tangential external magnetic field is simulated using a one-dimensional Ginzburg-Landau (GL) model. A symmetric solution to the GL equations is sought using a central difference approximation. This work will illustrate the behavior of superconducting electrons across thick and thin films for various magnetic field strengths. The Meissner effect for type-1 and type II superconductors is verified.</div>},
	author = {Criswelynn D Castro and Cristine Villagonzalo}
}
@inproceedings {cortez-spp-2007,
	title = {Analysis of transient heat conduction on aluminum and iron contacts},
	booktitle = {Proceedings of the 25th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2007},
	month = {24{\textendash}26 Oct 2007},
	pages = {SPP-2007-115},
	address = {University of the Phlippines Los Ba{\~n}os, Laguna},
	abstract = {<div>In this paper, transient temperature profiles for material configurations involving Al-Fe layers in contact with continuous and single pulse heat sources are obtained. A stable and convergent solution is achieved by imposing a suitable Courant-Friedrichs-Lewy (CFL) condition. The magnitude of the effective thermal conductivity of composite materials as compared to their homogeneous counterparts is demonstrated to rise as the number of contacts increases.</div>},
	author = {April D Cortez and Cristine Villagonzalo}
}
@inproceedings {capili-spp-2007,
	title = {Calculation of the average occupation probability for spanning cluster formation in site percolation},
	booktitle = {Proceedings of the 25th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2007},
	month = {24{\textendash}26 Oct 2007},
	pages = {SPP-2007-192},
	address = {University of the Phlippines Los Ba{\~n}os, Laguna},
	abstract = {<div>The average occupation probability for spanning cluster formation in square lattice using site percolation was calculated. Starting with an empty lattice, sites were generated using a pseudo-random number generator. A cluster labeling procedure was implemented for the detection of a spanning cluster. Our results show that the average occupation probability decreases linearly with inverse lattice edge length. The average occupation probability{\textquoteright}s extrapolated value of 0.608 as the lattice size approaches infinity is found to be a good estimate of the accepted value for the site percolation threshold in two dimensions equal to 0.598.</div>},
	author = {Micielle N Capili and Cristine Villagonzalo and Ronald S. Banzon}
}
@inproceedings {gammag-spp-2007,
	title = {The influence of the level broadening on the chemical potential and the heat capacity of two-dimensional electron system},
	booktitle = {Proceedings of the 25th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2007},
	month = {24{\textendash}26 Oct 2007},
	pages = {SPP-2007-179},
	address = {University of the Phlippines Los Ba{\~n}os, Laguna},
	abstract = {<div>Experiments on two dimensional electron systems (2DES) in a magnetic field show a broadening of an otherwise delta-like density of states. This broadening has been attributed as an effect of impurity scattering, electron-electron interactions, electron transfers into a reservoir of electronic states, or disorder, among others. As this broadening is usually present in 2DES, this work incorporates the Landau level broadening to have a more realistic electron system. Its effect on the chemical potential and the heat capacity is simulated. Under strong magnetic fields <em>B</em> and fixed electron concentration, both the chemical potential and the heat capacity are found to oscillate with <em>B</em> at constant temperature.</div>},
	author = {Rayda P Gammag and Cristine Villagonzalo}
}
@inproceedings {villagonzalo-spp-2007,
	title = {Reduction of the thermal conductivity in a 1D superlattice},
	booktitle = {Proceedings of the 25th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2007},
	month = {24{\textendash}26 Oct 2007},
	pages = {SPP-2007-022},
	address = {University of the Phlippines Los Ba{\~n}os, Laguna},
	abstract = {<div>Many technological applications require solids having very low thermal conductivity. This work demonstrates that this can be achieved for specific layer parameters in a superlattice. The phonon dispersion for thick or uneven layered superlattices is modeled using the harmonic approximation in one dimension (1D). The dispersion curves for superlattices with uneven thickness in layers are characterized by the disappearance of forbidden regions. These gaps are usually found in superlattices having uniform thickness. From these results and assuming a constant phonon mean free path, a marked decrease in the lattice contribution to the thermal conductivity is obtained as a function of layer thickness, mass ratio between the layers or asymmetry in its structure.</div>},
	author = {Cristine Villagonzalo and Jane Dizon}
}
@inproceedings {villagonzalo-spp-2006,
	title = {A finite temperature magnetic phase transition in a one-dimensional lattice},
	booktitle = {Proceedings of the 24th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2006},
	month = {25{\textendash}27 Oct 2006},
	pages = {SPP-2006-155},
	address = {Ateneo de Davao University, Davao City},
	abstract = {<div>The effects of a long-range interaction on the magnetization and the heat capacity of a one-dimensional (1D) magnetic ring is studied. In this work, the Monte Carlo method is employed to obtain the equilibrium state of a modified 1D Ising model in a ring geometry having a dipolar interaction of strength <em>G</em> between spins. For a ring consisting of ten spins, the critical temperature at which an ordered state is achieved increases with the absolute value of the dipolar strength. A peak in the heat capacity is observed at finite temperatures. Depending on the magnitude of <em>G</em>, a phase transition at finite temperature may occur.</div>},
	author = {Cristine Villagonzalo and Gina Rose Tongco}
}
@inproceedings {villagonzalo-spp-2006,
	title = {Heat capacity of logistic maps with bidirectional coupling},
	booktitle = {Proceedings of the 24th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2006},
	month = {25{\textendash}27 Oct 2006},
	pages = {SPP-2006-22},
	address = {Ateneo de Davao University, Davao City},
	abstract = {<div>The energy spectra and the heat capacity of two bidirectionally coupled logistics maps obtained from direct numerical simulations are presented. Varying the coupling strength changes the onset of chaos in both maps. In some cases, synchronization in the orbits and heat capacity oscillations of coupled maps are also observed. It is found that the logarithmic-periodic oscillations in the heat capacity of uncoupled logistic maps are no longer observed upon introduction of dissipative coupling.</div>},
	author = {Cristine Villagonzalo and Razel Flores}
}
@inproceedings {gammag-spp-2006,
	title = {Thermodynamic properties of a two-dimensional electron gas in the fractional quantum Hall regime},
	booktitle = {Proceedings of the 24th Samahang Pisika ng Pilipinas Physics Congress},
	year = {2006},
	month = {25{\textendash}27 Oct 2006},
	pages = {SPP-2006-69},
	address = {Ateneo de Davao University, Davao City},
	abstract = {<div>We investigate the chemical potential μ and specific heat <em>CV</em> of a two-dimensional electron gas at temperatures <em>T</em>\&lt;1 K subject to high magnetic fields <em>B</em> that correspond to filling factors <em>ν</em>\&lt;1. We found out that at <em>ν</em>=1/2, the chemical potential is pinned at ℏω_c/2 which is equal to the energy of the first Landau level. The μ(<em>T</em>) increases when <em>ν</em>\&gt;1/2 and decreases when <em>ν</em>\&lt;1/2. This behavior is dictated by the change in the number of states available to the electrons. We also found out that μ(<em>B</em>) converged at ℏω_c/2 and this trend has a negligible <em>T</em>-dependence. When this behavior of the chemical potential was used in calculating the specific heat at the fractional quantum Hall filling factors, <em>CV</em>(<em>T</em>) exhibited a humpy behavior as compared to its smoothly increasing trend at the integral quantum Hall regime. This trend is a direct consequence of the discontinuities in the chemical potential.</div>},
	author = {Rayda P Gammag and Cristine Villagonzalo}
}
@inproceedings {vasquez-spp-2005,
	title = {Constant perturbation in the tight-binding bands of a single-walled carbon nanotube},
	booktitle = {Proceedings of the 23rd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2005},
	month = {26{\textendash}28 Oct 2005},
	pages = {SPP-2005-019},
	address = {Central Philippine University, Iloilo City},
	abstract = {<div>Using a tight-binding calculation, we explain the origin of a large thermopower <em>S</em> in ropes of single-walled carbon nanotubes. We have shown that, by breaking the symmetric energy bands and increasing the overlap energy γ0, a larger <em>S</em> can be obtained. We have also found that a 10\% increase in γ0 resulted to a thermopower that is four times as large.</div>},
	author = {Louella Judy Vasquez and Cristine Villagonzalo}
}
@inproceedings {jaramillo-spp-2005,
	title = {Constrained percolation in two dimensions},
	booktitle = {Proceedings of the 23rd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2005},
	month = {26{\textendash}28 Oct 2005},
	pages = {SPP-2005-130},
	address = {Central Philippine University, Iloilo City},
	abstract = {<div>The effect of implanting square lattice sites of zero occupation probability at the center of a finite but large square lattice for site percolation was investigated. The introduction of symmetrical forbidden site lattices of different sizes inside 40x40 and 60x60 lattices yield percolation threshold values increasing with the size of the forbidden lattice. It was shown that topology in two dimensions dictates the value of the percolation threshold.</div>},
	author = {Christian Oliver Jaramillo and Cristine Villagonzalo and Ronald Banzon}
}
@inproceedings {delatorre-spp-2005,
	title = {Finite difference scheme for hyperbolic heat conduction with continuous and pulsed heat sources},
	booktitle = {Proceedings of the 23rd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2005},
	month = {26{\textendash}28 Oct 2005},
	pages = {SPP-2005-015},
	address = {Central Philippine University, Iloilo City},
	abstract = {<div>The temperature distribution in a region of a one dimensional (1D) semi-infinite slab as defined by the hyperbolic heat equation, is modelled by utilizing the finite-difference approximation for partial-differential equations. The distribution is observed for both pulsed and continuous heat sources for a homogenous medium. The stability of the numerical scheme is determined via the Von-Neumann stability criteria.</div>},
	author = {de la Torre, Agatha and Cristine Villagonzalo}
}
@inproceedings {gammag-spp-2005,
	title = {Magnetic field effects on the chemical potential and specific heat of a 2DEG},
	booktitle = {Proceedings of the 23rd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2005},
	month = {26{\textendash}28 Oct 2005},
	pages = {SPP-2005-040},
	address = {Central Philippine University, Iloilo City},
	abstract = {<div>We investigate the effect of the filling factors on the chemical potential and the specific heat of a two-dimensional electron gas in an external magnetic field. At integer and half-integer filling factors <em>ν</em> the chemical potential is found to jump with an increment of ℏω_c from one Landau level to the next as the magnetic field is increased. In between these special <em>ν</em>{\textquoteright}s the chemical potential tends to infinitesimally increase or decrease in a stepwise manner. This behavior of the chemical potential, in turn, caused the specific heat to vary its behavior with temperature.</div>},
	author = {Rayda P Gammag and Cristine Villagonzalo}
}
@inproceedings {vasquez-spp-2004,
	title = {Diffusive thermoelectric transport in carbon nanotubes},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {119},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {This paper aims to study the diffusion thermoelectric power of singe-walled carbon nanotubes (SWNTs). This was obtained using a linear response method and a semiclassical treatment of metallic-like systems. The energy and the density of states were generated using the tight binding model for both electrons and holes. The SWNT, exhibiting a diffusive transport, is shown to have a negative thermoelectric power for an energy-dependent relaxation time. This indicates that electrons are dominant and demonstrates a metallic-like behavior.},
	author = {Louella Judy Vasquez and Cristine Villagonzalo}
}
@inproceedings {tongco-spp-2004,
	title = {Dipolar interaction in a 1-D Ising ring},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {31},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {<div>The effects of a long-range interaction in a one-dimensional (1D) magnetic ring is studied. In this work, the Monte Carlo method is employed to obtain the equilibrium state of a modified quasi-1D Ising model having a dipolar interaction in a ring geometry. At low temperatures (<em>T</em>), the behavior of the magnetization (<em>M</em>) as a function of <em>T</em> is shown to depend on the strength of the dipolar interaction <em>(G)</em>. For a ring consisting of 10 spins, ferromagnetic states occur in the ground state (<em>T</em>=0) for a small value of <em>G</em>. Antiferromagnetic states are favored when <em>G</em> is large even at a finite temperature. The magnetization saturates at high <em>T</em>. Also, an effective lowering of the total energy is achieved when smaller values of <em>G</em> approach zero and when large value of <em>G</em> is increased further. Similar effects are observed for rings having a larger number of spins.</div>},
	author = {Gina Rose Tongco and Cristine Villagonzalo}
}
@article {Tongco2004,
	title = {Dipolar interaction in a one-dimensional Ising ring},
	journal = {Science Diliman},
	volume = {16},
	year = {2004},
	pages = {8-9},
	abstract = {<div>As synthesis of low-dimensional magnetic systems become possible, the need for theoretical understanding of their behavior arises. In this work, the authors explore a one-dimensional magnetic structure with the spins having dipolar interaction.</div>},
	url = {http://journals.upd.edu.ph/index.php/sciencediliman/article/view/85},
	author = {Gina Rose Tongco and Cristine Villagonzalo}
}
@inproceedings {marco-spp-2004,
	title = {Effect of the hopping parameter on the current through a two-atom molecular wire},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {27},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {<div>We utilize the Hubbard model approach in analyzing the effect of the hopping parameter in the transport properties of a two-atom molecular device. A root finding method was utilized to implement self-consistent calculations to obtain the transport properties and our results agree with that in the literature. Variation of the hopping parameter shows that its value is proportional to the threshold voltage in the presence of an on-site Coulomb interaction <em>U</em>.</div>},
	author = {Ronaldo Marco Jr and Ronald S. Banzon and Cristine Villagonzalo}
}
@article {Morales2004,
	title = {Enhancement of Fe magnetic moments in Fe/Co (001) multilayers},
	journal = {Science Diliman},
	volume = {16},
	year = {2004},
	pages = {6-7},
	abstract = {<div>In order to investigate the electronic and magnetic properties of a bcc Fe/Co (001) multilayer, we have performed electronic structure calculations employing the total energy full-potential linear muffin tin orbital method. The magnetic moments of the layers are calculated. Based on these results, the magnetization profiles and the microscopic origin of the enhancement of Fe moments in the multilayers of the same Co content but with different interface qualities are reported. Large enhancement of magnetic moment is observed in the Fe monolayer located at the interface, and an even greater increase is obtained for the multilayer with one monolayer of intermixing between Fe and Co layers. The Co atoms were found to have similar magnetic moments in the bulk and at the interface.</div>},
	url = {http://journals.upd.edu.ph/index.php/sciencediliman/article/view/84},
	author = {Marienette B Morales and Cristine Villagonzalo}
}
@inproceedings {morales-spp-2004,
	title = {Enhancement of Fe magnetic moments in Fe/Co(001) multilayer},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {72},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {<div>In order to investigate the electronic and magnetic properties of a bcc Fe/Co (001) multilayer, we have performed electronic structure calculations employing the total energy full-potential linear muffin tin orbital method (FP-LMTO). The magnetic moments of the layers are calculated. Based on these results, the magnetization profiles and the microscopic origin of the enhancement of Fe moments in the multilayers of the same Co content but with different interface quality are reported. Large enhancement of magnetic moment is observed in Fe monolayer located at the interface, and even greater increase is obtained for the multilayer with one monolayer of intermixing between Fe and Co layers. The Co atoms were found to have similar magnetic moments in the bulk and at the interface.</div>},
	author = {Marienette Morales and Cristine Villagonzalo}
}
@article {Marco2004,
	title = {Hopping in a molecular wire},
	journal = {Science Diliman},
	volume = {16},
	year = {2004},
	pages = {10-11},
	abstract = {<div>We utilize the Hubbard model approach in analyzing the effect of the hopping parameter in the transport properties of a two-atom molecular device.</div>},
	url = {http://journals.upd.edu.ph/index.php/sciencediliman/article/view/86},
	author = {Ronaldo Marco Jr and Ronald Banzon and Cristine Villagonzalo}
}
@inproceedings {pari{\~n}as-spp-2004,
	title = {Pre-collision dynamics of short 1D anharmonic chains},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {118},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {<div>This paper investigates the dynamics of short one-dimensional lattices having a temperature gradient. An anharmonic potential using the Fermi-Pasta-Ulam (FPU) model was assumed. Three different systems were studied: a uniform chain, a chain with a single impurity, and a diatomic chain. Numerical simulations using the fourth order Runge-Kutta algorithm was used to solve the equations of motion. An impurity in the chain serves as a {\textquotedblleft}thermal barrier{\textquotedblright}. In the diatomic lattice, an energy excitation is formed even before any particle collision occurs. Here the energy distribution at different times will be shown.</div>},
	author = {Ricelle Pari{\~n}as and Cristine Villagonzalo}
}
@inproceedings {gammag-spp-2004,
	title = {The specific heat of a 2DEG in a magnetic field},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {154},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {<div>We investigate the specific heat of a two-dimensional electron gas (2DEG) and how it varies with the system{\textquoteright}s temperature and chemical potential in the presence of an external magnetic field. In this work, we take into account the dependence of the chemical potential with temperature. We show that the specific heat of the system exhibits a two-peak behavior at low T regardless of the field{\textquoteright}s magnitude if μ~T^2 and that this is consistent with results in the literature where the broadening of Landau levels was considered. The underlying processes may be different but the effect on CV is similar.</div>},
	author = {Rayda P Gammag and Cristine Villagonzalo}
}
@inproceedings {dizon-spp-2004,
	title = {Thermal conductivity in a one-dimensional superlattice},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {155},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {<div>We investigate the phonon thermal conductivity in a one-dimensional superlattice for transport perpendicular to the layers by utilizing analytic methods. A reduction on the thermal conductivity was observed when the layer thickness of the superlattice was increased. Difference in the mass ratio of the atoms in the chain and force constant ratio on alternating layers produced variations in the thermal conductivity. Higher difference of either of these two parameters results in the decrease of the thermal conductivity.</div>},
	author = {Jane Dizon and Cristine Villagonzalo}
}
@inproceedings {muldera-spp-2004,
	title = {Thermal efficiency of a quasi-one-dimensional conductor},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {113},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {<div>How does a system{\textquoteright}s dimension affect its thermal efficiency? In this work, we will show that if conduction is constrained in a single direction in a three dimensional conductor, it will become more efficient in conducting thermal energy as its width is reduced. Herein, we investigate the electric and thermal properties of a quasi-one-dimensional conductor using the free electron gas model in the relaxation time approximation scheme. This is demonstrated in the wire width dependence of the electric conductivity, thermal conductivity and figure of merit. We will also show the dependence of these transport properties to the system{\textquoteright}s chemical potential.</div>},
	author = {Joselito Muldera and Cristine Villagonzalo}
}
@inproceedings {sandagon-spp-2004,
	title = {Variational Monte Carlo of a 2D quantum well system},
	booktitle = {Proceedings of the 22nd Samahang Pisika ng Pilipinas Physics Congress},
	year = {2004},
	month = {25{\textendash}27 Oct 2004},
	pages = {153},
	address = {Bohol Tropics Resort, Tagbilaran City},
	abstract = {<div>The convergence of a Variational Monte Carlo method with Gaussian importance sampling (VMC-GIS), as applied to a two-dimensional time-independent quantum well, was improved by optimizing the standard deviation σ of the Gaussian probability distribution. The fastest convergence in attaining the ground state is obtained by using VMC-GIS with σ=0.6 .</div>},
	author = {Ryan Sandagon and Ronald S. Banzon and Cristine Villagonzalo}
}
@inproceedings {labora-spp-2003,
	title = {Ab-initio structural investigation of Fe/V (001) multilayers},
	booktitle = {Proceedings of the 21st Samahang Pisika ng Pilipinas Physics Congress},
	year = {2003},
	month = {22{\textendash}25 Oct 2003},
	pages = {SPP-2003-1D-02},
	address = {University of San Carlos, Cebu City},
	abstract = {The structural properties of Fe/V multilayers are studied using electronic structure calculations. A first principle numerical calculation is implemented using the Full Potential Linear Muffin-Tin (FPLMTO) method based on the Density Functional Theory (DFT). The number of Fe and V in a Fe/V supercell is varied while the total atoms is kept constant. The energetically equilibrium structure of this multilayer is presented. Two types of multilayer structures, the Fe-on-V and the V-on-Fe, are investigated and compared to bulk values and reported experimental results. Results show that the total energy for Fe-on-V multilayer gives a lower value compared to V-on-Fe multilayer signifying that the V-on-Fe is a more stable structure than Fe-on-V, a result which is consistent with experimental observations. The simple case of the effect of intermixing is presented and its total energy and spin moment is compared to the bulk value.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2003-1D-02},
	author = {Marites Labora and Cristine Villagonzalo}
}
@inproceedings {morales-spp-2003,
	title = {Magnetization and magnetic anisotropy of Fe/Co multilayers: A first-principles calculation},
	booktitle = {Proceedings of the 21st Samahang Pisika ng Pilipinas Physics Congress},
	year = {2003},
	month = {22{\textendash}25 Oct 2003},
	pages = {SPP-2003-PB-05},
	address = {University of San Carlos, Cebu City},
	abstract = {The magnetic properties of bcc Fe/Co (001) multilayers are calculated from first-principles by means of the full potential linear muffin tin orbital method (FLMTO). A significant enhancement in the magnetic moments of each atom in the multilayer is observed. The magnetization and magnetic anisotropy energy MAE are described as a function of Fe/Co layer thickness ratio. For comparison with the obtained multilayer values, the magnetization and MAE of constituent atoms, bcc Co bulk are also verified. The results for the magnetization and MAE agree well with experimental results.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2003-PB-05},
	author = {Marienette B Morales and Cristine Villagonzalo}
}
@inproceedings {villagonzalo-spp-2003,
	title = {Numerical investigations on the growth textures of Co/Cu multilayers},
	booktitle = {Proceedings of the 21st Samahang Pisika ng Pilipinas Physics Congress},
	year = {2003},
	month = {22{\textendash}25 Oct 2003},
	pages = {SPP-2003-3D-01},
	address = {University of San Carlos, Cebu City},
	abstract = {A first-principles investigation on the growth textures of cobalt/copper multilayers is presented. Through an ab-initio local-density approximation in a full-potenatial linear muffin-tin orbital (FPLMTO) electronic structure method, the different growth textures of Co/Cu multilayers are compared. Through total energy minimization, the face-centered cubic (111) configuration is determined to be the preferred textural orientation as compared to the (110) and (001) directions. Physically, this suggests that (111) configuration is more favorable to grow. Included in this report are results of intermixing Co/Cu interfaces which verifies that Co and Cu are immiscible in the bulk. This work provides a benchmark for thermal transport studies in multilayers.},
	url = {https://proceedings.spp-online.org/article/view/SPP-2003-3D-01},
	author = {Cristine Villagonzalo and Arun K Setty and Bernard R Cooper}
}
@article {R{\"o}mer2003167,
	title = {Thermoelectric properties of disordered systems},
	journal = {J. Phys. Soc. Jpn. Suppl. A},
	volume = {72},
	number = {SUPPL. A},
	year = {2003},
	pages = {167-168},
	abstract = {<div>The electronic properties of disordered systems have been the subject of intense study for several decades. Thermoelectric properties, such as thermopower and thermal conductivity, have been relatively neglected. A long standing problem is represented by the sign of the thermoelectric power. In crystalline semiconductors this is related to the sign of the majority carriers, but in non-crystalline systems it is commonly observed to change sign at low temperatures. In spite of its apparent universality this change has been interpreted in a variety of ways in different systems. We have developed a Green{\textquoteright}s function recursion algorithm based on the Chester-Thelling-Kubo-Greenwood formula for calculating the kinetic coefficients on long strips or bars. From these we can deduce the electrical conductivity, the Seebeck and Peltier coefficients and the thermal conductivity, as well as the Lorenz number. We present initial results for 1D systems. We observe a Lorentzian distribution for the thermopower which is modified by the presence of inelastic scattering. This could give rise to non-negligible quantum fluctuations in macroscopic systems at low temperatures.</div>},
	doi = {10.1143/JPSJS.72SA.167},
	author = {Rudolf A R{\"o}mer and Angus MacKinnon and Cristine Villagonzalo}
}
@article {Villagonzalo200016446,
	title = {Behavior of the thermopower in amorphous materials at the metal-insulator transition},
	journal = {Phys. Rev. B},
	volume = {62},
	number = {24},
	year = {2000},
	note = {cited By 4},
	pages = {16446-16452},
	abstract = {<div>We study the behavior of the thermal transport properties in three-dimensional disordered systems close to the metal-insulator transition within linear response. Using a suitable form for the energy-dependent conductivity σ, we show that the value of the dynamical scaling exponent for noninteracting disordered systems such as the Anderson model of localization can be reproduced. Furthermore, the values of the thermopower <em>S</em> have the right order of magnitude close to the transition as compared to the experimental results. A sign change in the thermoelectric power <em>S</em> {\textemdash} as is often observed in experiments {\textemdash} can also be modeled within the linear response formulation using modified experimental σ data as input.</div>},
	doi = {10.1103/PhysRevB.62.16446},
	author = {Cristine Villagonzalo and Rudolf A R{\"o}mer and Michael Schreiber and Angus MacKinnon}
}
@article {Villagonzalo1999179,
	title = {Thermoelectric transport properties in disordered systems near the Anderson transition},
	journal = {Eur. Phys. J. B},
	volume = {12},
	number = {2},
	year = {1999},
	note = {cited By 12},
	pages = {179-189},
	abstract = {<div>We study the thermoelectric transport properties in the three-dimensional Anderson model of localization near the metal-insulator transition (MIT). In particular, we investigate the dependence of the thermoelectric power <em>S</em>, the thermal conductivity <em>K</em>, and the Lorenz number <em>L</em>0 on temperature <em>T</em>. We first calculate the <em>T</em> dependence of the chemical potential μ from the number density <em>n</em> of electrons at the MIT using an averaged density of states obtained by diagonalization. Without any additional approximation, we determine from the behavior of <em>S</em>, <em>K</em> and <em>L</em>0 at low <em>T</em> as the MIT is approached. We find that the d.c. conductivity and <em>K</em> decrease to zero at the MIT as <em>T</em> -\&gt; 0 and show that <em>S</em> does not diverge. Both <em>S</em> and <em>L</em>0 become temperature independent at the MIT and depend only on the critical behavior of the conductivity.</div>},
	doi = {10.1007/s100510050994},
	author = {Cristine Villagonzalo and Rudolf A R{\"o}mer and Michael Schreiber}
}
@article {Villagonzalo1998394,
	title = {Low temperature behavior of the thermopower in disordered systems near the Anderson transition},
	journal = {Ann. Phys. (Leipzig)},
	volume = {7},
	number = {5-6},
	year = {1998},
	note = {cited By 5},
	pages = {394-399},
	abstract = {<div>We investigate the behavior of the thermoelectric power <em>S</em> in disordered systems close to the Anderson-type metal-insulator transition (MIT) at low temperatures. In the literature, we find contradictory results for <em>S</em>. It is either argued to diverge or to remain a constant as the MIT is approached. To resolve this dilemma, we calculate the number density of electrons at the MIT in disordered systems using an averaged density of states obtained by diagonalizing the three-dimensional Anderson model of localization. From the number density we obtain the temperature dependence of the chemical potential necessary to solve for <em>S</em>. Without any additional approximation, we use the Chester-Thellung-Kubo-Greenwood formulation and numerically obtain the behavior of <em>S</em> at low <em>T</em> as the Anderson transition is approached from the metallic side. We show that indeed <em>S</em> does not diverge.</div>},
	doi = {10.1002/(SICI)1521-3889(199811)7:5/6<394::AID-ANDP394>3.0.CO;2-Q},
	author = {Cristine Villagonzalo and Rudolf A R{\"o}mer}
}
