La propagación de haces de luz en cristales fotovoltaicos piroeléctricos fotorefractivos se modela mediante una generalización específica de la ecuación de Schrödinger no lineal (GNLSE). Utilizamos una aproximación variacional (VA) para predecir la propagación de entradas de ondas solitarias en los cristales, encontrando que las ecuaciones de VA involucran una función especial de dilogaritmo. La VA predice que existen solitones y respiradores, y el criterio de Vakhitov-Kolokolov predice que los solitones son soluciones estables. Las simulaciones directas de la GNLSE subyacente corroboran la existencia de tales modos estables. Las soluciones numéricas producen tanto respiradores regulares como aquellos que presentan pulsaciones (modulaciones de largo período de oscilaciones rápidas). En este último caso, la transformada de Fourier de las oscilaciones de amplitud revela un espectro casi discreto que caracteriza la dinámica de las pulsaciones. Las soluciones numéricas de otro tipo demuestran la división espontánea del pulso de entrada en dos o varios pulsos secundarios.