Miroslav Kolesik
 Professor, Optical Sciences
 Member of the Graduate Faculty
Contact
 (520) 6214602
 Meinel Optical Sciences, Rm. 538
 Tucson, AZ 85721
 kolesik@acms.arizona.edu
Bio
No activities entered.
Interests
No activities entered.
Courses
202425 Courses

The Beam Propagation Method
OPTI 547 (Fall 2024) 
Thesis
OPTI 910 (Fall 2024)
202324 Courses

Nonlinear lightmatter interac
OPTI 551 (Spring 2024) 
Thesis
OPTI 910 (Spring 2024) 
The Beam Propagation Method
OPTI 547 (Fall 2023) 
Thesis
OPTI 910 (Fall 2023)
202223 Courses

Dissertation
OPTI 920 (Spring 2023) 
Dissertation
PHYS 920 (Spring 2023) 
Nonlinear lightmatter interac
OPTI 551 (Spring 2023) 
Dissertation
OPTI 920 (Fall 2022)
202122 Courses

Dissertation
OPTI 920 (Spring 2022) 
Dissertation
PHYS 920 (Spring 2022) 
Master's Report
OPTI 909 (Spring 2022) 
Physical Optics II
OPTI 330 (Spring 2022) 
Dissertation
OPTI 920 (Fall 2021) 
Dissertation
PHYS 920 (Fall 2021) 
The Beam Propagation Method
OPTI 547 (Fall 2021) 
Thesis
OPTI 910 (Fall 2021)
202021 Courses

Dissertation
OPTI 920 (Spring 2021) 
Master's Report
OPTI 909 (Spring 2021) 
Nonlinear lightmatter interac
OPTI 551 (Spring 2021) 
Dissertation
OPTI 920 (Fall 2020) 
The Beam Propagation Method
OPTI 547 (Fall 2020) 
Thesis
OPTI 910 (Fall 2020)
201920 Courses

Master's Report
OPTI 909 (Spring 2020) 
Physical Optics I
OPTI 210 (Spring 2020) 
The Beam Propagation Method
OPTI 547 (Fall 2019)
201819 Courses

Directed Graduate Research
OPTI 792 (Spring 2019) 
Physical Optics I
OPTI 210 (Spring 2019) 
Directed Graduate Research
OPTI 792 (Fall 2018) 
Physical Optics I
OPTI 310 (Fall 2018) 
The Beam Propagation Method
OPTI 547 (Fall 2018)
201718 Courses

Dissertation
OPTI 920 (Spring 2018) 
Independent Study
OPTI 599 (Spring 2018) 
Dissertation
OPTI 920 (Fall 2017) 
Mathematical Optics Lab
OPTI 512L (Fall 2017) 
Physical Optics I
OPTI 310 (Fall 2017)
201617 Courses

Independent Study
OPTI 599 (Summer I 2017) 
Dissertation
OPTI 920 (Spring 2017) 
Independent Study
OPTI 599 (Spring 2017) 
Thesis
OPTI 910 (Spring 2017) 
Dissertation
OPTI 920 (Fall 2016) 
Physical Optics I
OPTI 310 (Fall 2016) 
Thesis
OPTI 910 (Fall 2016)
201516 Courses

Dissertation
OPTI 920 (Spring 2016) 
The Beam Propagation Method
OPTI 547 (Spring 2016)
Scholarly Contributions
Chapters
 Panagiotopoulos, P., Whalen, P. T., Kolesik, M., & Moloney, J. V. (2016). Numerical Simulation of UltraShort Laser Pulses. In Laser Filamentation(pp 185213). Springer International Publishing.
Journals/Publications
 Heinz, J., & Kolesik, M. (2023). Open boundary conditions for the simulation of leaky modes. Optics Continuum, 2(1), 143154.
 Kolesik, M. (2023). Assessment of tightbinding models for highharmonic generation in zinc blende materials. Optics Letters, 48(12), 31913194.
 Kolesik, M. (2023). Which part of the Brillouin zone contributes most to the highharmonic radiation?. arXiv preprint arXiv:2309.07281.
 Kolesik, M., & Moloney, J. V. (2023). Numerical discreteness and dephasing in highharmonic calculations in solids. Physical Review B, 108(11), 115433.
 Panagiotopoulos, P., Heinz, J., Kolesik, M., & Moloney, J. V. (2023). Higher order longwave infrared Bessel beams generating tubular plasma waveguides in the atmosphere. Optics Continuum, 2(2), 449455.
 Gu, J., & Kolesik, M. (2022). FullBrillouinzone calculation of highorder harmonic generation from solidstate media. Physical Review A, 106(6), 063516.
 Gu, J., & Kolesik, M. (2022). FullBrillouinzone calculation of higherorder harmonic radiation from solidstate media. Phys. Rev. A, 106, 063516.
 Hastings, M. G., Panagiotopoulos, P., Kolesik, M., Hasson, V., Tochitsky, S., & Moloney, J. V. (2022). Fewcycle 10 $\mu$m multiterawatt pulse selfcompression in a gasfilled multipass cell: a numerical experiment. JOSA B, 39(1), 266272.
 Kolesik, M., Panagiotopoulos, P., & Moloney, J. V. (2022). Nonlinear localization of high energy long wave laser pulses in fully correlated 3D turbulence. Optics Letters, 47(7), 17821785.
 Yusofsani, S., & Kolesik, M. (2022). Exact energy eigenstates of the CoulombStark Hamiltonian. Physical Review A, 105(3), 032203.
 Gu, J., Schweinsberg, A., Vanderhoef, L., Tripepi, M., Valenzuela, A., Wolfe, C., Ensley, T. R., Chowdhury, E., & Kolesik, M. (2021). Random quasiphasematching in polycrystalline media and its effects on pulse coherence properties. Optics Express, 29(5), 74797493.
 Heinz, J., & Kolesik, M. (2021). Transparent boundary conditions for discretized nonHermitian eigenvalue problems. International Journal of Modern Physics C, 32(10), 2150138.
 Le, M., Salehi, F., Railing, L., Kolesik, M., & Milchberg, H. (2021). 15 MeV quasimonoenergetic electrons at 1 kHz with circularly polarized fewcycle pulses. Bulletin of the American Physical Society, 66.
 Panagiotopoulos, P., Kolesik, M., Hasson, V., Tochitsky, S., & Moloney, J. V. (2021). Numerical study of spatial propagation dynamics and energy delivery of TW squareaperture CO 2 laser pulses in the atmosphere. JOSA B, 38(4), 12141221.
 Panagiotopoulos, P., Kolesik, M., Tochitsky, S., & Moloney, J. V. (2021). Generation of long homogeneous plasma channels with high power longwave IR pulsed Bessel beams. Optics letters, 46(21), 54575460.
 Salehi, F., Le, M., Railing, L., Kolesik, M., & Milchberg, H. M. (2021). Laseraccelerated, lowdivergence 15mev quasimonoenergetic electron bunches at 1 khz. Physical Review X, 11(2), 021055.
 Vasilyev, S., Gu, J., Mirov, M., Barnakov, Y., Moskalev, I., Smolski, V., Peppers, J., Kolesik, M., Mirov, S., & Gapontsev, V. (2021). Lowthreshold supercontinuum generation in polycrystalline media. JOSA B, 38(5), 16251633.
 Yusofsani, S., & Kolesik, M. (2021). Beyond FowlerNordheim model: harmonic generation from metallic nanostructures. The European Physical Journal Special Topics, 19.
 Gu, J., Hastings, M. G., & Kolesik, M. (2020). Simulation of harmonic and supercontinuum generation in polycrystalline media. JOSA B, 37(5), 15101517.
 Heinz, J., Reiff, R., Joyce, T., Becker, A., Jaro\'nBecker, A., & Kolesik, M. (2020). Nonlinear polarization and ionization in O 2: metastable electronic state model. Optics Express, 28(17), 2502225036.
 Kilen, I., Kolesik, M., Hader, J., Moloney, J. V., Huttner, U., Hagen, M. K., & Koch, S. W. (2020). Propagation Induced Dephasing in Semiconductor HighHarmonic Generation. Physical review letters, 125(8), 083901.
 Panagiotopoulos, P., Hastings, M. G., Kolesik, M., Tochitsky, S., & Moloney, J. V. (2020). Multiterawatt femtosecond 10 $\mu$m laser pulses by selfcompression in a CO 2 cell. OSA Continuum, 3(11), 30403047.
 Tolliver, J., Zahedpour, S., Wahlstrand, J. K., Milchberg, H. M., & Kolesik, M. (2020). Nonlinearity and ionization in Xe: experimentbased calibration of a numerical model. Optics Letters, 45(20), 57805783.
 Yusofsani, S., & Kolesik, M. (2020). Quantum tunneling time: Insights from an exactly solvable model. Physical Review A, 101(5), 052121.
 Eftekhar, M. A., SanjabiEznaveh, Z., LopezAviles, H. E., Benis, S., AntonioLopez, J. E., Kolesik, M., Wise, F., AmezcuaCorrea, R., & Christodoulides, D. N. (2019). Accelerated nonlinear interactions in gradedindex multimode fibers. Nature communications, 10(1), 110.
 Kolesik, M., & Wright, E. M. (2019). Universal longwavelength nonlinear optical response of noble gases. Optics express, 27(18), 2544525456.
 Panagiotopoulos, P., Kolesik, M., Koch, S. W., Wright, E. M., Tochitsky, S., & Moloney, J. V. (2019). Control of the filament dynamics of 10 $\mu$m pulses via designer pulse trains. JOSA B, 36(10), G33G39.
 Panagiotopoulos, P., Kolesik, M., Tochitsky, S., Koch, S. W., & Moloney, J. V. (2019). Twostage filamentation of 10 $\mu$m pulses as a broadband infrared backlighter in the atmosphere. Optics letters, 44(12), 31223125.
 Ren, X., Wang, Y., Chang, Z., Welch, J., Bernstein, A., Downer, M., Brown, J., Gaarde, M., Couairon, A., Kolesik, M., & others, . (2019). Inline Spectral Interferometry in ShortwaveInfrared Laser Filaments in Air. Physical Review Letters, 123(22), 223203.
 Rosenow, P., Kolesik, M., Koch, S. W., & Moloney, J. V. (2019). Effective nonlinear rovibrational response of water vapor for efficient pulse propagation simulations. JOSA B, 36(2), 267274.
 Rosenow, P., Panagiotopoulos, P., Kolesik, M., Koch, S. W., & Moloney, J. V. (2019). Nonlinear rovibrational response in the propagation of longwavelength infrared pulses and pulse trains. JOSA B, 36(12), 34573463.
 Tochitsky, S., Welch, E., Panagiotopoulos, P., Polyanskiy, M., Pogorelsky, I., Kolesik, M., Koch, S. W., Moloney, J. V., & Joshi, C. (2019). Filamentation of longwave infrared pulses in the atmosphere. JOSA B, 36(10), G40G51.
 Tochitsky, S., Welch, E., Polyanskiy, M., Pogorelsky, I., Panagiotopoulos, P., Kolesik, M., Wright, E. M., Koch, S. W., Moloney, J. V., Pigeon, J., & others, . (2019). Megafilament in air formed by selfguided terawatt longwavelength infrared laser. Nature Photonics, 13(1), 4146.
 Werner, K., Hastings, M. G., Schweinsberg, A., Wilmer, B. L., Austin, D., Wolfe, C. M., Kolesik, M., Ensley, T. R., Vanderhoef, L., Valenzuela, A., & others, . (2019). Ultrafast midinfrared high harmonic and supercontinuum generation with n 2 characterization in zinc selenide. Optics express, 27(3), 28672885.
 Wright, E. M., Koch, S. W., Kolesik, M., & Moloney, J. V. (2019). Memory effects in the longwave infrared avalanche ionization of gases: a review of recent progress. Reports on Progress in Physics, 82(6), 064401.
 Jakobsen, P., Mansuripur, M., & Kolesik, M. (2018). Leakymode expansion of the electromagnetic field inside dispersive spherical cavity. Journal of Mathematical Physics, 59(3), 033501.
 Juhasz, D., Kolesik, M., & Jakobsen, P. K. (2018). Convergence and completeness for squarewell Stark resonant state expansions. Journal of Mathematical Physics, 59(11), 113501.
 Tochitsky, S., Welch, E., Polyanskiy, M., Pogorelsky, I., Panagiotopoulos, P., Kolesik, M., Wright, E., Koch, S., Moloney, J., Pigeon, J., & others, . (2018). Observation of arrest of Kerr selffocusing in a 10 $\mu$m filament in air at 1 TW/cm 2 clamped intensity. Bulletin of the American Physical Society.
 Tolliver, J., & Kolesik, M. (2018). True versus effective Kerr nonlinear response in optical filamentation. Optics express, 26(23), 3017230182.
 Wahlstrand, J. K., Zahedpour, S., Bahl, A., Kolesik, M., & Milchberg, H. M. (2018). Boundelectron nonlinearity beyond the ionization threshold. Physical review letters, 120(18), 183901.
 Bahl, A., Majety, V. P., Scrinzi, A., & Kolesik, M. (2017). Nonlinear optical response in molecular nitrogen: from abinitio calculations to optical pulse simulations. Optics letters, 42(12), 22952298.
 Bahl, A., Wahlstrand, J. K., Zahedpour, S., Milchberg, H. M., & Kolesik, M. (2017). Nonlinear optical polarization response and plasma generation in noble gases: Comparison of metastableelectronicstateapproach models to experiments. Physical Review A, 96(4), 043867.
 Eftekhar, M. A., Wright, L. G., Mills, M. S., Kolesik, M., Correa, R. A., Wise, F. W., & Christodoulides, D. N. (2017). Versatile supercontinuum generation in parabolic multimode optical fibers. Optics express, 25(8), 90789087.
 Eznaveh, Z. S., Eftekhar, M. A., Lopez, J. A., Kolesik, M., Sch\"ulzgen, A., Wise, F. W., Christodoulides, D. N., & Correa, R. A. (2017). Tailoring frequency generation in uniform and concatenated multimode fibers. Optics letters, 42(5), 10151018.
 Mansuripur, M., Kolesik, M., & Jakobsen, P. (2017). Leaky modes of solid dielectric spheres. Physical Review A, 96(1), 013846.
 Schuh, K., Kolesik, M., Wright, E. M., Moloney, J. V., & Koch, S. W. (2017). Selfchanneling of highpower longwave infrared pulses in atomic gases. Physical review letters, 118(6), 063901.
 Schuh, K., Panagiotopoulos, P., Kolesik, M., Koch, S. W., & Moloney, J. V. (2017). Multiterawatt 10 $\mu$m pulse atmospheric delivery over multiple Rayleigh ranges. Optics letters, 42(19), 37223725.
 Schuh, K., Rosenow, P., Kolesik, M., Wright, E. M., Koch, S. W., & Moloney, J. V. (2017). Nonlinear rovibrational polarization response of water vapor to ultrashort longwave infrared pulses. Physical Review A, 96(4), 043818.
 Wang, T., & Kolesik, M. (2017). On the manifestation of higherorder nonlinearities in a noble gas medium undergoing strong ionization. Optics letters, 42(20), 41954198.
 Bahl, A., Wright, E. M., & Kolesik, M. (2016). Nonlinear optical response of noble gases via the metastable electronic state approach. Physical Review A, 94(2), 023850.
 Brown, J. M., Jakobsen, P., Bahl, A., Moloney, J. V., & Kolesik, M. (2016). On the convergence of quantum resonantstate expansion. Journal of Mathematical Physics, 57(3), 032105.
 Brown, J. M., Shanor, C., Wright, E. M., & Kolesik, M. (2016). Carrierwave shape effects in optical filamentation. Optics letters, 41(5), 859862.
 Mansuripur, M., Kolesik, M., & Jakobsen, P. (2016). Leaky Modes of Dielectric Cavities. Spintronics IX, edited by H.J. Drouhin, J.E. Wegrowe, and M. Razeghi, Proceedings of SPIE, 9931, 99310B~1:20.
 Panagiotopoulos, P., Couairon, A., Kolesik, M., Papazoglou, D. G., Moloney, J. V., & Tzortzakis, S. (2016). Nonlinear plasmaassisted collapse of ringAiry wave packets. Physical Review A, 93(3), 033808.
 Panagiotopoulos, P., Kolesik, M., & Moloney, J. V. (2016). Exploring the limits to energy scaling and distanttarget delivery of highintensity midinfrared pulses. Physical Review A, 94(3), 033852.
 Panagiotopoulos, P., Schuh, K., Kolesik, M., & Moloney, J. V. (2016). Simulations of 10 $\mu$m filaments in a realistically modeled atmosphere. JOSA B, 33(10), 21542161.
 Shanor, C., Ensley, T., Hagan, D. J., Van Stryland, E. W., Wright, E. M., & Kolesik, M. (2016). Numerical investigation of enhanced femtosecond supercontinuum via a weak seed in noble gases. Optics express, 24(13), 1511015119.
 Bahl, A., Brown, J., Wright, E., & Kolesik, M. (2015). Assessment of the metastable electronic state approach as a microscopically selfconsistent description for the nonlinear response of atoms. Optics letters, 40, 49874990.
 Brown, J. M., & Kolesik, M. (2015). Properties of Stark Resonant States in Exactly Solvable Systems. Advances in Mathematical Physics, 2015.
 Kolesik, M., Brown, J., Teleki, A., Jakobsen, P., Moloney, J., & Wright, E. (2015). Metastable electronic states and nonlinear response for highintensity optical pulses: erratum. Optica, 2, 509509.
 Mills, M., Heinrich, M., Kolesik, M., & Christodoulides, D. (2015). Extending optical filaments using auxiliary dress beams. Journal of Physics B: Atomic, Molecular and Optical Physics, 48, 094014.
 Panagiotopoulos, P., Whalen, P., Kolesik, M., & Moloney, J. V. (2015). Carrier field shock formation of longwavelength femtosecond pulses in singlecrystal diamond and air. JOSA B, 32, 17181730.
 Panagiotopoulos, P., Whalen, P., Kolesik, M., & Moloney, J. V. (2015). Super high power midinfrared femtosecond light bullet. Nature Photonics, 9, 543548.
 Andreasen, J., Bahl, A., & Kolesik, M. (2014). Spatial effects in supercontinuum generation in waveguides. Optics express, 22(21), 2575625767.
 Bahl, A., Teleki, A., Jakobsen, P. K., Wright, E. M., Kolesik, M., & others, . (2014). Reflectionless beam propagation on a piecewise linear complex domain. Journal of Lightwave Technology, 32(22), 36703676.
 Kolesik, M., & Moloney, J. (2014). Modeling and simulation techniques in extreme nonlinear optics of gaseous and condensed media. Reports on Progress in Physics, 77(1), 016401.
 Kolesik, M., Brown, J., Moloney, J., & Faccio, D. (2014). Historydependent effects in subcyclewaveform strongfield ionization. Physical Review A, 90(3), 033414.
 Kolesik, M., Brown, J., Teleki, A., Jakobsen, P., Moloney, J., & Wright, E. (2014). Metastable electronic states and nonlinear response for highintensity optical pulses. Optica, 1(5), 323331.
 Roger, T., Majus, D., Tamosauskas, G., Panagiotopoulos, P., Kolesik, M., Genty, G., \.e\fi, I., Dubietis, A., & Faccio, D. (2014). Extreme events in resonant radiation from threedimensional light bullets. Phys. Rev. A, 90, 033816.
 Scheller, M., Chen, X., Ariunbold, G. O., Born, N., Moloney, J., Kolesik, M., & Polynkin, P. (2014). Raman conversion in intense femtosecond Bessel beams in air. Physical Review A, 89(5), 053805.
 Scheller, M., Mills, M. S., Miri, M., Cheng, W., Moloney, J. V., Kolesik, M., Polynkin, P., & Christodoulides, D. N. (2014). Externally refuelled optical filaments. Nature Photonics, 8(4), 297301.
 Schuh, K., Kolesik, M., Wright, E., & Moloney, J. (2014). Simple model for the nonlinear optical response of gases in the transparency region. Optics letters, 39(17), 50865089.
 Whalen, P., Panagiotopoulos, P., Kolesik, M., & Moloney, J. V. (2014). Extreme carrier shocking of intense longwavelength pulses. Physical Review A, 89(2), 023850.
 Andreasen, J., & Kolesik, M. (2013). Coreconfined beam propagation method for guiding and leaky structures. Journal of Lightwave Technology, 31(18), 29993005.More infoAbstract: We propose a beam propagation method variant, which only requires the computational grid to span the core region of a waveguiding structure. Optical fields outside are accounted for through semitransparent boundary conditions. The method is fast, simple to implement, and accurate over a wide range of wavelengths. It is applicable to both guiding and leaky structures, with sufficiently strong contrast between refractive indices in the core, substrate, and cladding. © 19832012 IEEE.
 Andreasen, J., & Kolesik, M. (2013). Midinfrared femtosecond laser pulse filamentation in hollow waveguides: A comparison of simulation methods. Physical Review E  Statistical, Nonlinear, and Soft Matter Physics, 87(5).More infoAbstract: This work compares computational methods for laser pulse propagation in hollow waveguides filled with rare gases at high pressures, with applications in extreme nonlinear optics in the midinfrared wavelength region. As the wavelength of light λ=2π/k increases with respect to the transverse size R of a leaky waveguide, the loss of light out of the waveguide upon propagation, in general, increases. The now standard numerical approach for studying such structures is based on expansion of the propagating field into approximate leaky waveguide modes. We compare this approach to an improved method that resolves the electric field in real space and correctly captures the energy loss through the waveguide wall. The comparison reveals that the expansionbased approach overestimates losses that occur in nonlinearly reshaped pulsed waveforms. For a modest increase in computational effort, the alternate method offers a physically more accurate model to describe phenomena (e.g., extreme pulseselfcompression) in waveguides with smaller values of kR. © 2013 American Physical Society.
 Chen, X., Polynkin, P., & Kolesik, M. (2013). Raman effect in selffocusing of fewcycle laser pulses in air. Optics Letters, 38(12), 20172019.More infoPMID: 23938962;Abstract: Selffocusing of ultrashort pulses in air is investigated by means of numerical simulations. The role of the vibrational Raman effect and its dependence on pulse chirp is studied, with results shedding new light on the interpretation of the measurements of the critical selffocusing power. We also discuss computational modeling issues important specifically for fewcycle pulses. © 2013 Optical Society of America.
 Mills, M. S., Kolesik, M., & Christodoulides, D. N. (2013). Dressed optical filaments. Optics Letters, 38(1), 2527.More infoPMID: 23282826;Abstract: In this Letter we show that by appropriately providing an auxiliary "dress" beam one can extend the longevity of an optical filament by almost one order of magnitude. These optical dressed filaments can propagate substantially further by judiciously harnessing energy from their secondary beam reservoir. This possibility is theoretically investigated in air when the filament is dressed with a conically convergent annular Gaussian beam. © 2012 Optical Society of America.
 Mills, M. S., Scheller, M., Polynkin, P., Kolesik, M., & Christodoulides, D. N. (2013). Propagation dynamics of dressed optical filaments. CLEO: QELS_Fundamental Science, CLEO:QELS FS 2013, QW1E.2.More infoAbstract: We demonstrate that the longevity of an optical filament can be drastically increased by prudently surrounding it with a low intensity annular wavefront. Here, we systematically study the dynamics and robustness of such dressed filaments. © OSA 2013.
 Polynkin, P., & Kolesik, M. (2013). Critical power for selffocusing in the case of ultrashort laser pulses. Physical Review A  Atomic, Molecular, and Optical Physics, 87(5).More infoAbstract: We attempt to evaluate the applicability of the concept of the critical power for selffocusing, originally developed for intense quasicontinuouswave (cw) beams, to the case of ultraintense and ultrashort laser pulses propagating in air. Our results show that, unlike in the cw case, no particular value of peak pulse power can be viewed as a sharp demarcation line between linear and nonlinear propagation regimes. Our analysis further reveals the important role played by chromatic dispersion in the propagation dynamics of the laser pulse. © 2013 American Physical Society.
 Andreasen, J., & Kolesik, M. (2012). Nonlinear propagation of light in structured media: Generalized unidirectional pulse propagation equations. Physical Review E  Statistical, Nonlinear, and Soft Matter Physics, 86(3).More infoAbstract: Unidirectional pulse propagation equations [UPPE, Phys. Rev. E10.1103/PhysRevE.70.036604 70, 036604 (2004)] have provided a theoretical underpinning for computeraided investigations into dynamics of highpower ultrashort laser pulses and have been successfully utilized for almost a decade. Unfortunately, they are restricted to applications in bulk media or, with additional approximations, to simple waveguide geometries in which only a few guided modes can approximate the propagating waveform. The purpose of this work is to generalize the directional pulse propagation equations to structures characterized by strong refractive index differences and material interfaces. We also outline a numerical solution framework that draws on the combination of the bulkmedia UPPE method with singlefrequency beampropagation techniques. © 2012 American Physical Society.
 Brown, J. M., Lotti, A., Teleki, A., & Kolesik, M. (2011). Exactly solvable model for nonlinear lightmatter interaction in an arbitrary timedependent field. Physical Review A  Atomic, Molecular, and Optical Physics, 84(6).More infoAbstract: Exact analytic expressions are derived for the dipole moment and nonlinear current of a onedimensional quantum particle subject to a shortrange attractive potential and an arbitrary timedependent electric field. An efficient algorithm for the current evaluation is described and a robust implementation suitable for numerical simulations is demonstrated. © 2011 American Physical Society.
 Couairon, A., Brambilla, E., Corti, T., Majus, D., de, O., & Kolesik, M. (2011). Practitioner's guide to laser pulse propagation models and simulation. European Physical Journal: Special Topics, 199(1), 576.More infoAbstract: The purpose of this article is to provide practical introduction into numerical modeling of ultrashort optical pulses in extreme nonlinear regimes. The theoretic background section covers derivation of modern pulse propagation models starting from Maxwell's equations, and includes both envelopebased models and carrierresolving propagation equations. We then continue with a detailed description of implementation in software of Nonlinear Envelope Equations as an example of a mixed approach which combines finitedifference and spectral techniques. Fully spectral numerical solution methods for the Unidirectional Pulse Propagation Equation are discussed next. The modeling part of this guide concludes with a brief introduction into efficient implementations of nonlinear medium responses. Finally, we include several workedout simulation examples. These are miniprojects designed to highlight numerical and modeling issues, and to teach numericalexperiment practices. They are also meant to illustrate, first and foremost for a nonspecialist, how tools discussed in this guide can be applied in practical numerical modeling. © 2011 EDP Sciences and Springer.
 Rubino, E., Belgiorno, F., Cacciatori, S. L., Clerici, M., Gorini, V., Ortenzi, G., Rizzi, L., Sala, V. G., Kolesik, M., & Faccio, D. (2011). Experimental evidence of analogue Hawking radiation from ultrashort laser pulse filaments. New Journal of Physics, 13.More infoAbstract: Curved spacetimes and, in particular, event horizons of astrophysical black holes are expected to excite the quantum vacuum and give rise to an emission of quanta known as Hawking radiation. Remarkably, many physical systems may be considered analogous to black holes and as such hold promise for the detection of Hawking radiation. In particular, recent progress in the field of transformation optics, i.e. the description of optical systems in terms of curved spacetime geometries, has led to a detailed description of methods for generating, via superluminal dielectrics, a blocking horizon for photons. Our measurements highlight the emission of photons from a moving refractive index perturbation induced by a laser pulse that is in quantitative agreement with the Hawking model. This opens an intriguing and readily accessible observation window into quantum field theory in curved spacetime geometries. © IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.
 Zang, L., Kang, M. S., Kolesik, M., Scharrer, M., & Russell, P. (2010). Dispersion of photonic Bloch modes in periodically twisted birefringent media. Journal of the Optical Society of America B: Optical Physics, 27(9), 17421750.More infoAbstract: We investigate the polarization evolution and dispersive properties of the eigenmodes of birefringent media with arbitrarily twisted axes of birefringence. Analytical and numerical methods based on a transfer matrix approach are developed and used to study specifically helically twisted structures and the Bloch modes of periodically twisted media, as represented in particular by structural "rocking" filters inscribed in highly birefringent photonic crystal fibers. The presence of periodically twisted birefringence axes causes the group velocity dispersion curves to separate strongly from each other in the vicinity of the anticrossing wavelength, where the interpolarization beatlength equals an integer multiple of the rocking period. The maximum separation between these curves and the bandwidth of the splitting depend on the amplitude of the rocking angle. We also show that suitably designed adiabatic transitions, formed by chirping the rocking period, allow a broadband conversion between a linearly polarized fiber eigenmode and a single Bloch mode of a uniform rocking filter. The widely controllable dispersive properties provided by rocking filters may be useful for manipulating the phasematching conditions in nonlinear optical processes such as fourwave mixing, supercontinuum generation, and the generation of resonant radiation from solitons. © 2010 Optical Society of America.
 Clerici, M., Jedrkiewicz, O., Faccio, D., Averchi, A., Lotti, A., Rubino, E., Valiulis, G., Tartara, L., Degiorgio, V., Couairon, A., Kolesik, M., & Trapani, P. D. (2008). Controlled X wave formation in bulk quadratic and cubic nonlinear media. Conference Proceedings  Lasers and ElectroOptics Society Annual MeetingLEOS, 874875.
 Buendía, G., Rikvold, P. A., Kolesik, M., Park, K., & Novotny, M. A. (2007). Nanostructure and velocity of fielddriven solidonsolid interfaces moving under a phononassisted dynamic. Physical Review B  Condensed Matter and Materials Physics, 76(4).More infoAbstract: The nanoscopic structure and the stationary propagation velocity of (1+1) dimensional solidonsolid interfaces in an Ising latticegas model, which are driven far from equilibrium by an applied force, such as a magnetic field or a difference in (electro)chemical potential, are studied by an analytic nonlinearresponse approximation together with kinetic Monte Carlo simulations. Here, we consider the case that the system is coupled to a twodimensional phonon bath. In the resulting dynamic, transitions that conserve the system energy are forbidden, and the effects of the applied force and the interaction energies do not factorize (a socalled hard dynamic). In full agreement with previous general theoretical results, we find that the local interface width changes dramatically with the applied force. However, in contrast with other hard dynamics, this change is nonmonotonic in the driving force. Results are also obtained for the force dependence and anisotropy of the interface velocity, which also show differences in good agreement with the theoretical expectations for the differences between soft and hard dynamics. However, significant differences between theory and simulation are found near two special values of the driving force, where certain transitions allowed by the solidonsolid model become forbidden by the phononassisted dynamic. Our results show that different stochastic interface dynamics that all obey detailed balance and the same conservation laws nevertheless can lead to radically different interface responses to an applied force. Thus, they represent a significant step toward providing a solid physical foundation for kinetic Monte Carlo simulations. © 2007 The American Physical Society.
 Faccio, D., Averchi, A., Kolesik, M., Couairon, A., Polesana, P., Tamosauskas, G., Dubietis, A., Trapani, P. D., & Piskarskas, A. (2007). Pulse compression and X wave generation by CrossPhaseModulation induced spatiotemporal reshaping. Conference on Lasers and ElectroOptics Europe  Technical Digest.
 Buendía, G., Rikvold, P. A., & Kolesik, M. (2006). Fielddriven solidonsolid interfaces moving under a stochastic Arrhenius dynamics: Effects of the barrier height. Journal of Molecular Structure: THEOCHEM, 769(13), 207210.More infoAbstract: We present analytical results and kinetic Monte Carlo simulations for the mobility and microscopic structure of solidonsolid (SOS) interfaces driven far from equilibrium by an external force, such as an applied field or (electro)chemical potential difference. The interfaces evolve under a specific stochastic dynamics with a local energy barrier (an Arrhenius dynamics), known as the transition dynamics approximation (TDA). We calculate the average height of steps on the interface, the average interface velocity, and the skewness of the interface as functions of the driving force and the height of the energy barrier. We find that the microscopic interface structure depends quite strongly on the barrier height. As the barrier becomes higher, the local interface width decreases and the skewness increases, suggesting increasing shortrange correlations between the step heights. © 2006 Elsevier B.V. All rights reserved.
 Buendía, G., Rikvold, P. A., & Kolesik, M. (2006). Microstructure and velocity of fielddriven solidonsolid interfaces moving under stochastic dynamics with local energy barriers. Physical Review B  Condensed Matter and Materials Physics, 73(4).More infoAbstract: We study the microscopic structure and the stationary propagation velocity of (1+1) dimensional solidonsolid interfaces in an Ising latticegas model, which are driven far from equilibrium by an applied force, such as a magnetic field or a difference in (electro)chemical potential. We use an analytic nonlinearresponse approximation [P. A. Rikvold and M. Kolesik, J. Stat. Phys. 100, 377 (2000)] together with kinetic Monte Carlo simulations. Here we consider interfaces that move under Arrhenius dynamics, which include a microscopic energy barrier between the allowed Ising or latticegas states. Two different dynamics are studied: the standard onestep dynamics (OSD) [H. C. Kang and W. Weinberg, J. Chem. Phys. 90, 2824 (1992)] and the twostep transitiondynamics approximation (TDA) [T. AlaNissila, J. Kjoll, and S. C. Ying, Phys. Rev. B 46, 846 (1992)]. In the OSD the effects of the applied force and the interaction energies in the model factorize in the transition rates (soft dynamics), while in the TDA such factorization is not possible (hard dynamics). In full agreement with previous general theoretical results we find that the local interface width under the TDA increases dramatically with the applied force. In contrast, the interface structure with the OSD is only weakly influenced by the force, in qualitative agreement with the theoretical expectations. Results are also obtained for the force dependence and anisotropy of the interface velocity, which also show differences in good agreement with the theoretical expectations for the differences between soft and hard dynamics. Our results confirm that different stochastic interface dynamics that all obey detailed balance and the same conservation laws nevertheless can lead to radically different interface responses to an applied force. © 2006 The American Physical Society.
 Teipel, J., Türke, D., Giessen, H., Killi, A., Morgner, U., Lederer, M., Kopf, D., & Kolesik, M. (2005). Diodepumped, ultrafast, multioctave supercontinuum source at repetition rates between 500 kHz and 20 MHz using Yb:glass lasers and tapered fibers. Optics Express, 13(5), 14771485.More infoPMID: 19495023;Abstract: We present a compact, all diodepumped supercontinuum source based on a SESAM modelocked Yb:glass oscillator at 1040 nm and a tapered fiber. The oscillator has a repetition rate of 20 MHz, a pulse duration of 200 fs, and a maximum pulse energy of about 15 nJ. This system delivers an 1100 nm broad spectrum with an output power of more than 100 mW. Decreasing the repetition rate to 500 kHz by cavitydumping results in a supercontinuum with a high pulse energy of about 50 nJ. Furthermore, using the frequencydoubled output of this laser at 520 nm with 300 fs pulse duration resulted in supercontinua in the nearUV and visible spectral region. We compare the experimental spectra with theoretical simulations. © 2005 Optical Society of America.
 Kolesik, M., Novotny, M. A., & Rikvold, P. A. (2003). Extreme longtime dynamic Monte Carlo simulations for metastable decay in the d = 3 Ising ferromagnet. International Journal of Modern Physics C, 14(1), 121131.More infoAbstract: We study the extreme longtime behavior of the metastable phase of the threedimensional Ising model with Glauber dynamics in an applied magnetic field and at a temperature below the critical temperature. For these simulations, we use the advanced simulation method of projective dynamics. The algorithm is described in detail, together with its application to the escape from the metastable state. Our results for the field dependence of the metastable lifetime are in good agreement with theoretical expectations and span more than 50 decades in time.
 Rikvold, P. A., & Kolesik, M. (2003). Microstructure and velocity of fielddriven Ising interfaces moving under a soft stochastic dynamic. Physical Review E  Statistical, Nonlinear, and Soft Matter Physics, 67(6 2), 066113/1066113/8.More infoAbstract: A study was conducted on the dependence of the local structure of driven interfaces on the applied field and temperature and on the form of the stochastic dynamics under which they move. The local interface structure was of interest because it was this, rather than the largescale scaling behavior, which determines such important interface properties as mobility and chemical reactivity. It was found that the results for the Ising model with the soft Glauber dynamic differ relatively little from the solidonsolid (SOS) model with soft dynamics studied.
 Rikvold, P. A., & Kolesik, M. (2002). Microstructure and velocity of fielddriven solidonsolid interfaces: Analytic approximations and numerical results. Physical Review E  Statistical, Nonlinear, and Soft Matter Physics, 66(6), 066116/1066116/12.More infoPMID: 12513356;Abstract: The local structure of a solidonsolid interface in a twodimensional kinetic Ising ferromagnet or attractive latticegas model with singlespinflip Glauber dynamics, which is driven far from equilibrium by an applied field or chemical potential, is studied by an analytic meanfield, nonlinearresponse theory [P. A. Rikvold and M. Kolesik, J. Stat. Phys. 100, 377 (2000)], and by dynamic Monte Carlo simulations. The probability density of the height of an individual step in the surface is obtained, both analytically and by simulation. The width of the probability density is found to increase dramatically with the magnitude of the applied field, with close agreement between the theoretical predictions and the simulation results. Excellent agreement between theory and simulations is also found for the field dependence and anisotropy of the interface velocity. The joint distribution of nearestneighbor step heights is obtained by simulation. It shows increasing correlations with increasing field, similar to the skewness observed in other examples of growing surfaces. © 2002 The American Physical Society.
 Rikvold, P. A., & Kolesik, M. (2002). Soft versus hard dynamics for fielddriven solidonsolid interfaces. Journal of Physics A: Mathematical and General, 35(9), L117L123.More infoAbstract: Analytical arguments and dynamic Monte Carlo simulations show that the microscopic structure of fielddriven solidonsolid interfaces depends strongly on the details of the dynamics. For nonconservative dynamics with transition rates that factorize into parts dependent only on the changes in interaction energy and field energy, respectively (soft dynamics), the intrinsic interface width is field independent. For nonfactorizing rates, such as the standard Glauber and Metropolis algorithms (hard dynamics), it increases with the field. Consequences for the interface velocity and its anisotropy are discussed.
 Novotny, M. A., Rikvold, P. A., Kolesik, M., Townsley, D. M., & Ramos, R. A. (2000). Simulations of metastable decay in two and threedimensional models with microscopic dynamics. Journal of NonCrystalline Solids, 274(1), 356363.More infoAbstract: We present a brief analysis of the crossover phase diagram for the decay of a metastable phase in a simple dynamic latticegas model of a twophase system. We illustrate the nucleationtheoretical analysis with dynamic Monte Carlo (MC) simulations of a kinetic Ising lattice gas on square and cubic lattices. We predict several regimes in which the metastable lifetime has different functional forms, and we provide estimates for the crossovers between the different regimes. In the multidroplet regime, the KolmogorovJohnsonMehlAvrami (KJMA) theory for the time dependence of the orderparameter decay and the twopoint density correlation function allows extraction of both the order parameter in the metastable phase and the interfacial velocity from the simulation data.
 Rikvold, P. A., & Kolesik, M. (2000). Analytic approximations for the velocity of fielddriven Ising interfaces. Journal of Statistical Physics, 100(12), 377403.More infoAbstract: We present analytic approximations for the field, temperature, and orientation dependences of the interface velocity in a twodimensional kinetic Ising model in a nonzero field. The model, which has nonconserved order parameter, is useful for ferromagnets, ferroelectrics, and other systems undergoing orderdisorder phase transformations driven by a bulk freeenergy difference. The solidonsolid (SOS) approximation for the microscopic surface structure is used to estimate mean spinclass populations, from which the mean interface velocity can be obtained for any specific singlespinflip dynamic. This linearresponse approximation remains accurate for higher temperatures than the singlestep and polynuclear growth models, while it reduces to these in the appropriate lowtemperature limits. The equilibrium SOS approximation is generalized by meanfield arguments to obtain fielddependent spinclass populations for moving interfaces, and thereby a nonlinearresponse approximation for the velocity. The analytic results for the interface velocity and the spinclass populations are compared with Monte Carlo simulations. Excellent agreement is found in a wide range of field, temperature, and interface orientation.
 Novotny, M. A., Kolesik, M., & Rikvold, P. A. (1999). Slow forcing in the projective dynamics method. Computer Physics Communications, 121, 330333.More infoAbstract: We provide a proof that when there is no forcing the recently introduced projective dynamics Monte Carlo algorithm gives the exact lifetime of the metastable state, within statistical uncertainties. We also show numerical evidence illustrating that for slow forcing the approach to the zeroforcing limit is rather rapid. The model studied numerically is the 3dimensional 3state Potts ferromagnet.
 Kolesik, M., Novotny, M. A., & Rikvold, P. A. (1998). Projection method for statics and dynamics of lattice spin systems. Physical Review Letters, 80(15), 33843387.More infoAbstract: A method based on Monte Carlo sampling of probability flows projected onto the subspace of one or more slow variables, is proposed for the investigation of dynamic and static properties of lattice spin systems. We illustrate the method by applying it, with projection onto the orderparameter subspace, to the 3dimensional 3state Potts model in equilibrium, and to metastable decay in a 3dimensional 3state kinetic Potts model.
 Novotny, M. A., Kolesik, M., & Rikvold, P. A. (1998). Magnetization switching in singledomain ferromagnets. Journal of Magnetism and Magnetic Materials, 177181(PART 2), 917918.More infoAbstract: A model for singledomain uniaxial ferromagnetic particles with high anisotropy, the Ising model, is studied. Recent experimental observations have been made of the probability that the magnetization has not switched, PNOT. Here an approach is described in which it is emphasized that a ferromagnetic particle in an unfavorable field is, in fact, a metastable system, and the switching is accomplished through the nucleation and subsequent growth of localized droplets. Nucleation theory is applied to finite systems to determine the coercivity as a function of particle size and to calculate PNOT. Both of these quantities are modified by different boundary conditions, magnetostatic interactions, and quenched disorder. © 1998 Elsevier Science B.V. All rights reserved.
 Kolesik, M., Novotny, M. A., & Rikvold, P. A. (1997). Monte Carlo simulation of magnetization reversal in Fe sesquilayers on W(110). Physical Review B  Condensed Matter and Materials Physics, 56(18), 1179111796.More infoAbstract: Iron sesquilayers grown at room temperature on W(110) exhibit a pronounced coercivity maximum near a coverage of 1.5 atomic monolayers. On lattices which faithfully reproduce the morphology of the real films, a kinetic Ising model is utilized to simulate the domainwall motion. Simulations reveal that the dynamics is dominated by the secondlayer islands, which act as pinning centers. The simulated dependences of the coercivity on the film coverage, as well as on the temperature and the frequency of the applied field, are very similar to those measured in experiments. Unlike previous micromagnetic models, the presented approach provides insight into the dynamics of the domainwall motion and clearly reveals the role of thermal fluctuations.
 Kolesik, M., Novotny, M. A., & Rikvold, P. A. (1997). Monte Carlo simulation of magnetization reversal via domainwall motion in Fe sesquilayers on W(110). Materials Research Society Symposium  Proceedings, 492, 313318.More infoAbstract: Iron sesquilayers are ultrathin films with coverages between one and two atomic monolayers. They consist of an almost defectfree monolayer with compact islands of a second atomic layer on top. This variation of the film thickness results in a strong interaction between domain walls and the island structure. It makes these systems an ideal laboratory to study the dynamics of domain walls driven by weak external fields. We present computer simulations which provide insight into the role of the thermally activated nucleation processes by which a driven domain wall overcomes the obstacles created by the islands.
 Kolesik, M., Richards, H. L., Novotny, M. A., Rikvold, P. A., & Lindgård, P. (1997). Magnetization switching in nanoscale ferromagnetic grains: Simulations with heterogeneous nucleation. Journal of Applied Physics, 81(8 PART 2B), 56005602.More infoAbstract: We present results obtained with various types of heterogeneous nucleation in a kinetic Ising model of magnetization switching in singledomain ferromagnetic nanoparticles. We investigate the effect of the presence of the system boundary and make comparison with simulations on periodic lattices. We also study systems with bulk disorder and compare how two different types of disorder influence the switching behavior. © 1997 American Institute of Physics.
 Richards, H. L., Kolesik, M., Lindgård, P., Rikvold, P. A., & Novotny, M. A. (1997). Effects of boundary conditions on magnetization switching in kinetic Ising models of nanoscale ferromagnets. Physical Review B  Condensed Matter and Materials Physics, 55(17), 1152111540.More infoAbstract: Magnetization switching in highly anisotropic singledomain ferromagnets has been previously shown to be qualitatively described by the droplet theory of metastable decay and simulations of twodimensional kinetic Ising systems with periodic boundary conditions. In this paper we consider the effects of boundary conditions on the switching phenomena, A rich range of behaviors is predicted by droplet theory: the specific mechanism by which switching occurs depends on the structure of the boundary, the particle size, the temperature, and the strength of the applied field. The theory predicts the existence of a peak in the switching field as a function of system size in both systems with periodic boundary conditions and in systems with boundaries. The size of the peak is strongly dependent on the boundary effects. It is generally reduced by open boundary conditions, and in some cases it disappears if the boundaries are too favorable towards nucleation. However, we also demonstrate conditions under which the peak remains discernible. This peak arises as a purely dynamic effect and is not related to the possible existence of multiple domains. We illustrate the predictions of droplet theory by Monte Carlo simulations of twodimensional Ising systems with various system shapes and boundary conditions.
 Kolesik, M., & Suzuki, M. (1995). Accurate estimates of 3D Ising critical exponents using the coherentanomaly method. Physica A: Statistical Mechanics and its Applications, 215(12), 138151.More infoAbstract: An analysis of the critical behavior of the threedimensional Ising model using the coherentanomaly method (CAM) is presented. Various sources of errors in CAM estimates of critical exponents are discussed, and an improved scheme for the CAM data analysis is tested. Using a set of meanfield type approximations based on the variational series expansion approach, accuracy comparable to the most precise conventional methods has been achieved. Our results for the critical exponents are given by α = 0.108(5), β = 0.327(4), γ = 1.237(4) and δ = 4.77(5). © 1995.
 Kolesik, M., & Suzuki, M. (1995). Critical exponents of the 3D antiferromagnetic threestate Potts model using the coherentanomaly method. Physica A: Statistical Mechanics and its Applications, 216(4), 469477.More infoAbstract: The antiferromagnetic threestate Potts model on the simplecubic lattice is studied using the coherentanomaly method (CAM). The CAM analysis provides the estimates for the critical exponents which indicate the XY universality class, namely α = 0.011, β = 0.351, γ = 1.309 and δ = 4.73. This observation corroborates the results of the recent Monte Carlo simulations, and disagrees with the proposal of a new universality class. © 1995.
 Kolesik, M., & Suzuki, M. (1995). On the lowtemperature ordering of the 3D antiferromagnetic threestate Potts model. Journal of Physics A: General Physics, 28(23), 65436555.More infoAbstract: The antiferromagnetic threestate Potts model on the simplecubic lattice is studied using Monte Carlo simulations. The ordering in a medium temperature range below the critical point is investigated in detail. Two different regimes have been observed: the socalled broken sublatticesymmetry phase dominates at sufficiently low temperatures, while the phase just below the critical point is characterized by an effectively continuous order parameter and by a fully restored rotational symmetry. However, the latter phase is not the permutationally sublattice symmetric phase recently predicted by the cluster variation method.
 Aamaj, L., Percus, J. K., & Kolesík, M. (1994). Twodimensional onecomponent plasma at coupling Γ=4: Numerical study of pair correlations. Physical Review E, 49(6), 56235627.More infoAbstract: We consider a classical twodimensional onecomponent plasma of charged particles in a circularly symmetric neutralizing background, at a coupling constant of Γ=e2/kBT=4. The numerical results, based technically on a successive increase of the number of particles and on a Van der Monde determinantal representation of Boltzmann factors, strongly indicate a Gaussiantype falloff of the truncated bulk chargecharge correlations, similarly as in the exactly solvable Γ=2 case. © 1994 The American Physical Society.
 Kolesík, M. (1994). Nofreeends method for lattice animals and vertex models with arbitrary number of states. Physica A: Statistical Mechanics and its Applications, 202(34), 529539.More infoAbstract: A new nofreeends method for generating power series expansions for vertex models with arbitrary number of states is described. It is based on firstorder differential equations providing a recurrent relation connecting the freeends part of the series expansion in the (n + 1)th order with the nth and (n  1)th order of the complete expansion. The number of necessary nofreeends graphs can be further reduced for symmetric systems by eliminating the nodes which have two neighbors and their incident bonds are in different states. © 1994.
 Kolesík, M., & Šamaj, L. (1993). Evidence for the nonuniversality of a 3D vertex model. Physics Letters A, 177(1), 8792.More infoAbstract: By combining a variational series expansion and the coherent anomaly method we observe a strong evidence for the nonuniversal critical behaviour of the symmetric vertex model formulated on the 2D square and 3D diamond lattices. © 1993.
 Kolesík, M., & Šamaj, L. (1993). Series expansion and CAM study of the nonuniversal behavior of the symmetric 16vertex model. Journal of Statistical Physics, 72(56), 12031226.More infoAbstract: The phase diagram and the critical indices are investigated for the symmetric 16vertex model on the square lattice by combining a variational series expansion and the coherent anomaly method. The nonuniversal critical exponents smoothly interpolate between two exactly solvable cases, namely the Baxter eightvertex model and the Ising model. © 1993 Plenum Publishing Corporation.
 Šamaj, L., & Kolesík, M. (1993). Selfduality of the O(2) gauge transformation and the phase structure of vertex models. Physica A: Statistical Mechanics and its Applications, 193(1), 157168.More infoAbstract: For the symmetric twostate vertex model formulated on a lattice with an arbitrary coordination number q, we construct a variational series expansion of the free energy with a free gauge parameter playing the role of the variational variable. In the lowest order of the variational series expansion we obtain the Bethe approximation. Its analytical treatment provides a new method of searching for the selfdual manifolds for lattices of higher coordination number q and gives some information about the internal structure of the selfdual manifolds where the first and secondorder phase transitions take place. The results are systematically improved by considering higherorder terms in the variational series expansion. © 1993.
 Kolesík, M., Tunega, D., & Sobolev, B. P. (1992). Lowfrequency vibrational modes in fluorite based superionic conductors  Raman study. Solid State Ionics, 58(34), 237242.More infoAbstract: The lowfrequency Raman spectra of various fluoritestructured superionic conductors are presented. The observed lowfrequency peaks are interpreted in terms of resonant vibrational states which are (quasi) localized in the vicinity of dopant clusters. © 1992.
 Šamaj, L., & Kolesík, M. (1992). Mapping of the symmetric vertex model onto the Ising model for an arbitrary lattice coordination. Physica A: Statistical Mechanics and its Applications, 182(3), 455466.More infoAbstract: We investigate the mapping of the twostate symmetric vertex model onto the Ising model in a field for a lattice with general coordination number q. The analysis is based on a combination of decoration and generalized weakgraph transformations. It is shown that the mapping is restricted to the manifold formed by the intersection of q3 hypersurfaces in the space of symmetric vertex weights. We also derive dual relations among the model parameters and the free energies of the related symmetric vertex and Ising models. The choices of the energies of the vertex configurations for which the mapping can be performed throughout the entire temperature range are briefly discussed. © 1992.
 Kolesik, M., & Sobolev, B. P. (1991). Raman scattering study of the Ba_{1x}R_{x}F_{2+x} superionic conductors. Solid State Ionics, 47(34), 325329.More infoAbstract: The nonstoichiometric solid solutions Ba1xRxF2+x (where R = Ce, Pr, Gd, Tb; and 0.10≤x≤0.30) are studied by Raman spectroscopy. Besides the main Raman peak, two additional bands are observed on its highfrequency side. A new interpretation for the extrapeaks, based on the cluster structure, is proposed. The observed vibrations seem to be connected with 8:12: (0, 1, 2) clusters. © 1991.
 Kolesík, M., & Šamaj, L. (1991). On correlation functions of twostate vertex models on the honeycomb lattice. Physica A: Statistical Mechanics and its Applications, 179(1), 145157.More infoAbstract: A family of twostate vertex models on the honeycomb lattice is solved exactly using a generalized weakgraph transformation technique. Two systems are analyzed in detail: the "spinflip" symmetric model, which is integrable in the whole temperature range, and the symmetric model, solvable at one specific temperature β*. This temperature turns out to be significant from the point of view of edgeedge correlations, namely they vanish at β*. © 1991.
 Kolesík, P. V., & Kolesík, M. (1989). Injuriousness of lentil gall midge Contarinia lentis Aczél (Diptera, Cecidomyiidae) and its distribution in Czechoslovakia. Anzeiger für Schädlingskunde Pflanzenschutz Umweltschutz, 62(8), 150156.More infoAbstract: There has been an outbreak of lentil gall midge C. lentis on lentil in Czechoslovakia since 1984. Its extent and cause are given below. Infestation symptomthe gallis described. The level of damage is defined as percent rate of the amount of galls and flower buds containing eggs of C. lentis from the amount of pods, flowers, galls, flower buds containing eggs and flower buds not containing eggs. As estimation is taken mostly at the end of lentil vegetation the level of damage is simplified to percent rate of the amount of galls from the amount of galls and pods. The level of damage of small parts of lentil stands can achieve 80%. Single plants with 95% level of damage were found. Proportional dependence of thousandkernel weight on the level of damage was found out. Two functions of dependence of yield on the level of damage were compared: a simple linear one and a quadratic one which included the dependence of thousandkernel weight on the level of damage. The generally valid linear function of dependence of yield y on the level of damage x is recommended: y=0.01 e (100x) where e means yield at 0% damage level. Recent distribution area of this species in Czechoslovakia is given. In 19851988 the average level of damage of lentil stands from outbreak centre ranged from 3.08 to 52.05%. The average yield losses in Czechoslovakia were 5.49% in 1985, 4.98% in 1986, 2.87% in 1987 and 3.97% in 1988. Concentration of lentil fields in a small area was recognized as the main cause of the outbreak of C. lentis in Czechoslovakia. © 1989 Verlag Paul Parey.
Proceedings Publications
 Hastings, M. G., Panagiotopoulos, P., Kolesik, M., & Moloney, J. V. (2023). The existence of multioctave spanning conical emission from ultrafast LWIR pulse filamentation. In Nonlinear Frequency Generation and Conversion: Materials and Devices XXII, 12405.
 Vasilyev, S., Vodopyanov, K., Kolesik, M., & Mirov, S. (2023). Using Crossphase Modulation to Transfer Coherence Between a Cr: Zns Frequency Comb and its Optical Pump. In 2023 Conference on Lasers and ElectroOptics (CLEO).
 Railing, L., Le, M., Salehi, F., Gu, J., Kolesik, M., & Milchberg, H. M. (2022). Time Resolved Measurements of the Polarization State of Supercontinuum Generated in a Monatomic Gas. In Frontiers in Optics.
 Rudenko, A., Hagen, M. K., Hader, J., Kolesik, M., Koch, S. W., & Moloney, J. V. (2022). Maxwellsemiconductor Bloch simulations of highharmonic generation in finite thickness semiconductor slabs. In Ultrafast Phenomena and Nanophotonics XXVI, 11999.
 Gu, J., Vasilyev, S., Mirov, M., & Kolesik, M. (2021). Modeling Harmonic and Supercontinuum Generation in Polycrystalline Materials. In CLEO: Science and Innovations.
 Hastings, M. G., Panagiotopoulos, P., Hasson, V., Kolesik, M., & Moloney, J. V. (2021). Filamentationfree Selfcompression of LWIR Pulses in a CO2 Gasfilled Multipass Cell: A Numerical Study. In CLEO: Applications and Technology.
 Heinz, J., & Kolesik, M. (2021). Transparent Boundary Conditions for QuasiBound States in Open Systems. In Frontiers in Optics.
 Panagiotopoulos, P., Kolesik, M., Hasson, V., Tochitsky, S., & Moloney, J. V. (2021). Simulations on the propagation dynamics of TW squareaperture CO2 laser pulses in the atmosphere. In CLEO: Science and Innovations.
 Valenzuela, A., Schweinsberg, A., Gu, J., Kolesik, M., Ensley, T., Vanderhoef, L., Tripepi, M., Wolfe, C., & Chowdhury, E. (2021). Spatial structure and simulations of midwave infrared ultrashort pulse laser frequency conversion in polycrystalline and amorphous optical materials. In Nonlinear Frequency Generation and Conversion: Materials and Devices XX, 11670.
 Vasilyev, S., Sander, M. Y., Gu, J., Smolski, V., Moskalev, I., Mirov, M., Barnakov, Y., Peppers, J., Kolesik, M., Mirov, S., & others, . (2021). Vector Solitons in a Kerrlens Modelocked Laser Oscillator. In Advanced Solid State Lasers.
 Schweinsberg, A., Tripepi, M., Talisa, N., Vanderhoef, L., Wolfe, C., Hastings, M., Kolesik, M., Valenzuela, A., & Chowdhury, E. (2020). Symmetry Properties of Harmonics Generated in SingleCrystal Chalcogenides Using Ultrafast MidInfrared Pulses. In 2020 Conference on Lasers and ElectroOptics (CLEO).
 Tripepi, M., Hastings, M., Schweinsberg, A., Vanderhoef, L., Wolfe, C., Ensley, T., Kolesik, M., Moloney, J., Chowdhury, E., & Valenzuela, A. (2020). Midwave infrared ultrashort pulse laser frequency conversion in single crystal, polycrystalline, and amorphous optical materials. In Nonlinear Frequency Generation and Conversion: Materials and Devices XIX, 11264.
 Ensley, T. R., Hastings, M. G., Werner, K., Schweinsberg, A., Tripepi, M., Talisa, N., Wilmer, B. L., Austin, D., Wolfe, C. M., Kolesik, M., & others, . (2019). Ultrafast Nonlinear Refraction Measurements of Transparent Materials in the MidInfrared for Modeling Harmonic and Supercontinuum Generation. In Nonlinear Optics.
 Hastings, M. G., Gu, J., & Kolesik, M. (2019). Simulation of harmonic and supercontinuum generation in polycrystalline media. In 2019 IEEE Research and Applications of Photonics in Defense Conference (RAPID).
 Hastings, M. G., Moloney, J. V., Kolesik, M., Werner, K., Austin, D., Chowdhury, E., Schweinsberg, A., Wilmer, B. L., Wolfe, C. M., Vanderhoef, L., & others, . (2019). Modeling Harmonic Generation in Poly crystalline ZnSe. In 2019 Conference on Lasers and ElectroOptics (CLEO).
 Panagiotopoulos, P., Rosenow, P., Tochitsky, S., Koch, S. W., Wright, E. W., Kolesik, M., & Moloney, J. V. (2019). 10 micron filaments in the atmosphere: Modeling perspective. In 2019 IEEE Research and Applications of Photonics in Defense Conference (RAPID).
 Tochitsky, S. Y., Welch, E. C., Polyanskiy, M. N., Pogorelsky, I. V., Panagiotopoulos, P., Kolesik, M., Wright, E. M., Koch, S. W., Moloney, J. V., & Joshi, C. (2019). Longwave Infrared Megafilaments in Air. In Nonlinear Optics.
 Werner, K., Tripepi, M., Schweinsberg, A., Hastings, M., Vanderhoef, L., Wolfe, C., Ensley, T., Wilmer, B., Kolesik, M., Moloney, J., & others, . (2019). Generation of broad spectral components from midwave infrared ultrashort pulse laser propagation through ZnSe and ZnS. In Nonlinear Frequency Generation and Conversion: Materials and Devices XVIII, 10902.
 Eftekhar, M. A., SanjabiEznaveh, Z., AntonioLopez, J. E., Aviles, H. L., Benis, S., Kolesik, M., Sch\"ulzgen, A., Wise, F. W., Correa, R. A., & Christodoulides, D. N. (2018). Accelerating nonlinear interactions in tapered multimode fibers. In CLEO: QELS\_Fundamental Science.
 Eznaveh, Z. S., Eftekhar, M. A., Lopez, J. A., Kolesik, M., Aviles, H. L., Wise, F. W., Christodoulides, D. N., & Correa, R. A. (2018). Blueenhanced supercontinuum generation in a gradedindex fluorinedoped multimode fiber. In 2018 Optical Fiber Communications Conference and Exposition (OFC).
 Panagiotopoulos, P., Rosenow, P., Schuh, K., Kolesik, M., Wright, E. M., Koch, S. W., & Moloney, J. V. (2018). Simulation of LWIR TW ultrashort pulses over kilometer ranges in the atmosphere. In Ultrafast Bandgap Photonics III, 10638.
 Bahl, A., Kolesik, M., Wahlstrand, J., & Milchberg, H. (2017). Theoryexperiment comparison of a quantum based lightmatter interaction model for optical filamentation. In Frontiers in Optics.
 Bahl, A., Wahlstrand, J. K., Milchberg, H. M., & Kolesik, M. (2017). Experimenttheory comparison and verification of metastable electronic state description of nonlinear optical response in atoms and molecules. In Nonlinear Optics.
 Moloney, J. V., Schuh, K., Panagiotopoulos, P., Kolesik, M., & Koch, S. W. (2017). Long range robust multiterawatt MWIR and LWIR atmospheric light bullets. In Ultrafast Bandgap Photonics II, 10193.
 Rafailov, M. K., Freericks, J. K., Abdurazakov, O., Kemper, A. F., Yabana, K., Hagan, D. J., Reichert, M. C., Zhao, P., Pattanaik, H. S., Van, S., & others, . (2017). Ultrafast Bandgap Photonics II. In Proc. of SPIE Vol, 10193.
 Eftekhar, M. A., Mills, M. S., Wright, L. G., Kolesik, M., Correa, R. A., Wise, F. W., & Christodoulides, D. N. (2016). Versatile supercontinuum generation in parabolic multimode optical fibers. In Lasers and ElectroOptics (CLEO), 2016 Conference on.
 Kolesik, M., Brown, J., & Bahl, A. (2016). Modeling of ultrafast laser pulse propagation. In SPIE Defense+ Security.
 Mansuripur, M., Kolesik, M., & Jakobsen, P. (2016). Leaky modes of dielectric cavities. In SPIE Nanoscience+ Engineering.
 Parto, M., Eftekhar, M. A., Miri, M., AmezcuaCorrea, R., Li, G., & Christodoulides, D. N. (2016). ZeroDGD multicore optical fibers. In CLEO: Applications and Technology.
 Faccio, D., Roger, T., Petev, M., Clerici, M., Morandotti, R., Legare, F., Majus, D., Tamosauskas, G., Dubietis, A., Couairon, A., & others, . (2014). Resonant Radiation from Collapsing Light Pulses and Spatiotemporal Light Bullets. In Laser Science, LTh3H1.
 Faccio, D., Roger, T., Petev, M., Clerici, M., Morandotti, R., Legare, F., Majus, D., Tamosauskas, G., Dubietis, A., Couairon, A., & others, . (2014). Spatiotemporal Light Bullets in Bulk Media. In Nonlinear Photonics, NM3A1.
 Kolesik, M. (2014). Nonperturbative timedomain modeling of lightmatter interactions for computer simulation in extreme nonlinear optics. In Laser Science, LTh3H4.
 Roger, T., Majus, D., Tamo\vsauskas, G., Dubietis, A., Genty, G., Kolesik, M., & Faccio, D. (2014). Extreme Events in Resonant Radiation from 3dimensional light bullets. In CLEO: QELS\_Fundamental Science, FW3D4.
Presentations
 Mansuripur, M., Kolesik, M., & Jakobsen, P. (2016, August). Spin and orbital angular momenta of trapped electromagnetic fields in leaky optical cavities. SPIE Optics & Photonics Symposium. San Diego, California.
 Kolesik, M. (2015, July). Firstprinciple Models of Lightmatter Interactions for Largescale Computer Simulations in Extreme Nonlinear Optics. 2015 Progress in Electromagnetic Research Symposium on Photonics. Prague: PIERS.
 Kolesik, M. (2015, October). Lightmatter interactions and optical filamentation.. NonLinear Optics contractor meeting. Arlington.
 Kolesik, M. (2015, October). Mathematical Properties and Physical Interpretation of Quantum Metastable States. Nonlinear Optics contractor meeting. Arlington.
 Kolesik, M. (2014, 03/10/14  03/14/14). Invited: Nonperturbative lightmatter interactions in extreme nonlinear optics. Workshop on Mathematical Methods and Models in Laser Filamentation. Montreal, Canada.
 Kolesik, M. (2014, 05/17/14  05/22/14). Invited: Pulse propagation and nonperturbative lightmatter interaction modeling for extreme nonlinear optics. The Nonlinear Meeting 201. HeriotWatt University, Edinburgh, UK.
 Kolesik, M. (2014, 08/17/14  08/21/14). Invited: Ultrafast lightmatter coupling in condensed and gaseous nonlinear media. SPIE Optics + Photonics 2014. SanDiego, U.S.A..
 Kolesik, M. (2014, 09/18/14 — 09/24/14). Invited: NonPerturbative LightMatter Interaction Model for Optical Filamentation. COFIL 2014, 5th International Symposium on Filamentation. Shanghai, China.
 Kolesik, M. (2014, 10/19/14  10/23/14). Invited: Nonperturbative timedomain modeling of lightmatter interactions for computer simulation in extreme nonlinear optics. Frontiers in Optics/Laser Science Conference. Tucson, U.S.A..